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EXECUTIVE  SUMMARY 


One  of  the  basic  problems  to  be  solved  by  geotechnical  engineers  in 
regions  where  earthquake  hazards  exist  is  to  estimate  the  site -specific 
dynamic  response  of  a  layered  soil  deposit  for  level-ground  conditions.  The 
computer  program  described  and  provided  in  this  report,  WESHAKE,  may  be  used 
to  accomplish  this  task.  WESHAKE  is  an  adaptation  of  the  original  computer 
program,  SHAKE,  written  at  the  University  of  California  at  Berkeley  by 
Schnabel,  Lysmer,  and  Seed  (1972).  WESHAKE  was  created  and  has  been 
continually  modified  by  WES  to  keep  pace  with  state-of-the-art  technology  and 
provide  a  user-friendly  interface. 

The  WESHAKE  package  consists  of  this  report  and  a  single  floppy  disk 
that  contains  the  executable  program,  data  base  files,  a  plotting  program,  and 
example  input  and  output  files.  It  is  imperative  that  the  user  of  this 
program  have  a  competent  understanding  of  the  problem  statement,  basic 
assumptions ,  and  mathematical  formulation  used  by  the  original  authors  of 
SHAKE.  Documentation  of  the  original  program,  SHAKE,  is  not  duplicated  herein 
nor  is  this  report  intended  to  be  a  primer  on  dynamic  site  response  analysis. 
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PREFACE 


This  study  Is  sponsored  by  the  Headquarters,  US  Army  Corps  of  Engineers 
(USACE)  under  the  Numerical  Model  Maintenance  Program  (NMMP) .  This  program 
provides  for  the  maintenance,  documentation,  and  corrections  of  existing 
computer  models  (programs)  that  have  existing,  or  the  potential  for, 
widespread  usage  among  Corps  personnel.  The  program  also  provides  for  user 
consultation  with  WES  authors.  Mr.  Richard  Davidson,  USACE,  is  the  Technical 
Monitor  for  this  particular  model. 

This  report  and  accompanying  software,  the  first  in  the  report  series  on 
Geotechnical  Earthquake  Engineering  Software  (GEES),  contains  the  information 
necessary  to  run  the  program  WESHAKE.  WESHAKE  is  a  wave  equation  solver  for 
one -dimensional  problems  and  has  been  widely  used  to  solve  problems  in 
earthquake  engineering  for  USACE  projects. 

The  purpose  of  establishing  GEES  is  to  provide  a  set  of  easy  to  use  and 
understandable  tools  that  can  support  the  needs  of  the  district  and  division 
engineers  in  evaluating  the  dynamic  response  effects  of  earthquakes  on 
foundations,  earth  structures,  and  soil-structure  systems  at  specific  sites 
throughout  the  world.  GEES  will  also  allow  USACE  to  establish  and  maintain 
consistency  of  programs  among  offices  and  a  center  for  validation  studies. 

The  WES  Principal  Investigator  was  Mr.  David  W.  Sykora,  Earthquake 
Engineering  and  Seismology  Branch  (EESB) ,  Earthquake  Engineering  and 
Geosciences  Division  (EEGD) ,  Geotechnical  Laboratory  (GL) ,  WES.  Mr.  Ronald  E. 
Wahl,  Soil  and  Rock  Mechanics  Division  (SRMD),  GL,  initiated  the  study  and 
provided  technical  assistance.  Mr.  Michael  K.  Sharp,  Engineering  Geophysics 
Branch  (EGB) ,  EEGD,  GL,  has  provided  useful  additions  and  suggestions  to  ^he 
core  computer  program  over  the  past  four  years.  Dr.  David  C.  Wallace, 

Illinois  State  University,  performed  most  of  the  programming  for  thi  .  study 
while  at  WES  during  the  summers  of  1991  and  1992  under  the  US  Army  Summer 
Faculty  Research  and  Engineering  Program  (SFREP)  provided  through  the  US  Army 
Research  Office.  Messrs.  Willie  McGeehee  and  Daniel  Habeeb  and  Ms.  Jennifer 
Davis,  EESB,  drafted  figures,  made  copies,  and  helped  to  prepare  the  final 
report.  Dr.  Mary  Ellen  Hynes  was  Chief,  EESB,  during  the  course  of  this  study 
and  provided  direct  technical  oversight. 

Overall  direction  at  WES  was  provided  by  Dr.  A.  G.  Franklin,  Chief, 

EEGD,  and  Dr.  William  F.  Marcuson  III,  Chief,  GL.  Ms.  Mary  K.  Vincent, 
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Chief,  Office  of  Technical  Programs  and  Plans,  was  the  overall  WES  program 
manager  of  the  NMHP. 

At  the  time  of  publication  of  this  report.  Director  of  WES  was  Dr. 
Robert  W.  Whalin.  Commander  and  Deputy  Director  was  COL  Leonard  G.  Hassel, 
EN. 
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CONVERSION  FACTORS,  NON-SI  to  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non- SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows; 


Multinlv 

Abbreviation 

Bv 

To  Obtain 

feet  per  second 

fps 

0.3048 

metres  per  second 

kips  per  cubic 

kef 

16.03 

megagrams  per  cubic 

foot 

metre 
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USAGE  GEOTECHNICAL  EARTHQUAKE  ENGINEERING  SOFTWARE 
WESHAKE  FOR  PERSONAL  COMPUTERS 


PART  I :  INTRODUCTION 

Background 

1.  One  of  the  basic  problems  to  be  solved  by  geotechnical  engineers  in 
regions  where  earthquake  hazards  exist  is  to  estimate  the  site-specific 
dynamic  response  of  a  layered  soil  deposit  under  a  level  ground  surface.  This 
problem  is  commonly  referred  to  as  a  site-specific  response  analysis  or  soil 
amplification  study  (although  motions  may  be  deamplified) .  The  solution  of 
this  problem  allows  the  geotechnical  engineer  to  evaluate  the  potential  for 
liquefaction,  to  conduct  the  first  analytical  phase  of  seismic  stability 
evaluations  for  slopes  and  embankments,  to  calculate  site  natural  periods,  to 
assess  ground  motion  amplification,  and  to  provide  structural  engineers  with 
various  parameters,  primarily  response  spectra,  for  design  and  safety 
evaluations  of  structures. 

2.  A  site-specific  response  analysis  is  critical  to  US  Army  projects, 

both  Civil  Works  and  Military.  Department  of  the  Army,  Engineering  Regulation 

1110-2-1906  (1983)  provides  guidance  for  US  Army  Corps  of  Engineers  (USAGE) 

Civil  Works  projects.  For  "Embankments  and  Soil  Foundations"  projects  located 

within  seismic  zones  2,  3,  and  h  (refer  to  Figures  I-l  through  1-4),  the 

pseudo-static  method  of  analysis  is  superseded.  Rather: 

...appropriate  analytical  techniques  [shall  be  used]  to  evaluate 
liquefaction  potential  and/or  to  estimate  deformations,  beginning 
with  the  more  simplified  methods,  and  progressing  as  necessary  to 
more  rigorous,  sophisticated  procedures. 

US  Army  Technical  Manual  TM  5-809-10-1  (1986)  provides  seismic  design 
guidelines  for  essential  buildings  (Military  projects).  The  "analytical  soil- 
column  response"  method  represents  one  of  three  methods  that  can  be  used  to 
develop  site-specific  response  spectra  (refer  to  section  3-6  and  Appendix  C, 
section  C-3  of  TM  5-809-10-1). 

3.  The  computer  program  SHAKE  was  written  in  the  early  1970 's  by 
Schnabel,  Lysmer,  and  Seed  (1972)  to  conduct  analytical  site  response  analyses 
via  solution  of  the  wave  equation.  This  program  has  been  distributed  freely 
and  is  still  widely  used  by  the  profession  although  many  versions  of  the 
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program  have  been  modified  by  different  organizations  (e.g.,  GeoTech 
International,  Ltd.  1985).  This  program  has  been,  and  continues  to  be, 
successfully  validated  with  measured  earthquake  motions  and  site  response. 

The  US  Army  Waterways  Experiment  Station  (WES)  has  been  using  the  computer 
program  SHAKE  to  calculate  site  response  for  level-ground  soil  sites  for  more 
than  15  years,  including  use  on  a  number  of  USAGE  projects.  A  partial  list  of 
projects  are  provided  in  a  supplemental  Bibliography  at  the  end  of  this 
report . 

4.  WES  has  continually  made  adaptations  to  SHAKE  as  the  use  f  r  each 
new  project  required.  The  original  version  for  use  on  a  pe  jonal  computer  was 
obtained  from  the  University  of  California  at  Berkeley  (UCB)  around  1985. 

This  program  at  WES  is  now  called  WESHAKE  to  reflect  the  numerous  changes  that 
have  been  made  to  keep  pace  with  state-of-the-art  technology,  to  provide  for 
needs  of  USAGE  users,  and  to  provide  a  user-friendly  interface.  These 
adaptations  facilitate  transfer  technology  to,  and  wide-spread  use  among, 

USAGE  personnel. 

Purpose 

5.  The  purpose  of  this  report  and  software  is  to  provide  USAGE  district 
and  division  engineers  a  means  to  calculate  the  horizontal  site  response  of 
level-ground  soil  sites  caused  by  vertically-propagating,  horizontally- 
polarized  shear  waves  that  can  be  used  to  solve  a  variety  of  site  response 
problems .  Data  bases  have  also  been  created  to  attempt  to  reduce  the  amount 
of  effort  in  preparing  cr  input  file  to  solve  a  particular  problem.  It  is  the 
intention  of  the  authors  to  create  a  user  interface  that  is  convenient  to  use 
and  requires  little,  if  any,  guidance  from  the  user's  manual. 

6.  WESHAKE  is  part  of  the  Geotechnical  Earthquake  Engineering  Software 
(GEES)  library  that  was  est, '-lished  to  prepare,  validate,  and  raair*-'>in 
programs  used  to  evaluate  the  dynamic  response  effects  of  earthquv.  _  on 
foundations,  earth  structures,  and  soil-structure  systems  and  establish  a 
center  for  free  distribution  and  support.  Three  programs  in  the  GEES  library 
are  currently  being  supportc  ,  all  by  the  Numerical  Model  Maintenance  Program 
(NMMP);  WESHAKE,  WESRISK,  and  WESFUJSH . 
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Intended  Users 


7.  It  is  imperative  that  the  user  of  WESHAKE  have  a  competent 
understanding  of  the  problem  statement,  basic  assumptions,  and  mathematical 
formulation  used  by  the  authors  of  SHAKE.  The  simplicity  of  the  user 
interface  should  not  be  associated  with  the  minimum  level  of  capability  of  the 
user.  The  interface  is  designed  to  allow  the  first- time  or  occasional  user 
with  a  rapid  and  easy  means  to  run  the  program.  Documentation  of  the  original 
program,  SHAKE,  is  not  duplicated  herein. 

Merits  of  WESHAKE 

8.  WESHAKE  is  considered  to  be  an  improvement  over  SHAKE  for  all  users 
for  the  following  reasons: 

a .  Shear  modulus  can  be  defined  by  shear  wave  velocity  or  K2 
individually  for  each  layer. 

b .  Recent  sets  of  shear  modulus  degradation  and  damping  curves  are 
available  in  a  data  base . 

c .  The  relationships  for  normalized  shear  moduli  and  damping  ratio  can 
be  specified  separately  for  each  material. 

d.  The  experience  of  WES  engineers  has  been  interjected  into  the  data 
input  requirements  and  in  the  discussions  of  the  report. 

e .  An  option  for  interactive  plotting  of  earthquake  motions  is 
included. 

WESHAKE  should  benefit  the  first- time  or  occasional  user  for  the  following 
reasons : 

a .  The  menus  and  displays  allow  the  user  to  create  an  input  file  and 
conduct  the  analysis  without  need  to  reference  the  user's  manual  or  have 
knowledge  of  the  format  of  input  fields. 

b .  Some  parameters  that  take  on  typical  values  have  been  assximed  for 
the  initial  runs. 

c .  Only  parameters  required  for  analysis  are  part  of  the  query 
sequences . 

d.  Several  error  checks  have  been  incorporated  including  the 
establishment  of  reasonable  bounds  for  some  data  values. 

9.  Despite  the  many  cosmetic  changes  made  as  part  of  the  study,  the 
core  structure  of  subroutines,  algorithms,  and  data  flow  have  not  changed 
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significantly.  Only  minor  differences  exist  between  formats  used  for  SHAKE 
and  WESHAKE.  These  differences  are  documented  in  this  report. 

Report  Organization  and  Suggested  Use 

10.  This  report  is  intended  to  serve  as  a  guide  to  the  first- time  or 
occasional  user  of  WESHAKE.  Information  for  the  data  files  are  collected  in 
an  interactive  mode  with  the  computer  providing  menus  and  information 
requests.  Example  input  screens  and  responses  are  provided  in  this  manual  to 
give  the  user  a  clearer  understanding  of  the  basic  requirements.  Output 
screens,  files,  and  reports  are  also  provided  to  illustrate  the  relationship 
between  the  inputs  and  the  outputs  as  well  as  the  interpretation  of  the 
results.  Instructions  for  installation,  hardware  requirements,  listings  of 
data  bases  and  example  files ,  and  file  formats  are  provided  in  appendices  to 
this  report. 

11.  Some  basic  background  information  about  the  analytical  solution  to 
site  response  analysis  and  aspects  of  WESHAKE  that  are  different  from  SHAKE 
are  presented  in  Part  II  of  this  report.  The  salient  features  of  WESHAKE  and 
validation  are  presented  in  Part  III  of  this  report.  The  use  and  execution  of 
WESHAKE  are  described  in  Part  IV  of  this  report.  The  first  section  of  Part  IV 
focuses  on  the  preprocessing  stage  of  collecting  information.  The 
preprocessing  stage  is  divided  into  two  primary  sections:  mandatory  actions 
and  user  options.  The  mandatory  actions  are  used  to  define  the  soil  column, 
identify  an  earthquake  motion,  and  assign  the  location  of  object  motion.  User 
options  allow  for  variations  in  input  parameters  and  the  calculation  and 
preparation  of  specialized  forms  of  output.  Other  sections  of  Part  IV 
describe  the  analysis  stage  (execution)  of  WESHAKE  and  error  checks.  An 
example  problem  is  used  throughout  this  report  to  assist  the  user. 

12.  Part  V  of  this  report  presents  discussion  on  how  to  model  and  solve 
more  difficult  problems  that  may  not  have  been  addressed  previously.  It  is 
considered  to  be  a  supplement  to  the  other  sections  and  is  not  a  necessary 
part  of  learning  how  to  use  WESHAKE  to  solve  basic  problems.  The  reader 
should  be  familiar  with  the  topics  discussed  in  Part  V  and  refer  to  the 
appropriate  discussions  as  needed. 

13 .  The  complete  package  for  WESHAKE  includes  this  report  and  a  floppy 
disk  with  the  following  files: 
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VESHAKE.EXE 
SHEARDB 
DAMFDB 
EARTHQ 
MOTION.EXE 
EXAMPLE.DAT 
EXAMPLE.  EXT 
EXAMPLE. SPF 
READ  ME. TXT 


Executable  shell  containing  wave  propagation  code 

Data  base  of  shear  modulus  relationships 

Data  base  of  damping  relationships 

Data  base  of  accelerograms 

Executable  accelerogram  plotting  program 

Example  specification  file 

Companion  to  specification  file 

Soil  profile  file 

List  of  files  and  file  sizes 


Questions,  comments,  and  requests  for  updates  should  be  directed  to: 


U.S.  Army  Engineer  Waterways  Experiment  Station 
ATTN:  CEWES-GG-H  (Mr.  David  W.  Sykora) 

3909  Halls  Ferry  Road 
Vicksburg,  MS  39180-6199 


Voice:  (601)  634-3551 
FAX:  (601)  634-3453 
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Figure  I-l,  Seismic  zones  within  contiguous  48  United  States  (Department  of 
the  Army  1983;  adapted  from  Algermissen  et  al.  1982) 


Figure  1-2.  Seismic  zones  within  California  (Department  of  the 
Army  1983;  adapted  from  Algermissen  et  al.  1982) 
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Figure  1-3.  Seismic  zones  within  Alaska  (Department  of  the  Army  1983; 
adapted  from  Algermlssen  et  al.  1982) 
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Figure  1-4.  Seismic  zones  within  Hawaii  (Department  of  the  Army  1983; 
adapted  from  Algermlssen  et  al.  1982) 


PART  II:  WESHAKE  BACKGROUND  AND  THEORY 

Procedure  of  Site  Response  Analysis 

14.  A  site  response  analysis,  sometimes  referred  to  as  a  soli 
amplification  analysis,  Involves  the  determination  of  components  of  ground 
motion  for  design  or  seismic  evaluation.  Typically,  as  in  this  study,  that 
determination  is  made  for  a  "free-field"  response  --  the  response  at  the 
ground  surface  of  an  Ideal  soil  deposit  (horizontal  layers  extending  to 
infinity)  to  a  spatially-uniform  motion  applied  at  the  base.  The  conceptual 
relationship  between  free-field  response  with  respect  to  two  primary  control 
points  --  rock  outcrop  and  base  rock  --  in  a  site  response  analysis  is  shown 
in  Figure  II -1.  The  motions  at  these  three  points,  as  well  as  any  other  point 
in  the  vertical  profile,  are  unique.  Design  earthquakes  are  frequently 
specified  as  corresponding  to  a  rock  outcrop.  Mathematical  expressions 
(transfer  functions)  are  then  used  to  find  the  equivalent  motion  for  the  base 
rock  and  then  the  seismic  waves  are  propagated  through  the  soil  colvimn  to 
determine  the  free-field  motion. 

15.  The  determination  of  site-specific  earthquake  response  of  soil 
deposits,  then,  generally  involves  three  basic  steps: 

a.  Selection  of  earthquake  motions,  usually  corresponding  to  rock 
outcrop . 

b.  Idealization  of  stratigraphy  and  selection  of  material 
properties . 

c.  Calculation  and  evaluation  of  site  response. 

The  third,  and  final,  step  of  a  site-specific  earthquake  response  analysis  is 
the  subject  of  this  study. 

16.  Different  techniques  are  available  to  determine  site-specific 
response  of  soil  sites  to  earthquake  motions.  Analytical  formulations  include 
the  wave  equation  and  shear  beam  analogies  (both  continuous  formulations)  and 
a  lumped  mass  model  analysis  (discrete  formulation).  Initial  formulations  for 
site  specific  calculations  using  the  wave  equation  were  reported  in  the  U.S. 
by  Roesset  and  Whitman  (1969)  and  Roesset  (1970)  and  have  been  enhanced  since. 
A  number  of  computer  codes  are  available  to  solve  the  wave  equation  in  one, 
two,  or  three  dimensions.  This  report  summarizes  the  code  WESHAKE  which 
evolved  from  the  code  SHAKE,  described  in  the  next  section. 
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SHAKE 


17.  SHAKE  was  developed  at  the  University  of  California  at  Berkeley 
(Schnabel,  Lysmer,  and  Seed  1972)  and  written  in  FORTRAN  IV  to  run  on  a  CDC 
6400  computer.  It  has  since  been  adapted  to  run  on  a  number  of  computer 
platforms  including  personal  computers  by  various  sources.  SHAKE  is  widely 
used  by  the  geotechnical  earthquake  engineering  profession  for  the  calculation 
of  site  response  for  horizontal  motions. 

18 .  Several  investigators  have  reported  close  comparisons  between  the 
results  using  SHAKE  and  the  measured  horizontal  response  from  strong-motion 
instruments  triggered  during  earthquakes  at  site  periods  less  than  2  sec 
(e.g.,  Seed  et  al.  1987,  Idriss  1990,  and  Seed,  Dickenson,  and  Idriss  1991). 
The  experience  of  these  investigators  suggest  that  for  calculated  site  periods 
greater  than  4  sec,  motions  are  likely  to  be  significantly  affected  by  two- 
dimensional  effects  and  surface  wave  energy  and  are  not  well  represented  with 
SHAKE.*  The  user  should  be  cautious  when  using  values  at  periods  greater 
than  2  sec. 

19.  SHAKE  was  developed  to  calculate  the  horizontal  response  caused  by 
an  earthquake  at  any  depth  of  a  soil  profile.  The  approach  and  algorithms 
incorporated  in  the  program  are  simple,  straight  forward,  and  adequate  for  the 
purpose  intended  as  clearly  evident  through  the  prolific  publication  of 
results  and  favorable  comparisons  with  measured  response  (e.g..  Seed  et  al. 
1987,  Idriss  1990,  and  Seed,  Dickenson,  and  Idriss  1991).  The  simplicity 
associated  with  SHAKE  is  attributed  to  some  basic  assumptions  regarding  the 
cyclic  behavior  of  materials  and  geometry  of  the  problem.  The  basic 
assumptions  used  in  the  formulation  of  are: 

a.  The  soil  layers  are  horizontal  and  extend  to  infinity. 

b.  The  ground  surface  is  level. 

£.  Each  soil  layer  is  completely  defined  by  the  shear  modulus  and 
damping  as  a  function  of  strain,  the  thickness,  and  unit  weight. 

d.  The  non-linear  cyclic  material  behavior  is  adequately  represented 
by  the  linear  visco-elastic  (Voigt)  constitutive  model  and 
implemented  with  the  equivalent- linear  method. 

e.  The  incident  earthquake  motions  are  spatially-uniform, 
horizontally-polarized  shear  waves,  and  propagate  vertically. 


*  Personal  communication,  Prof.  Raymond  Seed,  University  of  California  at 
Berkeley,  23  September  1991. 
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In  general,  assumptions  (a),  (b) ,  and  (c)  used  to  derive  this  model  would  seem 
to  significantly  limit  the  applicability  of  this  method.  Past  studies  have 
shown,  however,  that  reasonable  results  are  obtained  for  a  much  broader 
spectrum  of  in  situ  conditions.  The  equivalent- linear  constitutive  model, 
assumption  (d) ,  is  described  later  in  this  section.  The  last  assumption,  (e) 
above,  narrows  the  focus  to  a  simple  class  of  problems,  but,  is  a  common 
assumption  for  this  type  of  problem. 

20.  It  is  important  to  realize  that  the  formulation  of  SHAKE  for  wave 
propagation  is  based  on  a  total  stress  analysis.  The  materials  are  considered 
to  be  continue  and  pore  water  pressures  are  non-existent.  The  calculation  of 
shear  modulus  using  values  of  Kj  does  involve  the  determination  of  mean 
effective  stress  using  the  depth  of  the  water  table  and  the  unit  weight  of 
water . 

Formulation  and  iteration  scheme 

21.  The  one -dimensional  wave  equation  model  (Kanai  1951)  was  used  to 
develop  SHAKE.  This  model  has  proven  to  be  effective  despite  the  simplicity 
and  number  of  assumptions  involved.  The  solution  algorithm  involves  the 
complex  response  technique  and  the  Fast  Fourier  Transform  (Cooley  and  Tukey 
1965) .  The  general  formulation  of  the  wave  equation  is  not  unique  to 
horizontally-polarized  shear  wave  motion;  the  equation  can  also  be  solved  for 
the  vertical  propagation  of  compression  waves. 

22.  In  general,  soil  is  a  non-linear  material  that  exhibits  hysteretic 
behavior  under  cyclic  loading.  An  example  of  the  stress -strain  behavior  is 
shown  in  Figure  II-2a.  Soil  is  difficult  to  model  accurately  for  cyclic 
response;  exact  representations  are  unavailable.  The  constitutive  model 
incorporated  into  SHAKE  is  linear  with  simulated  nonlinear  effects  to  account 
for  dependency  of  moduli  on  shear  strain.  The  method  used  to  implement  the 
linear  visco-elastic  model,  called  the  equivalent -linear  method,  was  proposed 
by  Seed  and  Idriss  (1970)  and  is  widely  used  in  geotechnical  earthquake 
engineering  studies. 

23.  The  basic  components  of  the  equivalent -linear  method  are  the 

maximum  shear  modulus,  ,  moist  unit  weight,  y  .  and  ratio  of  critical 

damping,  P  .  The  property  ,  which  corresponds  to  the  linear-elastic, 

continuum  material  property  (Lam4  1852),  can  be  calculated  from  low-strain 
seismic  shear  wave  velocity  using; 
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(1) 


G-x 


where 

g  -  gravitational  constant  of  earth 
V,  -  shear  wave  velocity 

When  using  the  shear  modulus  coefficient  in  lieu  of  V,,  the  following  equation 
is  used  (Seed  and  Idriss  1970): 


=  1000  (K,)^  (o')®  *  (2) 

where 

o„'  -  mean  effective  stress,  in  psf 
G„««  is  in  psf 

The  shear  modulus,  G,  of  a  soil  remains  constant  during  cyclic  loading  at  very 
small  shear  strains  (defined  as  G„.„  ) .  As  the  shear  strains  increase  above 
some  threshold  value,  generally  accepted  to  be  about  10**  percent  or  less,  G 
decreases.  The  equivalent- linear  method  uses  secant  shear  moduli  that  are 
adjusted  during  each  iteration  to  account  for  this  degradation  of  shear 
modulus.  Damping  is  input  by  using  complex  moduli,  G*  ,  and  hysteretic 
damping  (which  is  independent  of  frequency)  as  reported  by  Udaka  and  Lysmer 
(1973): 


G*  -  G  (l-2P*+2iPv'r^) 


(3) 


where 

i  =  ^-1 

Damping  increases  as  shear  strain  increases. 

24.  The  character  of  the  relationships  between  normalized  shear  modulus 
versus  shear  strain  and  damping  ratio  versus  shear  strain  was  addressed  in 
studies  at  the  University  of  Kentucky  in  the  late  1960 's  (Hardin  and  Drnevich 
1972a;  1972b).  Seed  and  Idriss  (1970)  and  Schnabel  (1973)  used  the  results  of 
these  studies  to  derive  the  equivalent  linear  model  and  the  first  set  of 
relationships  provided  with  SHAKE.  The  general  shapes  of  these  relationships 
are  shown  in  Figure  II-3.  Relationships  provided  with  WESHAKE  are  presented 
in  Part  III  of  this  report. 

25.  An  example  of  the  iterative  procedure  for  the  equivalent- linear 
method  is  shown  in  Figure  II-2b  and  described  below.  Assuming  shear  wave 
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propagation,  the  model  is  initiated  with  an  arbitrary  value  of  shear  modulus, 

,  chosen  to  be  less  than,  or  equal  to,  G _ _  For  the  first  cycle  of 

loading,  the  stress -strain  relation  is  linear  between  the  two  levels  of  shear 
strain,  +  ,  with  a  slope  of  Gj^  .  The  ordered  pair  (G^  ,  )  comes  from 

the  appropriate  modulus  degradation  curve  as  discussed  in  Part  III  of  this 
report  and  shown  schematically  in  Figure  II -2b.  Maximum  shear  strains  are 
obtained  from  the  solution  of  the  wave  equation.  The  ratio  of  effective  shear 
strain  to  maximum  shear  strain,  PRMUL,  (assumed  to  be  0.65)  is  used  to  obtain 
a  new  value  of  shear  modulus ,  G2  ,  from  the  appropriate  modulus  curve .  A  new 
value  of  P  is  also  obtained. 

26.  This  process  is  repeated  until  the  difference  in  moduli  and  damping 
for  two  successive  iterations  are  within  a  prescribed  tolerance,  ERR  (5 
percent  is  assumed) .  The  number  of  iterations  required  by  the  computer  to 
solve  the  problem  is  a  function  of  how  close  the  initial  estimates  are  to  the 
final  results;  the  closer  the  two  sets  of  values,  the  fewer  the  number  of 
iterations  required  (and  proportionally,  more  time  to  solve). 
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Figure  Il-l.  Three  primary  control  points  for  site  response  analysis 
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equivalent- linear  soil  model  for  a  constant  stress  state 
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Standard  relationships  between  normalized  shear  modulus  and 
damping  ratio  versus  shear  strain 


PART  III:  PERSONAL  COMPUTER  IMPLEMENTATION 


27.  The  computer  program  WESHAKE  is  primarily  a  menu  driven  shell 
system  with  various  other  requests  for  information  from  the  user  to  quickly 
and  conveniently  guide  the  user  to  execution  of  the  core  computer  code.  This 
program  is  linked  with  overlays  and  external  data  bases  to  reduce  memory  size 
during  execution.  Information  about  loading  the  program  and  companion  files 
is  presented  in  Appendix  A  along  with  discussions  about  hardware  requirements, 
array  limitations,  memory  requirements,  and  program  run  times. 

Changes  from  SHAKE 

28.  Over  the  20  years  that  SHAKE  has  been  used,  more  knowledge  has 
become  available  with  regard  to  the  specification  of  inputs  to  the  program  and 
significant  advances  have  been  made  in  computer  technology.  As  these  findings 
and  advances  have  been  made,  WES  has  refined  and  adapted  SHAKE.  WESHAKE  can 
best  be  described  as  a  program  that  performs  pre-  and  post-processing 
functions  and  uses  the  original  program,  called  SHAKEl,  as  a  core  for 
calculations.  SHAKE,  then,  is  essentially  embedded  in  WESHAKE  which 
facilitates  the  implementation  of  newer  versions  or  formulations  of  SHAKE. 

29.  The  pre-  and  post-processing  routines  were  written  from  scratch 
during  FY91  and  version  1.0  was  finalized  and  distributed  in  FY92.  Only  a  few 
minor  changes  were  made  to  SHAKEl  to  reflect  changes  in  the  state-of-the-art. 
The  changes  made  to  SHAKEl  were  purposely  kept  to  a  minimum  to  retain  the 
intent  of  that  program's  authors  and  minimize  the  amount  of  re-validation 
required.  The  options  to  SHAKE  are  also  numbered  differently  in  the  interface 
to  reduce  the  size  of  user  menus.  The  cross  reference  between  option  numbers 
is  provided  in  Table  III-l. 

30.  The  most  significant  change  made  to  SHAKEl  was  in  the  subroutine 

SOILIN  which  is  used  to  read,  calculate,  and  store  the  initial  values  of  shear 
wave  velocity  and  shear  modulus.  First,  the  initial  value  of  shear  modulus 
for  each  layer  used  in  the  iteration  scheme  for  the  equivalent  linear 
constitutive  model  is  fixed  to  equal  Fixing  the  initial  value 

corresponding  to  will  Increase  the  time  to  iterate  to  the  solution  but 
not  significantly  when  using  personal  computers  with  processor  speeds  greater 
than  about  20  MHz.  Then,  rather  than  having  fixed  variable  inputs  of  and 
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Kj  for  clay  and  sand,  respectively,  SHAKEl  allows  the  user  to  input  Kj  OR  V, 
for  any  soil;  is  no  longer  used.  V,  input  is  still  required  for  rock, 
however . 

31.  Other  significant  changes  include  the  separation  of  shear  modulus 
and  damping  relationships  for  soil  types.  Previously,  a  set  of  shear  modulus 
and  damping  relationships  was  selected.  Now,  each  relationship  can  be  chosen 
separately.  This  is  described  in  more  detail  in  the  next  section.  WESHAKE 
has  no  provisions  for  the  use  of  sublayers.  This  original  option  does  not 
appear  to  be  useful.  In  addition,  error  checks  were  added  and  will  be 
described  later. 


Shear  Modulus  and  Damping  Ratio  Data  Bases 

32 .  Two  separate  data  bases  are  used  to  contain  options  for  specifying 
the  variation  of  normalized  shear  modulus  and  damping  ratio  with  effective 
shear  strain  --  SHEAROB  and  DAMPDB,  respectively.  Both  data  bases  are  used 
independently  to  specify  material  relationships  for  each  layer  allowing 
versatility.  (In  previous  versions  of  SHAKE  and  WESHAKE  the  selection  of  one 
material  relationship  required  the  use  of  a  specific  damping  ratio 
relationship.)  A  description  of  the  data  bases,  including  format,  are 
presented  in  Appendix  B. 

33.  Eleven  relationships  for  normalized  shear  modulus  and  damping  ratio 
exist  in  the  two  data  bases  --  one  for  rock,  one  for  gravel,  three  for  sand, 
and  six  relations  for  cohesive  soils.  These  relationships  are  shown  in 
Figures  III-l  and  III-2  for  normalized  shear  modulus  and  Figures  III-3  and 
III -4  for  damping  ratio.  The  sources  of  the  relationships  are  summarized  in 
Table  III -2.  The  criteria  for  selection  of  appropriate  relationships  from 
the  data  bases  are  based  on  the  general  soil  classification  and  plasticity 
index  (for  fine-grained  soils). 

34.  Three  different  relationships  are  available  in  each  of  the  data 
bases  for  sands  --  average,  lower  bound,  and  upper  bound  (refer  to  Figures 
III-l  and  III-3).  These  relationships  are  based  on  the  results  of  laboratory 
studies  (e.g.,  Hardin  and  Drnevich  1972a;  Seed  and  Idriss  1970,  Wong  et  al. 
1974;  Seed  et  al.  1986;  and  Hynes  1988)  Some  investigators  (e.g.  Idriss  1990) 
tend  to  use  the  upper  bound  relationship  for  sands .  The  general 
recommendation  of  WES  is  to  use  the  lower  bound  sand  relationship  for  gravelly 


26 


sands  or  sandy  gravel  (rather  than  the  "gravel"  curve) ,  use  the  average 
relationship  for  clean  sands ,  and  use  the  upper  bound  relationship  for  clayey 
and  silty  sands. 

35.  The  results  of  independent  laboratory  studies  by  Sun,  Golesorkhi, 
and  Seed  (1988),  Zen  and  Higuchi  (1984),  and  Vucetic  and  Dobry  (1991)  have 
shown  that  for  cohesive  soils,  the  selection  on  an  appropriate  shear  modulus 
curve  should  be  based  on  the  plasticity  index.  The  set  of  relationships 
proposed  by  Sun,  Golesorkhi,  &  Seed  (1988)  are  included  in  the  data  base.  A 
discussion  related  to  the  choice  of  this  set  is  contained  in  Part  V  of  this 
report.  Note  that  the  range  of  relationships  for  cohesive  soils  overlaps  the 
range  of  relationships  for  sands . 

Earthquake  Record  Data  Base 

36.  Lists  of  earthquake  records  (accelerograms)  in  the  WESHAKE  data 
base  are  separated  into  measured  records  and  synthetic  records  and  are 
presented  in  Tables  III-3  and  III-4,  respectively.  More  information  on  the 
data  base  is  provided  in  Appendix  C,  including  data  base  sjmtax  and  plots  of 
the  variations  of  acceleration  with  time  and  velocity  spectra  for  all  of  the 
records . 

37.  The  data  base  used  in  WESHAKE  has  been  limited  initially  to  22 
records  representing  both  measured  and  synthetic  earthquake  records.  Most 
importantly,  this  data  base  purposely  Includes  only  records  corresponding  to 
rock  outcrop  (refer  to  Figure  II -1).  The  user  may  define  an  earthquake  not 
contained  in  the  data  base,  preferably  a  project- specific  record. 

38 .  A  number  of  (measured)  earthquake  records  are  available  to  the 
user,  particularly  with  the  establishment  of  national  data  bases  for  strong 
motion  (e.g..  Row  1990;  Friberg  and  Jacob  1990).  The  user  should  be  aware 
that  the  vast  majority  of  earthquake  records  have  been  measured  at  points 
other  than  rock  outcrop,  the  most  desired  point  for  WESHAKE.  Most  of  the  U.S. 
records,  for  instance,  have  been  measured  by  instruments  housed  in  buildings, 
primarily  the  basements,  and  at  the  top  of  free-field  soil  sites  (well  above 
the  top  of  rock) .  The  presence  of  the  building  and  soil  overlying  bedrock  can 
greatly  affect  the  recorded  motion  (e.g.,  refer  to  the  presentation  under  the 
previous  sub-heading).  In  fact,  this  is  the  purpose  of  conducting  an  analytic 
site  response  study.  Free-field  rock  outcrop  records  (or  base  rock  records. 
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when  they  become  available)  should  be  used  to  properly  evaluate  site  response 
and  maintain  consistency  in  approach. 

39 .  The  variation  of  acceleration  with  time  for  each  earthquake  can  be 
viewed  on  the  computer  monitor  by  enabling  the  computer  program  MOTION 
(outside  of  the  WESHAKE  shell) .  The  plots  represent  unscuLed  (measured 
accelerations)  records  as  they  are  contained  in  the  database.  Details 
concerning  its  procedures  and  operations  are  included  in  Appendix  D. 

Input  and  Output  Files 

40.  WESHAKE  creates  and  manages  a  number  of  different  files  as  needed. 
The  use  of  input  files  may  go  relatively  unnoticed  to  the  user.  Output  files 
are  created  to  be  easily  read  by  commercial  plotting  software.  Both  input  and 
output  files  are  described  below. 

Input  files 

41 .  WESHAKE  makes  use  of  three  different  files  for  input  to  the  core 
program  SHAKEl:  a  specification  file  (with  *.DAT  filename  extension),  a  soil 
profile  file  (with  *.SPF  filename  extension),  and  a  companion  file  (with  *.EXT 
filename  extension).  The  specification  file  is  the  complete  input  file  to 
SHAKEl.  A  soil  profile  file  contains  information  regarding  the  stratigraphy 
and  material  properties  and  is  a  subset  of  the  specification  file.  Its  use  is 
optional  but  useful  when  a  number  of  specification  files  are  to  be  created 
with  the  same  soil  profile.  For  instance,  an  engineer  may  be  interested  in 
subjecting  a  particular  soil  column  to  a  number  of  different  earthquake 

records  scaled  to  the  same  value  of  a _  .  The  companion  file  is  merely  a 

storage  area  for  descriptions  of  the  soil  types  and  is  created  by  WESHAKE 
along  with  the  specification  or  soil  profile  file.  The  manipulation  of  these 
three  files  is  relatively  transparent  to  the  user.  Options  within  WESHAKE 
allow  the  user  to  list  existing  files  in  the  working  directory  to  avoid 
confusion. 

42.  Old  input  files  (from  SHAKE  or  WESHAKE)  can  be  modified  with 
minimal  effort  using  a  DOS  editor  and  used  with  WESHAKE  as  described  in  Part 
IV  and  Appendix  E.  These  changes  include  adding  a  line  to  the  option  defining 
the  earthquake  motion  and  deleting  the  end  of  input  option  and  any  existing 
end-of-file  markers.  Differences  in  format  for  input  files  to  SHAKEl  and 
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WESHAKE  are  marked  for  easy  identification  in  Appendix  E  to  facilitate 
changing  of  existing  input  files . 

Output  files 

43.  Several  different  files  are  created  as  a  result  of  running  WESHAKE. 
Some  of  these  files  are  created  as  a  consequence  of  mandatory  actions  (i.e., 
always  created)  whereas  other  files  are  specific  to  different  options 
selected.  A  summary  table  of  output  files  is  shown  in  Table  III-5.  Output 
files  are  in  ASCII  format  so  they  may  be  read  by  a  DOS  editor,  word  processing 
software,  or  sent  directly  to  a  printer. 

44.  Five  of  the  output  files  are  automatically  created  whenever  WESHAKE 
is  run:  GMOD,  DAMP,  EQIN,  STRESS,  and  OUTPUT.  The  first  three  of  these 
present  some  of  the  input  parameters  into  a  convenient  form  for  plotting  using 
programs  such  as  GRAPHER.  The  format  used  for  these  files  is  comma- separated 
variables.  The  primary  output  file  and  typically  the  largest  in  size  is  named 
OUTPUT.  This  file  summarizes  the  inputs  contained  in  the  specifications  file 
and  all  of  the  results  of  the  various  actions  and  options.  The  file  STRESS 
also  has  a  comma -separated  form  of  peak  effective  shear  stresses  and  shear 
strains  for  the  top  of  each  layer. 

Error  Checks 

45.  Some  error  checks  have  been  Implemented  into  version  1.0  of  WESHAKE 
in  an  attempt  to  assist  the  user  and  preclude  the  calculation  of  erroneous 
results.  On  the  most  basic  level,  a  system  of  recognizing  valid  responses  to 
the  menu  queries  is  used  so  that  many  unintentional  entries  made  will  not  be 
accepted  and  the  user  is  allowed  to  re-enter  the  value.  This  system  is 
intended  to  prevent  the  pre-processing  portion  of  the  program  from  exiting 
abruptly  and  consequently  losing  all  entries  made  up  to  that  point. 

46.  The  implementation  of  WESHAKE  also  has  established  bounds  for  some 
entries.  This  includes  parameters  like  the  coefficient  of  lateral  earth 
pressure,  unit  weight  of  soil  and  pore  fluid  (including  units  of  kef),  and 
shear  wave  velocity.  These  checks  are  Intended  to  help  the  user  find  typing 
errors  and  help  operators  unfamiliar  with  typical  values  for  some  geotechnical 
parameters . 

47.  The  final,  and  most  important,  set  of  error  checks  involve  the 
detection  of  potential  calculation  errors.  Two  primary  checks  are  now  in 
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place:  checks  on  comparisons  between  the  lengths  of  earthquake  records  and 

size  of  FFT  arrays  and  the  level  of  effective  shear  strain  used  to  calculate 
new  moduli  and  damping  ratios  (must  be  less  than  1.0  percent  --  the  extent  of 
definition  of  the  curves).  In  both  cases,  the  user  is  warned  immediately  on 
the  computer  screen. 


Validation  of  WESHAKE 

48.  Although  only  minor  changes  have  been  made  to  the  core  program 
(original  version  of  SHAKE)  the  validation  of  WESHAKE  was  considered  to  be 
necessary.  One  good  source  of  validation  is  a  comparison  of  results  with  the 
example  problem  used  by  Schnabel,  Lysmer,  and  Seed  (1972)  in  chapter  6  of  the 
original  SHAKE  manual.  Minor  interpretations  were  required  to  accommodate  the 
parameters  used  by  WESHAKE  (e.g.,  shear  wave  velocity  in  lieu  of  S„  for 
clays) .  A  summary  of  this  comparison  is  presented  in  Appendix  F  including  a 
soil  column  representing  the  example  problem.  Note  that  this  example 
contained  an  error  in  the  definition  of  the  normalized  shear  modulus 
relationship  for  sands.  A  discontinuity  exists  at  a  shear  strain  of  0.3316 
percent  (shear  strain  should  have  been  defined  as  0.0316).  This  error  was 
repeated  for  validation  to  maintain  consistency. 

49.  In  general,  the  results  between  the  two  are  consistent.  A  few 
comparisons  are  made  within  the  output  file  presented  in  Appendix  F.  Small 
differences  may  be  the  consequence  of  differences  in  machine  accuracy.  Other 
favorable  comparisons  have  been  made  with  other  problems  previously  analyzed 
at  WES. 
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Table  III-l 


Correlation  Between  Potion  Numbers  Used  in  SHAKE  and  WESHAKE 


INPUT/ACTION 

WESHAKE 

SHAKE 

General  information 

"Project  Information" 

"Initialization" 

Soil  column/soil  properties 

Mandatory  Action  1 

Options 

2  and  8 

Earthquake  motion 

Mandatory  Action  2 

Option 

1 

Point  of  excitation 

Mandatory  Action  3 

Option 

3 

Compute  motions 

Mandatory  Action  3 

Option 

4 

Compute  motions  in  sublayers 

User  Option  1 

Option 

5 

Print  object  motion 

User  Option  2 

Option 

6 

Change  object  motion 

tl 

Option 

7 

Compute  response  spectra 

User  Option  3 

Option 

9 

Increase  time  interval 

User  Option  4 

Option 

10 

Decrease  time  interval 

tt 

Option 

11 

Plot  Fourier  spectra  of 
object  motion 

User  Option  5 

Option 

12 

Plot  Fourier  spectra  of 
computed  motion 

It 

Option 

13 

Plot  time  history  of  object 
motion 

User  Option  6 

Option 

14 

Compute  amplification 
spectra 

User  Option  7 

Option 

15 

Compute  stress  or  strain 
history 

User  Option  8 

Option 

16 

Close  input  file 

User  Option  10 

Analysis  Option  1 

Option 

0 
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Table  I 

References  for  Material 


Material  Type 


Gravel 

Sand  (upper  bound,  average 
and  lower  bound) 

Clay  &  Silt: 

PI  -  5-10 
PI  -  10-20 
PI  -  20-40 
PI  -  40-80 
PI  >  80 

Mexico  City  Clay 


*  Confirmed  by  Seed  et  al  (1986) 
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References 


Schnabe 1  (19 73) 
Seed  et  al  (1986) 
Seed  &  Idriss  (1970)* 


Sun,  Golesorkhi, 
&  Seed  (1988) 


Seed  and  Idriss  (1970) 
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Table  III-5 

Summary  of  Output  Files 


File  Name 

Source 

Description 

GMOD 

Mandatory  actions 

Table  of  points  defining  the 
variations  of  shear  modulus  with 
shear  strain  representing  soils  in 
soil  column 

DAMP 

Mandatory  actions 

Table  of  points  defining  the 
variations  of  damping  ratio  with 
shear  strain  representing  soils  in 
soil  column 

EQIN 

Mandatory  actions 

Two-column  table  of  time  and 
acceleration  defining  input 
earthquake  record 

OUTPUT 

Mandatory  actions 

Primary  output  file;  contains  all 
pertinent  information;  printing  of 
earthquake  record  is  option 
[OUTKEY] 

STRESS 

Mandatory  actions 

Two-column  table  of  (top-of-layer) 
depth  and  corresponding  shear 
stress  and  shear  strain 

AMAX 

User  option  1 

Two -column  table  of  (top-of-layer) 
depth  and  corresponding  maximum 
accelerations 

VELSPEC 

User  option  3 

Sub-option  0  or  2 

Table  of  pseudo  velocity  spectra 
showing  periods  and  velocities  at 
specified  levels  of  damping 
(columns) 

ACCSPEC 

User  option  3 

Sub-option  1  or  2 

Table  of  pseudo  acceleration 
spectra  showing  periods  and 
velocities  at  specified  levels  of 
damping  (columns) 
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Standard  relationships  between  normalized  shear  modulus  and 
shear  strain  for  granular  soils  and  rock 


Figure  III-2.  Standard  relationships  between  normalized  shear  modulus  and  shear 
strain  for  cohesive  soils  (Sun,  Golesorhki,  and  Seed  1988) 
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Figure  III- 3.  Standard  relationships  between  damping  ratio  and  shear 
strain  for  granular  soils  and  rock 
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Figure  III -4.  Standard  relationships  between  damping  ratio  and  shear 
strain  for  cohesive  soils  (Sun,  Golesorkhi,  and  Seed  1988) 


PART  IV:  RUNNING  WESHAKE 


50.  The  execution  of  WESHAKE  basically  two  stages:  preprocessing  and 
analysis.  In  the  preprocessing  stage,  the  computer  prompts  the  user  for 
various  mandatory  information  --  the  type  of  earthquake  motion,  the  nuni..>er  of 
soil  layers,  the  type  of  soil  layers,  and  certain  parameters  needed  for  the 
analysis  --  and  user  options  for  additional  analysis.  Output  files  containing 
results  from  the  analysis  are  generated  by  the  mandatory  actions  and  the 
various  options.  These  output  files  are  in  ASCII  format  and,  therefore,  can 
be  viewed  on  a  DOS  editor  or  converted  into  a  word  processing  document.  The 
final  action  posed  in  the  analysis  stage  is  to  allow  the  user  to  generate 
printed  files  of  the  various  output  documents.  Figure  IV- 1  illustrates  the 
processing  cycle  for  WESHAKE.  A  flowchart  of  subroutines  defining  these  two 
stages  is  shown  in  Appendix  G. 

51.  The  format  for  the  various  user  responses  in  Part  IV  will  be  menus 
and  user  response  screens.  Example  responses  are  shown  in  bold  type;  variable 
names  are  enclosed  in  square  brackets.  To  start  WESHAKE,  type  VESHAKE  at  the 
DOS  prompt.  An  example  problem  used  throughout  this  Part  is  shown  in  Figure 
III-2.  The  first  screen  is; 


U.S.  ARMY  WATERWAYS  EXPERIMENT  STATION 
GEOTECHNICAL  LABORATORY 
EARTHQUAKE  ENGINEERING  AND  SEISMOLOGY  BRANCH 
Vicksburg,  Mississippi 

****************** 

WESHAKE  1.0 
****************** 


OPTIMIZED  FOR  USE  ON 
DOS -COMPATIBLE  PERSONAL  COMPUTERS 

by 

David  W.  Sykora  and  Dr.  David  C.  Wallace 
SUMMER  1992 

For  consultation  and  recent  updates,  call: 

(601)  634-3551 

PRESS  <ENTER>  TO  CONTINUE 
<ENTER> 

(CONTINUE  NEXT  SCREEN) 
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This  designation  is  necessary  to  Indicate  whether  a  companion  file  should  be 
sought  (If  "USING  WESHAKE”  is  selected,  the  companion  file  must  exist). 
Mandatory  Actions 

53.  The  pre-processing  stage  now  follows  with  mandatory  actions  which 
are  listed  in  the  SUMMARY  OF  MANDATORY  ACTIONS  screen; 


SUMMARY  OF  MANDATORY  ACTIONS 
**************************** 


********************************************************** 


**  FOLLOWING  THE  INPUT  OF  SOME  INITIAL  INFORMATION,  ** 
**  THE  FOLLOWING  MANDATORY  ACTIONS  WILL  BE  UNDERTAKEN  ** 


**  IN  THE  SPECIFIED  ORDER:  ** 
**  - - - - - -  ** 

**  1  CREATE  SOIL  COLUMN  ** 
**  2  SELECT  EARTHQUAKE  RECORD  ** 
**  3  SPECIFY  POINT  OF  EARTHQUAKE  EXCITATION  ** 


********************************•*"*'*********************** 


PRESS  <ENTER>  TO  CONTINUE 
<ENTER> 

(CONTINUE  NEXT  SCREEN) 


These  actions  are  automatically  enacted  by  the  program  through  the  use  of 
additional  menus  and  screens  as  shown  below.  A  project  title  may  be  entered 
on  the  following  screen; 


PROJECT  INFORMATION 
******************* 


ENTER  A  PROJECT  TITLE  [PTITLE] : 

-  MAXIMUM  OF  74  CHARACTERS  IN  LENGTH 
TITLE  OF  EXAMPLE  PROBLEM 

(CONTINUE  NEXT  SCREEN) 


Next,  mandatory  action  1  is  summarized; 


MANDATORY  ACTION  1:  CREATE  SOIL  COLUMN 
*************************************** 


INFORMATION  AND  PROPERTIES  DEFINING  EACH  SOIL 
LAYER  WILL  NOW  BE  REQUESTED  BEGINNING  AT  THE 
TOP  LAYER  AND  WORKING  DOWN  TO  ROCK  (BASE). 

A  MENU  OF  SOIL  TYPES  WILL  BE  DISPLAYED  FOR 

EACH  SOIL  LAYER  TO  SELECT  REPRESENTATIVE  RELATIONSHIPS 

FOR  SHEAR  MODULUS  AND  DAMPING  RATIO. 


(continued  next  page) 
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This  designation  is  necessary  to  indicate  whether  a  companion  file  should  be 
sought  (If  "USING  WESHAKE”  is  selected,  the  companion  file  must  exist). 
Mandatory  Actions 

53.  The  pre-processing  stage  now  follows  with  mandatory  actions  which 
are  listed  in  the  SUMMARY  OF  MANDATORY  ACTIONS  screen: 


SUMMARY  OF  MANDATORY  ACTIONS 
**************************** 

********************************************************* 

** 

FOLLOWING  THE  INPUT  OF  SOME  INITIAL  INFORMATION, 

** 

** 

THE  FOLLOWING  MANDATORY  ACTIONS  WILL  BE  UNDERTAKEN 

** 

** 

IN 

THE  SPECIFIED  ORDER: 

** 

** 

.  .  ..  . 

** 

** 

1 

CREATE  SOIL  COLUMN 

** 

** 

2 

SELECT  EARTHQUAKE  RECORD 

** 

** 

3 

SPECIFY  POINT  OF  EARTHQUAKE  EXCITATION 

** 

********************************************************* 

PRESS 

<ENTER>  TO  CONTINUE 

<ENTER> 

(CONTINUE  NEXT  SCREEN) 

These  actions  are  automatically  enacted  by  the  program  through  the  use  of 
additional  menus  and  screens  as  shown  below.  A  project  title  may  be  entered 
on  the  following  screen: 


PROJECT  INFORMATION 
★★★★★*★★*★★★***★★** 

ENTER  A  PROJECT  TITLE  [PTITLE]: 

-  MAXIMUM  OF  74  CHARACTERS  IN  LENGTH 
TITLE  OF  EXAMPLE  PROBLEM 

(CONTINUE  NEXT  SCREEN) 


Next,  mandatory  action  1  is  summarized: 


MANDATORY  ACTION  1:  CREATE  SOIL  COLUMN 
*************************************** 


INFORMATION  AND  PROPERTIES  DEFINING  EACH  SOIL 
LAYER  WILL  NOW  BE  REQUESTED  BEGINNING  AT  THE 
TOP  LAYER  AND  WORKING  DOWN  TO  ROCK  (BASE). 

A  MENU  OF  SOIL  TYPES  WILL  BE  DISPLAYED  FOR 

EACH  SOIL  LAYER  TO  SELECT  REPRESENTATIVE  RELATIONSHIPS 

FOR  SHEAR  MODULUS  AND  DAMPING  RATIO. 


(continued  next  page) 
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ONCE  THE  COLUMN  HAS  BEEN  CREATED,  A  SUMMARY 
TABLE  WILL  BE  DISPLAYED  FOR  REVIEW.  ERRORS 
MAY  BE  CORRECTED  AT  THAT  TIME  BY  RE-ENTERING 
ALL  VALUES  FOR  THAT  LAYER. 

St***********************************************-****** 

PRESS  <ENTER>  TO  CONTINUE 
<ENTER> 

(CONTINUE  NEXT  SCREEN) 


54.  The  next  screen  will  ask  general  information  about  the  soil  column 
to  prepare  for  input  of  properties  for  each  layer.  Note  that  the  unit  of 
measurement  for  the  total  unit  weight  of  soil  is  kips  per  cubic  foot.  The 
entry  corresponding  to  the  layer  number  for  the  water  table  is  required  input 
to  proceed  but  is  only  used  for  the  calculation  of  when  values  of  K2  are 
used.  A  unit  weight  of  pore  fluid  corresponding  to  water  (0.0624  kef)  is 
assumed  for  this  calculation. 


ENTER  THE  NUMBER  OF  SOIL  LAYERS  FOR  THE  COLUMN 
INCLUDING  ROCK  [ML] : 


ENTER  THE  LAYER  NUMBER  AT  THE  TOP  OF  WHICH  LIES 
THE  WATER  TABLE  [MWL] : 

5 


ENTER  THE  IDENTIFICATION  FOR  THE  SOIL  PROFILE: 

-  MAXIMUM  OF  43  CHARACTERS  IN  LENGTH 
IDENTIFICATION  OF  EXAMPLE  SOIL  PROFILE 

(CONTINUE  NEXT  SCREEN) 


55.  The  next  three  screens  query  the  user  for  specific  information 
about  each  soil  layer.  The  first  two  screens  have  menus  that  summarize  the 
shear  modulus  and  damping  ratio  data  bases,  respectively. 


DATA  BASE  OF  NORMALIZED  SHEAR  MODULUS  CURVES 


******************************************************* 

**  NUMBER  DESCRIPTION  ** 

4,4.  ^ 

'k'k 

1 

ROCK,  AVERAGE 

(SCHNABEL  1973) 

** 

** 

2 

GRAVEL,  AVERAGE 

(SEED  ET  AL  1986) 

** 

** 

3 

SAND,  LOWER  BOUND 

(SEED  & 

IDRISS  1970) 

** 

** 

4 

SAND,  AVERAGE 

(SEED  & 

IDRISS  1970) 

** 

** 

5 

SAND,  UPPER  BOUND 

(SEED  & 

IDRISS  1970) 

** 

** 

6 

CLAY/SILT,  PI-5 -10 

(SUN  et 

al  1988) 

** 

** 

7 

CLAY/SILT,  PI-10-20 

(SUN  et 

al  1988) 

** 

** 

8 

CLAY/SILT,  PI-20-40 

(SUN  et 

al  1988) 

** 

** 

9 

CLAY/SILT,  PI-40-80 

(SUN  et 

al  1988) 

** 

** 

10 

CLAY/SILT,  PI>80 

(SUN  et 

al  1988) 

** 

** 

11 

MEXICO  CITY  CLAY 

(SUN  et 

al  1988) 

** 

(continued  next  page) 
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******************************************************* 
ENTER  THE  NUMBER  FOR  LAYER  1  OF  5  LAYERS  [TYPE] : 

7 

(CONTINUE  NEXT  SCREEN) 


and 


DATA  BASE  OF  DAMPING  RATIO  CURVES 
******************************************************* 


** 

icic 

NUMBER 

DESCRIPTION 

ick 

** 

1 

ROCK.  AVERAGE 

(SCHNABEL  1973) 

** 

** 

2 

GRAVEL,  AVERAGE 

(SEED  ET  AL  1986) 

** 

** 

3 

SAND,  LOWER  BOUND 

(SEED  & 

IDRISS 

1970) 

** 

** 

4 

SAND,  AVERAGE 

(SEED  & 

IDRISS 

1970) 

** 

** 

5 

SAND,  UPPER  BOUND 

(SEED  & 

IDRISS 

1970) 

•k-k 

** 

6 

CLAY/SILT, LOWER  BOUND 

(SEED  & 

IDRISS 

1970) 

** 

7 

CLAY/SILT, AVERAGE 

(SEED  & 

IDRISS 

1970) 

** 

** 

8 

CLAY/SILT, UPPER  BOUND 

(SEED  & 

IDRISS 

1970) 

** 

******************************************************* 


ENTER  THE  NUMBER  FOR  LAYER  1  OF  5  LAYERS  (TYPE): 

7 

(CONTINUE  NEXT  SCREEN) 


Site -specific  shear  modulus  and  damping  relations  can  be  added  (or  deleted) 
from  the  soil  property  data  bases  (SHEARDB  and  DAMPDB)  using  a  DOS  editor. 
There  are  two  requirements  for  customizing  the  data  bases:  use  the  exact 
format  as  documented  in  Appendix  B  and  provide  a  unique  identification  number 
for  the  curve,  INUM.  Note  that  any  new  entries  to  the  data  base  will  not  be 
displayed  on  the  computer  screen  even  though  they  can  be  accessed. 

56.  The  third  screen  is  used  to  query  the  user  about  properties  for 
each  layer.  The  use  of  K2  or  V,  to  calculate  shear  modulus  was  described  in 
Part  II  of  this  report.  The  example  shown  is  for  V,  input.  If  Kj  input  is 
specified,  the  appropriate  prompts  for  initial  and  maximum  K2  ,  (K2)max  >  will 
be  displayed  and  a  value  of  coefficient  of  lateral  earth  pressure  will  be 
requested.  Only  shear  wave  velocities  are  allowed  to  define  shear  modulus  for 
rock. 

ENTER  THE  THICKNESS  (ft)  OF  LAYER  1  [HL] : 

5 

THE  INITIAL  ESTIMATES  FOR  DAMPING  RATIO  [B]  WILL  BE 

5  PERCENT  FOR  SOIL  AND  2  PERCENT  FOR  ROCK. 

ENTER  THE  TOTAL  UNIT  WEIGHT  (kef)  OF  LAYER  1  [W] : 

0.110 

(continued  next  page) 
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DESIGNATE  METHOD  OF  MODULUS  CALCULATION 
FOR  SHEAR  VELOCITY,  ENTER  VS 
FOR  K2,  ENTER  K2 

VS 

ENTER  THE  SHEAR  WAVE  VELOCITY  (fps) 

FOR  LAYER  1 : 

425 

(CONTINUE  NEXT  SCREEN) 


57 .  Once  the  soil  column  has  been  created  a  summary  table  is  displayed 
to  allow  the  user  to  check  the  information.  The  summary  table  for  the  example 
problem  is: 


LAYER  SOIL  SOIL  S;  SHEAR 

LAYER 

EARTH 

UNIT 

NUM 

NUM 

ID  D ;  DAMPING 

THICK 

PRESS 

WEIGHT 

STIFFNESS 

(ft) 

(kef) 

1 

9 

S:  CLAY/SILT, PI-40 -80 

5.0 

.110 

Vs  =  425. 

7 

D;  CLAY/SILT,  Average 

2 

6 

S;  CLAY/SILT, PI-  5-10 

8.0 

.115 

Vs  -  575. 

7 

D:  CLAY/SILT,  Average 

3 

4 

S :  SAND ,  Average 

10.0 

.38 

.125 

K2  -  35. 

4 

D:  SAND,  Average 

4 

3 

S;  SAND,  Lower  Bound 

15.0 

.130 

Vs  -  1100. 

4 

D:  SAND,  Average 

5 

1 

S :  ROCK ,  Average 

.150 

Vs  -  5000. 

1 

D;  ROCK,  Average 

IS 

THIS 

SOIL  COLUMN  CORRECT? 

ENTER 

0  FOR  NO 

1 

ENTER 

1  FOR  YES 

TO 

CONTINUE  ON  TO  MANDATORY  ACTION 

2,  PRESS  <ENTER> 

<ENTER> 

(CONTINUE  NEXT  SCREEN) 

If  a  correction  is  required,  the  user  may  enter  the  layer  nirmber  at  the  prompt 
and  then  all  information  for  the  designated  layer  number  must  be  input  again 
via  computer  prompt.  Once  the  correction(s)  has  been  made  the  computer  will 
again  display  the  soil  profile  and  ask  the  user  if  the  soil  profile  is 
correct.  The  process  will  repeat  until  the  column  is  accepted  by  the  user 
(indicated  by  entering  1). 

58.  This  marks  the  end  of  defining  soil  stratigraphy  and  soil 
properties.  If  the  user  specified  the  option  for  creation  of  a  specification 
file,  the  program  continues  directly  to  the  next  mandatory  action  which  is  the 
selection  of  an  earthquake  record.  If  the  user  specified  the  option  for 
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creation  of  a  soil  profile,  the  information  collected  up  to  this  point  is 
saved  first  before  proceeding  to  the  next  mandatory  action. 

59 .  The  following  menu  is  used  to  initiate  mandatory  action  2 . 
Earthquake  records  in  the  data  base  (EARTHQ)  are  displayed  for  selection: 


MANDATORY  ACTION  2 :  SELECT  EARTHQUAKE  RECORD 
******************************************** 
EARTHQUAKE  DATA  BASE 


** 

NO. 

MEASURED  RECORD 

NO. 

SYNTHETIC 

RECORD 

** 

** 

** 

1 

GOLDEN  GATE  1957 

14 

FOLSOM  RECORD  "A" 

** 

** 

2 

PARKFIELD  1966 

15 

FOLSOM  RECORD  "B" 

** 

** 

3 

CASTAIC  RIDGE  1971 

16 

NEW  MADRID 

500-YR  HI 

** 

4 

LAKE  HUGHES  #  4  1971 

17 

NEW  MADRID 

500 -YR  H2 

** 

** 

5 

SITKA  1972 

18 

NEW  MADRID 

1000-YR 

HI 

** 

6 

GILROY  #1  1974 

19 

NEW  MADRID 

1000 -YR 

H2 

** 

** 

7 

GILROY  #1  1979 

20 

NEW  MADRID 

5000-YR 

HI 

** 

** 

8 

SUPERSTITION  1979 

21 

NEW  MADRID 

5000-YR 

H2 

** 

** 

9 

SUPERSTITION  1981 

22 

RIRIE  DAM 

** 

** 

10 

GILROY  #1  1984 

** 

** 

11 

IVERSON  1985 

** 

** 

12 

SLIDE  MT  1985 

"kit 

** 

13 

HOLLISTER  AIRPORT  1989 

** 

** 

** 

** 

Other: 

** 

** 

25 

STOP  AT  THIS  POINT  AND  SAVE  SOIL  PROFILE  FILE 

** 

** 

26 

PICK  YOUR  OWN  EARTHQUAKE  MOTION 

** 

ENTER  EARTHQUAKE  NUMBER: 

1  (CONTINUE  ON  NEXT  SCREEN) 


Specific  parameters  describing  these  earthquakes  are  listed  in  Tables  III-3 
and  III -4.  Listings  of  the  data  files  and  plots  of  the  records  and  velocity 
spectra  are  provided  in  Appendix  C.  If  a  record  is  selected  (i.e.,  option  25 
is  not  selected)  then  specify  if  the  values  of  acceleration  should  be  echoed 
(printed)  in  the  OUTPUT  file: 


SPECIFY  WHETHER  RECORD  DATA  ARE  TO  BE  PRINTED 
IN  OUTPUT  FILE  [OUTKEY] : 

ENTER  0  FOR  NO 
ENTER  1  FOR  YES 

0 

(CONTINUE  ON  NEXT  SCREEN) 


60.  Earthquake  records  not  in  the  data  base  may  be  selected  by  entering 
30  and  then  entering  the  name  of  the  file  (which  must  reside  in  the  current 
working  directory) .  The  format  for  the  earthquake  record  must  follow  the 
syntax  for  the  earthquake  data  base  which  is  specified  in  Appendix  C.  This 
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format  has  as  a  second  line  "earthquake  characteristics."  This  may  be  left 
blank,  but  the  line  must  exist.  Also,  make  sure  that  no  end-of-file  markers 
exist  (typically  '‘Z) . 

61.  After  the  earthquake  record  has  been  selected,  the  earthquake 
characteristics  and  the  maximum  acceleration  are  displayed  for  convenience  and 
options  for  scaling  the  record  are  posed.  The  accelerogram  can  be  scaled 
linearly  by  specifying  a  scaling  factor  (all  values  of  acceleration  are 
multiplied  by  this  constant)  or  a  maximum  value  of  acceleration,  a,.,  (all 
values  of  acceleration  are  multiplied  by  the  ratio:  (a....)„.„  /  (a,...)„i,<  ). 

The  screen  for  this  option  is: 


EARTHQUAKE  MOTION  CHARACTERISTICS: 
*************************************************** 
"Mag-7.1,  Dis-69  km,  Amax-.63  g. ,  Rock  Outcrop" 

MAXIMUM  CUT-OFF  FREQUENCY:  50.0 
******************************************************** 

DO  YOU  WANT  TO  SCALE  THE  RECORD? 

ENTER  0  NO 
ENTER  1  YES 

1 

SCALING  RECORD  OPTION: 

ENTER  0  TO  USE  SCALING  FACTOR 

ENTER  1  TO  SET  NEW  MAXIMUM  ACCELERATION 

0 

ENTER  THE  SCALING  FACTOR  [XF] : 

2.0 

TO  CONTINUE  ON  TO  MANDATORY  ACTION  3,  PRESS  <ENTER> 
<ENTER> 

(CONTINUE  NEXT  SCREEN) 


The  default  scaling  factor  on  the  earthquake  motion  record  is  set  to  1.00  and 
the  default  maximum  acceleration  is  set  equal  to  the  peak  measured  value  (also 
accomplished  by  setting  XMAX  -  0.0).  The  complete  accelerogram  can  be  altered 
by  selecting  option  1  (YES).  If  option  2  (NO)  is  selected,  then  the 
parameters  are  set  to  the  earthquake  motion  record. 

62.  The  final  mandatory  action  is  used  to  define  the  point  in  the  soil 
profile  at  which  the  input  motion  corresponds.  The  previous  screen  and  Tables 
III-3  and  III-4  list  where  the  motion  corresponds  (typically  rock  outcrop). 

In  most  cases,  particularly  for  design  or  seismic  stability  analysis,  rock 
outcrop  motions  are  used.  To  use  the  motion  as  an  outcrop  motion,  select  the 
layer  number  corresponding  to  the  base  rock  (last  layer)  and  specify  an 
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outcropping  layer  (type  0).  For  base  rock  motions  (refer  to  Figure  II-l), 
again  use  the  base  rock  layer  number  and  then  type  a  1  for  the  motion  type. 

If  the  motion  corresponds  to  some  other  layer,  for  example  ground  surface 
(free  field),  use  the  layer  number  (which  corresponds  to  the  top  of  the  layer) 
and  a  1.  The  screen  for  these  actions  is: 


MANDATORY  ACTION  3:  SPECIFY  OBJECT  MOTION 
***************************************** 

******************************************************* 
SPECIFY  A  LAYER  FOR  THE  OBJECT  MOTION  [IN]: 

(ROCK  IS  LAYER  5) 


ENTER  THE  MOTION  TYPE  [INT] : 

0  -  FOR  OUTCROPPING  LAYER 
1  -  FOR  LAYER  WITHIN  SOIL  PROFILE 


(CONTINUE  ON  NEXT  SCREEN) 


63.  Parameters  defining  the  convergence  criteria  have  been  preset  in 
WESHAKE  for  the  convenience  of  the  user.  These  values  are  displayed,  as  shown 
below,  to  remind  the  user: 


THE  FOLLOWING  VALUES  HAVE  BEEN  ASSUMED: 

THE  MAXIMUM  NUMBER  OF  ITERATIONS  [ITMAX]  -  20 

THE  ACCEPTABLE  DIFFERENCE  BETWEEN  THE 
LAST-USED  MODULUS  AND  DAMPING  VALUES,  AND  THE 
STRAIN  COMPATIBLE  VALUES  [ERR]  -  5  PERCENT 

THE  RATIO  BETWEEN  EFFECTIVE  STRAIN  AND 
MAXIMUM  STRAIN  [PRMUL]  -  65  PERCENT 

TO  CONTINUE  ON  TO  USER  OPTION  MENU,  PRESS  <ENTER> 
<ENTER> 

(CONTINUE  NEXT  SCREEN) 


The  user  may  change  these  values  once  the  specifications  file  has  been 
completed  (see  USER  OPTIONS  MAIN  MENU  in  next  section)  by  exiting  WESHAKE, 
editing  the  specific  file  using  a  DOS  editor,  and  restarting  WESHAKE  (and 
using  an  existing  specification  file). 

64.  Once  these  steps  have  been  completed,  the  first  component  of  the 
specifications  file  required  to  run  WESHAKE  has  been  created.  This  file  can 
be  used  as  often  as  necessary.  The  format  of  the  specifications  file  for  all 
of  these  responses  is  provided  in  Appendix  E.  An  example  specifications  file 


48 


is  shown  in  Appendix  H  corresponding  to  the  example  problem  shown  in  Figure 
IV- 3.  WESHAKE  continues  by  proceeding  to  the  USER  OPTIONS  MAIN  MENU  which  is 
described  below. 

User  Options 

65 .  The  USER  OPTIONS  MAIN  MENU  is  encountered  following  mandatory 
actions  for  new  specifications  files  or  immediately  following  selection  of  an 
existing  specifications  file.  This  menu  is  still  part  of  the  pre-processing 
stage  and  is  used  to  identify  various  options  and  outputs  available  to  the 
user.  It  is  also  used  to  save  the  specifications  file  and  exit  WESHAKE  or 
proceed  to  the  analysis  stage.  Each  option  on  the  menu  will  have  its  own  set 
of  screens  that  will  prompt  the  user  for  information  needed  to  implement  that 
option.  The  menu  is; 


USER  OPTION  MAIN  MENU 
********************* 

******************************************************* 

**  OPTION  DESCRIPTION  ** 

_ _  _  _  _  _ 

** 

1 

COMPUTE  MOTION  IN  SPECIFIED  SUBLAYERS 

** 

** 

2 

PRINT  OR  CHANGE  OBJECT  MOTION 

** 

** 

3 

COMPUTE  RESPONSE  SPECTRA 

** 

** 

4 

INCREASE  OR  DECREASE  TIME  INTERVAL 

** 

** 

5 

PLOT  FOURIER  SPECTRUM  OF  OBJECT  OR 

** 

** 

COMPUTED  MOTION 

** 

** 

6 

PLOT  TIME  HISTORY  OF  OBJECT  MOTION 

** 

** 

7 

COMPUTE  AMPLIFICATION  SPECTRUM 

** 

** 

8 

COMPUTE  STRESS/STRAIN  HISTORY 

** 

I  ★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★iHr********************** 

** 

9 

SAVE  FILE  AND  RETURN  TO  DOS 

** 

** 

10 

SAVE  FILE  AND  PROCEED  TO  ANALYSIS  STAGE  ** 

ENTER  ONE 

OPTION  NUMBER: 

1 

(CONTINUE  NEXT 

SCREEN) 

66.  The  first  eight  options  listed  in  the  USER  OPTIONS  MAIN  MENU  are 
described  briefly  with  user  menus  in  the  following  subsections.  These  options 
may  be  used  in  any  order  (although  order  of  options  2  and  4  will  affect 
subsequent  results)  and  repeated  as  often  as  desired.  Theoretical  details 
about  these  options  may  be  found  in  the  manual  for  SHAKE  (Schnabel,  Lysmer, 
and  Seed  1972)  and  in  textbooks  on  soil  and  structural  dynamics.  The 
correspondence  between  WESHAKE  and  SHAKE  option  numbers  is  presented  in 
Table  III-l. 


67.  The  last  two  selections  in  the  USER  OPTIONS  MAIN  MENU  pertain  to 
actions  once  all  of  the  program  options  have  been  chosen.  Option  9,  is  used 


to  save  the  current  specifications  file  and  exit  WESHAKE  (without  placing  the 
file  termination  statement) .  Option  10  is  used  to  save  the  file  and  continue 
on  for  execution  (adding  termination  statement:  OPTION  0:  END  OF  INPUT). 
Specification  files  that  are  reused  (Option  4  in  SPECIFICATION  FILE 
DESIGNATION  menu)  are  automatically  stripped  of  the  termination  statement. 

68 .  Compute  Motions  in  Sublayers .  Option  1  provides  for  information 
regarding  the  horizontal  accelerations  at  layers  of  interest.  Either  tables 
of  peak  values  or  the  complete  time  record  of  acceleration  may  be  chosen.  The 
screen  for  the  OPTION  1  MENU  is  shown  on  the  next  page . 


OPTION  1  MENU 
************* 

THIS  OPTION  IS  USED  TO  SPECIFY  THE  EXTENT 
OF  OUTPUT  FOR  ACCELERATIONS  AT  THE  TOP  OF 
SPECIFIED  LAYERS 

******************************. V************************ 

ENTER  0  TO  GET  MAXIMUM  ACCELERATION  ONLY 
ENTER  1  TO  GET  VARIATION  OF  ACCELERATION  WITH 
TIME  AND  PEAK  ACCELERATIONS 

0 

THE  PRESENT  SOIL  COLUMN  HAS  5  LAYERS 

THE  OBJECT  MOTION  HAS  BEEN  ASSIGNED  TO  LAYER  5 

SPECIFY  AT  WHICH  LAYERS  OF  THE  SOIL  COLUMN  THE 
ACCELERATIONS  SHOULD  BE  CALCULATED: 

-  ONE  AT  A  TIME 

-  MAXIMUM  OF  15 

ENTER  THE  LAYER  NUMBER  [LL5(1)1: 

1 

SPECIFY  THE  MOTION  TYPE  [LT5(1)]: 

ENTER  0  FOR  OUTCROPPING 
ENTER  1  FOR  WITHIN  SOIL  PROFILE 

1 

TO  COMPUTE  ACCELERATIONS  AT  MORE  LAYERS 
ENTER  0  FOR  NO 
ENTER  1  FOR  YES 

0 

(RETURN  TO  USER  OPTIONS  MAIN  MENU) 


69.  Print/Change  Motion  in  Sublayer  Options.  Option  2  allows  the  user 
to  print  the  object  motion  in  the  OUTPUT  file  or  change  the  object  motion  for 
recalculation  of  WESHAKE.  (Note  that  the  option  for  printing  the  input  motion 
is  also  contained  in  Mandatory  Action  2  [OUTKEY].)  The  initial  screen  for 
this  option  is  (next  page): 
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OPTION  2  MENU 

ic'k'k’k'kic^'k'kic’k'k'Jc 


THIS  OPTION  IS  USED  TO  PRINT  OR  CHANGE  THE  OBJECT 
(EARTHQUAKE)  MOTION 

ENTER  1  TO  PRINT  OBJECT  MOTION 
ENTER  2  TO  CHANGE  OBJECT  MOTION 

1 

(CONTINUE  NEXT  SCREEN) 


If  the  print  motion  option  is  selected,  the  following  screen  is  shown: 


OPTION  2  SUBMENU  (PRINT  M..TION) 

St  St  St  St  St  St  St  St  St  St  St  St  St  St -it  St  St  St  St  St  St  St  St  st  St -st  St -st  St  St -it 


ENTER  0  TO  PRINT  MAXIMUM  ACCELERATION  IN  OUTPUT  FILE 
ENTER  1  TO  PUT  OBJECT  MOTION  IN  PUNCH  FILE 
ENTER  2  TO  DO  BOTH 

0 

(RETURN  TO  USER  OPTIONS  MAIN  MENU) 


The  PUNCH  file  contains  output  data  without  annotation  (as  compared  with  the 
OUTPUT  file).  Notice  that  in  WESHAKE,  the  option  of  printing  the  object 
motion  in  the  OUTPUT  file  is  also  provided  in  Mandatory  Action  2.  If  the 
change  motion  option  is  selected,  the  following  screen: 


OPTION  2  SUBMENU  (CHANGE  MOTION) 
******************************** 

THIS  IS  USED  TO  CHANGE  PARAMETERS  SET 
IN  MANDATORY  ACTION  2:  SELECT  EARTHQUAKE  RECORD 
******************************************************** 
THERE  ARE  5  LAYERS  IN  THE  SOIL  COLUMN. 

LAYER  5  IS  CURRENTLY  THE  POINT  OF  OBJECT  MOTION 

SPECIFY  THE  LAYER  FOR  OBJECT  MOTION  [LLl]: 

ENTER  0  TO  KEEP  THE  MOTION  AT  THE  SAME  LAYER 
OTHERWISE  ENTER  NEW  LAYER  NUMBER 

I 

IDENTIFY  THE  TYPE  OF  ABOVE  LAYER  [LTl]: 

ENTER  0  FOR  OUTCROPPING 
ENTER  1  FOR  WITHIN  SOIL  PROFILE 


ENTER  THE  TIME  STEP  OF  OBJECT  MOTION  [DTNEW] : 

0.02 

SCALING  OF  ACCELERATIONS  [XF] : 

ENTER  1.  FOR  MAXIMUM  VALUE  OF  RECORD 
OTHERWISE  ENTER  SCALING  FACTOR  (DECIMAL) 

1. 

(RETURN  TO  USER  OPTIONS  MAIN  MENU) 
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70.  Compute  response  spectrum.  Option  3  should  be  selected  to  calculate 
pseudo-velocity  or  acceleration  response  spectrvun.  A  response  spectrum  is  the 
response  of  an  equivalent  damped  single -degree -of -freedom  (SDOF)  system  to  the 
free -field  motion.  The  step-by-step  method  is  used  to  calculate  the  response 
spectrum  in  WFSHAKE.  For  response  acceleration,  absolute  rather  than  relative 
values  are  preferred  (Weigel  1970) .  The  velocity  spectrum  typically  is  used 
for  design  and  analysis  by  structural  engineers. 

71.  The  calculation  of  pseudo  response  spectrum  for  the  input 
earthquake  motion  alone  does  not  require  execution  of  the  program.  This 
option  can  be  inserted  manually  following  the  specification  of  the  input 
motion  and  where  the  motion  occurs  (mandatory  actions  2  &  3). 

72.  The  ratio  of  the  spectrum  of  free-field  ground  surface  acceleration 
spectrum  to  rock  outcrop  acceleration  spectrum  is  typically  desired.  The 
variation  of  this  ratio  with  period  at  five  levels  of  system  damping  will  be 
used  for  design  and  seismic  stability  evaluations. 

73.  The  OPTION  3  MENU  is  displayed  on  two  screens  below: 


OPTION  3  MENU 
************* 

THIS  OPTION  IS  USED  TO  COMPUTE  THE  RESPONSE  SPECTRUM. 
ACCELERATION  OR  VELOCITY  SPECTRUM  CAN  BE  CALCULATED, 

OR  BOTH,  FOR  ANY  OR  ALL  LAYERS. 

******************************************************* 

THERE  ARE  5  LAYERS  IN  THE  SOIL  COLUMN 

LAYER  5  IS  CURRENTLY  THE  POINT  OF  OBJECT  MOTION 

ENTER  THE  LAYER  NUMBER  FOR  ANALYSIS  [LLl] : 

1 

IDENTIFY  THE  TYPE  OF  ABOVE  LAYER  [LTl] ; 

ENTER  0  FOR  OUTCROPPING 

ENTER  1  FOR  WITHIN  THE  SOIL  PROFILE 

1 

ENTER  THE  NUMBER  OF  DAMPING  VALUES  DESIRED  [ND] : 

-  MAXIMUM  OF  5 

2 

SELECT  THE  PARAMETER(S)  OF  INTEREST  [KAV] : 

ENTER  0  FOR  VELOCITY  SPECTRUM 
ENTER  1  FOR  ACCELERATION  SPECTRUM 
ENTER  2  FOR  BOTH 

2 

(CONTINUE  NEXT  SCREEN) 


and  (next  page) : 
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OPTION  3  MENU  (continued) 

******************************************************* 
SELECT  TIME  PERIODS  FOR  COMPUTATIONS 
KPER  -  0 


9 

STEPS 

FROM 

.1 

SEC 

TO 

1. 

SEC 

5 

STEPS 

FROM 

1. 

SEC 

TO 

2. 

SEC 

4 

KPER 

STEPS 

-  1 

FROM 

2. 

SEC 

TO 

4. 

SEC 

18 

STEPS 

FROM 

.1 

SEC 

TO 

1. 

SEC 

10 

STEPS 

FROM 

1. 

SEC 

TO 

2. 

SEC 

8 

KPER 

STEPS 

-  2 

FROM 

2. 

SEC 

TO 

4. 

SEC 

38 

STEPS 

FROM 

.05 

SEC 

TO 

1. 

SEC 

20 

STEPS 

FROM 

1. 

SEC 

TO 

2. 

SEC 

30 

STEPS 

FROM 

2. 

SEC 

TO 

5. 

SEC 

KPER  “  3 

LOGARITHMIC  INCREMENTS  WITH  10 
STEPS  IN  EACH  LOG.  UNIT  FROM 
.1  TO  5. 

KPER  -  4 

LOGARITHMIC  INCREMENTS  WITH  25 
STEPS  IN  EACH  LOG.  UNIT  FROM 
.05  TO  10. 

******************************************************* 
ENTER  VALUE  OF  KPER: 

1 

ENTER  THE  2  VALUES  OF  DAMPING  RATIO  [ZLD] 

ON  SEPARATE  LINES  BELOW: 

-  DECIMAL  FORM 

-  ASCENDING  ORDER 

0.02 

0.05 

(RETURN  TO  USER  OPTIONS  MAIN  MENU) 


74.  Increase/decrease  time  interval.  The  time  interval  of  the 
earthquake  record  can  be  increased  or  decreased  and  the  analysis  rerun.  The 
first  screen  is: 


OPTION  4  MENU 
************* 


THIS  OPTION  IS  USED  TO  INCREASE  OR  DECREASE 
THE  TIME  INTERVAL  FOR  THE  OBJECT  MOTION  BY  USING 
DIFFERENT  NUMBERS  OF  POINTS  IN  THE  FOURIER  TRANSFORM 
****************************************************** 

ENTER  1  TO  INCREASE  THE  TIME  INTERVAL 
ENTER  2  TO  DECREASE  THE  TIME  INTERVAL 

1 

(CONTINUE  WITH  APPROPRIATE  SUBMENU  SCREEN) 


If  an  increase  in  time  is  selected,  the  following  submenu  appears  (next  page) 
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OPTION  4  SUBMENU  (INCREASE  TIME  INTERVAL) 
***************************************** 

ENTER  THE  MULTIPLE  FOR  TIME  INCREASE  [IFR]: 

-  MUST  BE  A  POWER  OF  2 
I 

(RETURN  TO  USER  OPTIONS  MAIN  MENU) 


If  a  decrease  in  time  is  selected,  the  following  submenu  appears: 


OPTION  4  SUBMENU  (DECREASE  TIME  INTERVAL) 
***************************************** 


ENTER  THE  MULTIPLE  FOR  TIME  DECREASE  [IFR]: 

-  MUST  BE  A  POWER  OF  2 

2 

(RETURN  TO  USER  OPTIONS  MAIN  MENU) 


75.  Plot  Fourier  spectrum  of  obiect  or  computed  motion.  The  Fourier 
spectrum  is  calculated  and  plotted  in  the  OUTPUT  file  using  this  option.  The 
first  screen  for  this  option  is: 


OPTION  5  MENU 
************* 

THIS  OPTION  IS  USED  TO  CALCULATE  AND  PLOT  THE 
FOURIER  SPECTRUM  OF  THE  OBJECT  OR  COMPUTED  MOTION 
************************************************** 
ENTER  1  FOR  OBJECT  MOTION 
ENTER  2  FOR  COMPUTED  MOTION 

1 


(CONTINUE  WITH  APPROPRIATE  SUBMENU  SCREEN) 


If  the  object  motion  is  selected,  the  following  submenu  is  shown: 


OPTION  5  SUBMENU  (OBJECT  MOTION) 
******************************** 


ENTER  0  TO  STORE  SPECTRUM  FOR  LATER  PLOTTING: 

-  LIMIT  2  SPECTRA 

ENTER  1  TO  PLOT  ALL  STORED  SPECTRA 

1 

ENTER  THE  NUMBER  OF  TIMES  SPECTRUM  IS  TO  BE 
SMOOTHED  [LNSW]: 

2 

ENTER  THE  NUMBER  OF  VALUES  TO  BE  PLOTTED  [N] : 

-  MAXIMUM  OF  2049 

100 

(RETURN  TO  USER  OPTIONS  MAIN  MENU) 
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The  plot  is  character- style  plot  with  limited  resolution.  The  number  of 
plotted  values  corresponds  to  the  first  N  values  in  the  array.  If  the 
computed  motion  is  selected,  the  following  submenu  is  shown: 


OPTION  5  SUBMENU  (COMPUTED  MOTION) 
********************************** 

ENTER  THE  SUBLAYER  NUMBER  [LL] : 

(ROCK  IS  LAYER  5) 

1 

SPECIFY  THE  TYPE  OF  MOTION  [LT] : 

ENTER  0  FOR  OUTCROPPING 

ENTER  1  FOR  WITHIN  THE  SOIL  PROFILE 


ENTER  0  TO  STORE  SPECTRUM  FOR  LATER  PLOTTING 

-  LIMIT  2  SPECTRA 

ENTER  1  TO  PLOT  ALL  STORED  SPECTRA 

1 

ENTER  THE  NUMBER  OF  TIMES  SPECTRUM  IS  TO  BE 
SMOOTHED  [NSW] ; 

2 

ENTER  THE  VALUES  TO  BE  PLOTTED  [LLL] 

-  MAXIMUM  OF  2049 

100 

(RETURN  TO  USER  OPTIONS  MAIN  MENU) 


76.  Plot  time  history  of  object  motion.  This  option  plots  the 
variation  of  accelerations  with  time  for  specified  sublayers  in  the  OUTPUT 
file.  The  plot  is  character- style  plot  with  limited  resolution.  The  screen 
for  this  option  is: 


OPTION  6  MENU 
************* 


THIS  OPTION  IS  USED  TO  COMPUTE  AND  PLOT  THE 
VARIATION  OF  ACCELERATION  WITH  TIME  OF  THE 
OBJECT  MOTION 

******************************************************* 
SPECIFY  THE  CODE  FOR  VALUES  TO  BE  PLOTTED  [NSKIP] : 

ENTER  0  TO  PLOT  EVERY  VALUE 
ENTER  1  TO  PLOT  EVERY  SECOND  VALUE 
ENTER  2  TO  PLOT  EVERY  THIRD  VALUE 
ENTER  3  TO  PLOT  EVERY  FOURTH  VALUE 
ETC. 


ENTER  THE  NUMBER  OF  VALUES  TO  BE  PLOTTED  [NN] : 

-  MAXIMUM  OF  2049 
2049 

(RETURN  TO  USER  OPTIONS  MAIN  MENU) 
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77.  Amplification  spectrum.  The  amplification  spectrum  between  the 
motions  at  two  layers  can  be  calculated  and  plotted  with  this  option.  The 
plot  is  character-style  plot  with  limited  resolution.  The  screen  is: 


OPTION  7  MENU 
************* 


THIS  OPTION  IS  USED  TO  COMPUTE  THE  AMPLIFICATION 
SPECTRUM  BETWEEN  ANY  TWO  LAYERS. 

THERE  ARE  5  LAYERS  IN  THE  PRESENT  COLUMN. 
******************************************************* 
ENTER  THE  NUMBER  OF  THE  FIRST  LAYER  [LIN]: 

1 

SPECIFY  THE  TYPE  OF  MOTION  FOR  THIS  LAYER  [LINT]: 

ENTER  0  FOR  OUTCROPPING 

ENTER  1  FOR  LAYER  WITHIN  SOIL  PROFILE 


ENTER  THE  NUMBER  OF  THE  SECOND  LAYER  [LOUT] : 

5 

SPECIFY  THE  TYPE  OF  MOTION  FOR  THIS  LAYER  [LOTP] : 

ENTER  0  FOR  OUTCROPPING 

ENTER  1  FOR  LAYER  WITHIN  SOIL  PROFILE 

0 

SELECT  THE  TYPE  OF  PLOTTING  [KP] ; 

ENTER  0  TO  STORE  SPECTRUM  FOR  LATER  PLOTTING 
-  LIMIT  2 

ENTER  1  TO  PLOT  ALL  SPECTRA  STORED  SINCE  LAST 
PLOTTING 

1 

THE  AMPLIFICATION  FACTOR  IS  COMPUTED  FOR  THE 
FIRST  200  FREQUENCIES  WITH  INTERVAL  DFA  (Hz) 

BEGINNING  AT  0. 

ENTER  THE  FREQUENCY  STEPS  [DFA]: 

0.050 

ENTER  AN  IDENTIFICATION  -  LIMIT  40  CHARACTERS: 

SITE  1 

(RETURN  TO  USER  OPTIONS  MAIN  MENU) 


78.  Compute  stress /strain  history.  The  variation  of  peak  shear  stress 
or  peak  shear  strain  with  time  at  the  top  of  a  layer  may  be  calculated  with 
Option  8.  This  option  allows  for  the  specification  of  two  layers  at  once.  If 
the  calculation  is  desired  for  more  than  two  layers,  the  option  must  be  called 
again.  The  results  of  this  option  are  sent  to  the  OUTPUT  file.  The  variation 
of  effective  shear  stress  or  strain  may  be  determined  by  multiplying  these 
peak  values  by  PRMUL.  Note  that  effective  shear  strains  at  the  mid-height  of 
layers  are  used  to  adjust  shear  modulus  and  damping  values  in  the  iterative 
process.  The  OPTION  8  MENU  is  (next  page): 
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OPTION  8  MENU 
************* 


THIS  OPTION  IS  USED  TO  COMPUTE  THE  VARIATION  OF 
SHEAR  STRESS  OR  SHEAR  STRAIN  AT  THE  TOP  OF 
EITHER  ONE  OR  TWO  LAYERS 

******************************************************** 
ENTER  THE  LAYER  NUMBER  CORRESPONDING  TO  THE 
FIRST  LAYER  [LLL] : 

2 

SPECIFY  THE  PARAMETER  OF  INTEREST  [LLGS] : 

ENTER  0  FOR  SHEAR  STRAIN 
ENTER  1  FOR  SHEAR  STRESS 


ENTER  THE  NUMBER  OF  VALUES  TO  BE  PLOTTED  [LNV] : 
-  LIMIT  2049 
2049 

ENTER  AN  IDENTIFICATION  FOR  THIS  PLOT 
EXAMPLE  PLOT 

SPECIFY  THE  SCALE  FOR  PLOTTING  ( i . e MAXIMUM 
VALUE  OF  ORDINATE)  [SK] : 

ENTER  0  FOR  MAXIMUM  VALUE  OF  DATA 
ENTER  1  TO  SPECIFY  MAXIMUM  VALUE 


ENTER  THE  LAYER  NUMBER  CORRESPONDING  TO  THE 
SECOND  LAYER  (OR  LEAVE  BLANK)  [LLL] : 

4 

(RETURN  TO  USER  OPTIONS  MAIN  MENU) 


79.  Summary  of  Options.  After  entering  desired  options  from  the  USER 
OPTIONS  MAIN  MENU  and  option  9,  a  specifications  file  has  been  written  and  the 
pre-processing  stage  is  complete.  The  latter  options  of  the  USER  OPTIONS  MAIN 
MENU  are  then  used  to  proceed  with  steps  to  run  the  program  (analysis  stage, 
described  in  the  next  section) .  Although  a  user  interface  has  been  written 
for  WESHAKE  to  allow  easy  selection  of  user  options,  the  user  may  find  that  a 
DOS  editor  is  more  versatile  and  quicker. 

Analysis  Stage 

80.  The  analysis  stage  consists  of  calculation  of  the  solution  through 
the  execution  of  the  modified  SHAKEl  subprogram  within  the  WESHAKE  package, 
creation  of  output  files,  and  printing  of  output  files  from  within  the  WESHAKE 
shell.  The  solution  process  is  started  by  selecting  option  10  from  the  USER 
OPTIONS  MAIN  MENU.  The  ANALYSIS  STAGE  MAIN  MENU  is  then  shown  (next  page): 
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ANALYSIS  STAGE  MAIN  MENU 


**  OPTION  DESCRIPTION  ** 

**  -  - - -  ** 

**  1  CLOSE  FILE  AND  RUN  WESHAKE  ** 

**  2  EXIT  WESHAKE  AND  RETURN  TO  DOS  ** 


****************************-k*-k**-k**i(*****-k*********** 

ENTER  ONE  OPTION  NUMBER: 


If  execution  Is  chosen  (Option  1)  the  program  begins  running.  The  status  of 
the  program  will  be  continually  shown  on  the  screen  along  with  times  of 
execution  for  each  option  and  notification  of  completion. 

81.  The  most  important  aspect  of  the  use  of  this  program  is  the 
evaluation  of  results.  The  results  are  contained  in  various  output  files  that 
are  automatically  created  as  a  consequence  of  using  the  various  mandatory 
actions  and  user  options.  The  various  output  files  were  described  previously 
(Table  III-5).  The  OUTPUT  file  corresponding  to  the  example  problem  is  shown 
in  Appendix  I . 

Print  menu 

82.  Once  all  operations  have  been  completed,  the  PRINT  MENU  will  be 
displayed.  The  output  files  created  and  summarized  in  Table  IV-1  may  be  sent 
directly  to  a  printer  on  the  (first)  parallel  port  (LPTl)  from  within  the 
WESHAKE  shell.  The  screen  for  the  PRINT  MENU  is: 


PRINT  MENU 
********** 

*************************************************** 

**  1  . 

PRINT  OUTPUT  FILE 

** 

**  2  - 

PRINT  STRESS/STRAIN  FILE 

** 

**  3  . 

PRINT  ACCELERATION  FILE 

** 

4  - 

PRINT  PUNCH  FILE 

** 

**  5  . 

PRINT  ACCELERATION  TIME  FILE 

**  6  - 

PRINT  VELOCITY  SPECTRUM  FILE 

** 

**  7  - 

PRINT  ACCELERATION  SPECTRUM  FILE 

** 

**  8  - 

FUTURE  USE 

** 

**  9  . 

FUTURE  USE 

** 

**10  - 

EXIT  WESHAKE 

** 

*************************************************** 

ENTER 

ONE  CHOICE: 

1 

(REPEAT  UNTIL  OPTION 

10  SELECTED) 
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Portions  of  the  OUTPUT  file  are  greater  than  80-characters  wide  so  it  is 
recommended  that  a  132 -character  printer  be  used  or  a  small  print  font. 


59 


0 


10 


f  20 

a. 

Ld 

Q 


30 


40 


V,  =  5,000  fps 


Figure  IV-2.  Example  problem  used  for  Part  IV 
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MANDATORY  ACTION  SOIL  PROFILE  FILE 


1.  Create  Soil  Column;  - 

1 — 

- > 

Option  2; 

Retrieve  Shear  and 

Soil  Column  Profile 

Damping  Curves  from 

- > 

Option  8: 

Data  Bases 

Modulus  and  Damping 

Curves 

a.  For  Option  3:  Creating  a  Soil  Profile  File 


MANDATORY  ACTIONS 


SPECIFICATION  FILE 


1.  Create  Soil  Column; 
Retrieve  Shear  and 
and  Damping  Curves 
from  Data  Bases 

2 .  ID  Earthquake  Motion 

-  User  Responses 

-  Earthquake 

Data  Base 

3.  Assign  Object  Motion 

-  User  Responses 


-> 

-> 

-> 

-> 


Option  2; 

Soil  Column  Profile 
Option  8: 

Modulus  and  Damping 
Curves 
Option  1: 

Earthquake  Motion 
Option  3: 

Object  Motion  and 
Type 

Option  4: 

Strain  Compatible 
Parameters  (DEFAULT) 


b.  For  Option  4:  Creating  a  Specification  File 


Figure  IV- 3.  Components  of  mandatory  actions  for 
SPECIFICATION  FILE  DESIGNATION  menu 
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PART  V:  ADVANCED  TOPICS 


83.  In  the  course  of  using  SHAKE,  engineers  have  found  ways  of 
obtaining  reasonable  and  meaningful  results  for  more  complicated  problems  than 
originally  intended  to  solve.  Some  of  these  cases  are  described  below  to 
enlighten  the  user  and  spawn  further  creativity. 

Application  of  Free-Field  Results 

84.  It  may  not  be  appropriate  to  directly  apply  the  free-field 
response  to  the  base  of  the  structure  for  a  niimber  of  reasons,  including; 

a.  The  depths  of  the  footings  most  likely  are  not  at  the  ground 
surface  and  motions  will  vary  with  depth. 

b.  The  weight  of  the  structure  acting  on  the  footings  will  affect  the 
motions  beneath  the  footings. 

c.  The  friction  acting  on  the  sides  of  the  footing  will  affect  the 
motions  acting  on  the  footing. 

d.  The  impedance  contrast  between  the  soil  and  foundation  is  normally 
quite  large. 

The  application  of  ground  motions  to  the  base  of  structures,  i.e.,  the 
consideration  of  points  such  as  those  listed,  is  commonly  referred  to  as 
dynamic  soil -structure  interaction  (DSSI)  and  can  be  considered  to  be 
distinctly  different  from  asynchronous  effects  produced  with  long  period 
waves . 

85.  Basic  design  approaches  for  dynamic  soil -structure  interaction 
have  recently  been  documented  by  Johnson  (1980)  and  Veletsos,  Prasad,  and  Tang 
(1988).  Evaluation  of  simple  foundation  systems  by  Veletsos,  Prasad,  and  Tang 
(1988)  suggests  the  following  rule  of  thumb:  for  a  range  of  lower  periods, 

DSSI  will  have  no  effect  on  the  response;  for  a  range  of  higher  periods,  DSSI 
will  reduce  the  maximum  response;  for  a  range  of  intermediate  periods,  DSSI 
might  increase  or  decrease  the  maximum  response. 

Estimation  of  Shear  Wave  Velocity 

86.  In  some  cases,  all  material  properties  necessary  for  the  use  of 
WESHAKE  may  not  be  available.  This  may  be  the  case  for  sites  with  very  deep 
soil  deposits  or  projects  with  limited  budgets.  There  is  no  substitute  for 
measured  properties  in  geotechnical  engineering  analysis,  in  this  case  the 


63 


primary  variable  being  shear  wave  velocity.  However,  it  may  be  possible  to 
evaluate  a  range  of  site  response  based  on  estimated  material  properties. 

87.  A  number  of  correlations  exist  In  the  literature  that  allow  the 
user  to  derive  shear  wave  velocities  (refer  to  summary  by  Sykora  and  Koester 
1988)  or  shear  modulus  coefficients  (Seed  et  al.  1986)  for  an  Intended  soil 
column.  However,  the  variations  among  proposed  correlations  between  shear 
wave  velocity  and  Standard  Penetration  Test  (SPT)  N- value  or  depth  are  large 
(Sykora,  1987a;  Sykora  and  Koester  1988).  The  correlation  between  (Kj)^^^,  and 
SPT  N-value  has  beevi  shown  to  produce  poor  results  (Sykora  1989).  The  user, 
therefore,  must  exercise  caution  and  good  engineering  judgement  when  using 
published  relationships  to  determine  shear  wave  velocity.  The  recommended 
procedure  is  to  use  a  data  base  of  shear  wave  velocities  (e.g.,  Sykora  1987b, 
Sykora  and  Koester  1991)  to  establish  a  reasonable  range  of  velocities  using 
known  site  conditions  and  compare  with  published  correlations. 

Studies  of  Modulus  and  Damning  Relationships 

88.  The  variation  of  material  properties,  namely  shear  modulus  and 
damping,  with  shear  strain  are  continually  under  study.  As  described  in 

Part  II,  the  original  investigations  for  this  phenomenon  were  conducted  in  the 
1960 's  and  early  1970's.  More  recently  in  the  1980 's  and  1990's,  studies  have 
been  completed  to  provide  more  information  on  the  behavior  of  fine-grained 
soils  and  gravels. 

89.  Studies  about  fi^e- grained  soils  by  investigators  at  different 
institutions  have  produced  options  for  different  sets  of  relationships.  For 
instance,  sets  of  shear  modulus  degradation  relationships  are  proposed  by  Zen 
and  Higuchi  (1984),  Sun,  Golesorkhi,  and  Seed  (1988),  and  Vucetic  and  Dobry 
(1991).  The  study  by  Vucetic  an(’  Dobry  is  intended  to  supersede  that  by  Sun, 
Golesorkhi,  and  Seed  in  that  new  data  were  combined  with  data  used  in  the  Sun 
study.  It  appears  as  though  the  new  modulus  curves  by  Vucetic  and  Dobry 
(1991)  are  indeed  more  comprehensive.  However,  the  new  damping  relationships 
are  at  issue.*  Apparently  some  of  the  new  damping  data  does  not  include 
large  strain  determinations  and  serves  to  unduly  weight  the  averaging 


*  Personal  communication,  Dr.  Joseph  Sun,  Woodward- Clyde  Consultants,  Inc., 
Oakland,  CA,  4  June  1992. 
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procedure.  Until  this  issue  is  sufficiently  resolved,  the  relationships  by 
Sun,  Golesorkhi,  and  Seed  (1988)  have  been  supplied  in  the  data  base. 


Hiph  Effective  Stresses 

90.  The  relationships  between  shear  modulus  and  shear  strain  for 
cohesionless  soils  may  be  significantly  influenced  by  the  effective  confining 
pressure  (Iwasaki,  Tatsuoka,  and  Takagi  1976).  The  influence  of  confining 
stress  on  cohesive  soils  is  not  conclusive  but  is  certainly  less  than  that  for 
cohesionless  soils  (Sun,  Golesorkhi,  and  Seed  1988). 

91.  A  simple  means  to  evaluate  the  potential  effect  of  high  confining 
stresses  is  to  select  a  modulus  curve  to  the  right  of  the  "best  estimate" 
curve.  A  rule  of  thumb  based  on  the  results  by  Seed  and  Idriss  (1970)  and 
Iwasaki,  Tatsuoka,  and  Takagi  (1976)  is  that  if  the  effective  vertical  stress 
at  the  center  of  the  layer  is  less  than  500  psf,  move  one  curve  to  the  left. 

If  the  effective  vertical  stress  is  between  500  and  2,000  psf,  use  the  best 
estimate  curve.  If  the  effective  vertical  stress  is  between  2,000  and  8,000 
psf,  move  over  one  curve  to  the  right.  If  the  effective  vertical  stress  is 
greater  than  8,000  psf,  move  over  two  curves  to  the  right. 

92 .  The  relationship  between  damping  ratio  and  shear  strain  may  also 
be  affected  by  confining  pressure  (Seed  et  al.  1986  and  Reeves  and  Castro 
1991).  This  finding  does  not  appear  to  be  applied  often  in  analyses  by  the 
profession,  however,  which  may  be  the  result  of  general  uncertainties  about 
normalized  relationships  for  damping  ratio.  Limitations  to  the  number  of  soil 
types  allowed  in  the  program  used  for  this  study  do  not  normally  facilitate 
involving  stress  adjusted  relationships  for  damping  ratio.  Potential 
variations  in  damping  ratio  from  stress  effects  ere  best  addressed  in  the 
sensitivity  analysis  for  damping  ratio  relationships. 

Multiple  Soil  Columns 

93.  Recent  standards  by  the  US  Nuclear  RegulatOi*y  Commission  (1989) 
and  American  Society  of  Civil  Engineers,  ASCE,  (1987)  provide  recommendations 
for  soil -structure  interaction  problems.  They  suggest  that  for  sites  that 
have  not  been  "well  investigated,"  a  representative  soil  column  should  be 
derived  for  the  site  of  interest  using  an  average  of  stratigraphy  and  measured 
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properties,  including  shear  wave  velocities.  The  site  response  should  then  be 
calculated  for  the  same  column  using  an  upper  and  lower  bound  of  shear  wave 
velocity  determined  by: 


Lower  bound: 


(G  _ 

max  (1+FACTOR) 


W 


Upper  bound:  (GBax)“*  =  FACTOR) 


(5) 


where 

ASCE  4-36:  FACTOR  -  coefficient  of  variation  not  to  be  less 
than  0 . 5 

NRC  3.7.2:  FACTOR  -  1.0 

to  account  for  potential  variations  of  material  properties. 

94.  These  standards  are  not  recommended  by  WES  for  a  site  response 
analysis  or  liquefaction  assessment  (particularly  the  use  of  a  single  column 
to  represent  a  site) .  The  use  of  lower  and  upper  bounds  to  determine  a 
potential  range  in  site  response  is  recommended  (refer  to  "Sensitivity 
Analysis"  below) .  There  are  several  reasons  for  this  recommendation, 
including; 

a .  The  "averaging"  of  stratigraphy,  including  the  total 
column  height,  across  a  heterogeneous  site  may  be  too 
subjective ; 

b ■  Past  experience  has  indicated  that  average  columns  may 
produce  unconservative  results  compared  to  the  range 
developed  with  the  collection  of  individual  columns ; 

c ■  The  standards  were  derived  primarily  for  soil -structure 
interaction  studies. 


In  the  end,  the  averaging  of  all  the  input  parameters  may  have  the  effect  of 
hiding  resonance  peaks  produced  by  layers  that  are  of  limited  extent  both 
laterally  and  vertically. 


Sensitivity  Analysis 

95.  The  sensitivity  of  various  inputs  to  WESHAKE  should  be  considered 
in  most  cases  of  site  response  analyses  to  evaluate  the  effects  of  potential 
variations  across  the  site  of  interest.  The  range  of  parameters  used  should 
be  a  function  of  measured  variations  or  perceived  uncertainties.  The  primary 
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inputs  to  be  considered  are;  depth  to  bedrock,  shear  modulus  and  damping  ratio 
curves,  and  maximum  shear  modulus  (function  of  V,  or  unit  weight  of 
soil) . 


Vertical  Response 

96.  SHAKE  has  also  been  used  by  some  investigators  to  estimate  the 
vertical  response  of  sites  to  earthquake  excitations.  The  method  used  for 
this  purpose  only  evaluates  the  vertical  propagation  of  compression  waves;  it 
does  not  include  the  effect  of  shear  waves.  Details  of  the  method  are 
presented  below. 

97.  Vertical  response  can  be  estimated  using  SHAKE  by  taking  steps  to 
match  strains  (and  therefore  matching  percentage  of  modulus  reduction  and 
damping  increase)  from  calculations  for  vertical  response  with  calculations 
for  horizontal  response.*  This  may  be  accomplished  by  using  some  arbitrary 
value  of  the  ratio  of  effective  strain,  PRMUL.  Normally  for  the  calculation 
of  horizontal  response,  PRMUL  is  0.65  (65  percent).  The  strains  calculated 
for  the  vertical  response  and  initial  PRMUL  are  then  compared  with  the  strains 
calculated  for  horizontal  response.  *The  procedure  is  repeated  by  varying 
PRMUL  until  the  variation  of  strain  with  depth  for  vertical  response  matched 
those  for  horizontal  response  as  closely  as  possible.  Values  of  PRMUL  most 
likely  will  be  greater  than  unity,  and  could  be  as  large  as  50.  A  factor 
greater  that  unity  means  that  the  constrained  modulus  degrades  much  faster 
than  the  shear  modulus.  The  selection  of  an  appropriate  value  of  effective 
strain  involves  subjective  decision  making.  Close  matches  should  not  be 
expected,  especially  for  near-surface  layers. 

98.  Inherent  in  the  above  procedure  is  the  assumption  that  the 
variation  of  normalized  constrained  modulus  is  similar  to  that  for  normalized 
shear  modulus.  At  large  strains,  this  assumption  may  lead  to  considerable 
errors.  Some  laboratory  tests  conducted  by  WES  for  a  particular  project  in 
the  late  1970 's  included  torsional  and  longitudinal  vibration  testing  (Curro 
and  Marcuson  1978).  The  results  of  these  tests  indicate  that  at  longitudinal 
strains  between  about  10'*  and  10'^  percent,  the  maximum  modulus  decreased 
rapidly  by  about  50  percent.  This  characteristic  "break"  in  the  normalized 


Personal  communication,  Prof.  John  Lysmer,  University  of  California  at 
Berkeley,  1  November  1990. 
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curve  differs  greatly  from  the  characteristic  smooth  shape  and  maximum  slopes 
of  the  normalized  shear  modulus  curves.  Therefore,  large  differences  are 
expected  between  calculated  vertical  response  and  actual  site  response. 

99.  For  vertical  response  calculations,  the  constitutive  model  for 
stiffness  is  assumed  to  be  the  same  except  that  the  maximum  constrained 
moduli ,  ,  is  used,  calculated  from  compression  waves; 

^  V|  (6) 

™  g 


where 

Vp  -  compression  wave  velocity 

100.  The  comparison  of  calculations  of  motions  using  SHAKE  with 
measured  vertical  response  have  not  been  reported  in  the  literature. 
Consequently,  validation  of  SHAKE  for  vertical  response  is  considered  to  be 
inadequate.  The  calculation  of  vertical  motions  using  this  program  are 
considered  to  merely  provide  qualitative  insight  into  vertical  site  behavior. 

Other  Uses 


101.  Elton,  Shie,  and  Hadj-Hamou  (1991)  showed  that  the  shear  stresses 
calculated  using  SHAKE  and  a  computer  program  used  to  calculate  the  two- 
dimensional  response,  FLUSH  (Lysmer  et  al.  1973),  were  in  good  agreement  for 
both  level-ground  sites  and  sloped  sites  except  near  the  surface.  For 
embankments,  the  agreement  improved  as  the  period  of  excitation  moved  farther 
from  the  natural  period  of  the  embankment.  Comparisons  made  of  other  ground 
motion  parameters,  however,  such  as  the  variation  of  acceleration  with  time 
and  velocity  or  acceleration  response  spectrum  showed  much  greater 
variation.* 


*  Personal  communication.  Prof.  David  Elton,  Auburn  University,  Auburn,  AL, 
May  1992. 
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Computer  Hardware  and  Software  Requirements 


Al.  The  computer  program  will  run  on  an  IBM  compatible  personal  computer, 
including  an  XT  model  using  an  8086  processor.  The  size  and  power  of  the 
computer  will  serve  to  better  enhance  the  performance  of  the  program.  The  speed 
of  execution  is  compared  for  an  example  problem  with  15  layers  and  a  few  output 
options  in  Table  A-1.  These  times  may  be  used  as  a  relative  basis  for 
comparison. 

A2 .  Some  minimum  requirements  exist  to  successfully  execute  WESHAKE .  They 
include:  465  Kbyte  RAM  and  version  3.1  of  DOS.  (The  DOS  3.1  (and  higher) 
commands  mem  and  ver  will  determine  the  memory  available  and  version, 
respectively.)  If  DOS  5.0  or  greater  is  installed,  some  lower  system  memory  can 
be  transferred  to  upper  memory  thus  freeing  up  more  RAM.  The  math  coprocessor 
option  is  mandatory  for  computers  with  8086,  80286,  and  80386  processors.  No 
graphics  requirements  exist  for  WESHAKE  version  1.0. 

Installation 


A3.  All  necessary  components  of  WESHAKE  are  contained  in  one  executable 
file,  WESHAKE.EXE,  and  the  shear  modulus,  damping  ratio  and  earthquake  data  bases 
(SHEARDB,  DAMPDB,  and  EARTHQ,  respectively)  that  can  be  copied  to  a  directory  on 
a  computer  hard  drive  or  run  directly  from  the  floppy  disk  drive.  It  is 
recommended  that  a  directory  be  created  on  a  hard  disk  specifically  for  this 
program  file  in  addition  to  various  files  that  may  be  created  by  the  user  as  a 
consequence  of  running  this  program.  To  create  a  directory,  for  example  WESHAKE, 
and  load  the  software  onto  a  hard  disk,  begin  at  the  start-up  prompt,  (C:\>), 
then  type : 

C:\>  MD  WESHAKE 

Recall  that  DOS  commands  are  not  case  sensitive.  Now  change  directories  by 
typing: 

C:\>  CD  WESHAKE 

If  the  command  PROMPT  $p$g  is  in  the  AUTOEXEC.BAT  file  of  your  computer,  the 
computer  prompt  will  be: 

C : \WESHAKE> 


A2 


Read  the  files  on  the  distribution  disk  by  inserting  the  WESHAKE  floppy  disk  into 
a  floppy  drive,  say  drive  A,  and  typing: 

C:\WESHAKE>  COPY  A:*.* 

To  begin  WESHAKE,  t3rpe: 

C:\WESHAKE>  WESHAKE 

The  use  of  WESHAKE  is  described  in  the  main  text  of  this  report. 


Limitations 


A4.  Certain  limitations  exist  within  the  program  to  minimize  the  overall 
size  of  the  program.  Limitations  can  be  modified  by  WES  to  accomodate  the  needs 
of  the  users.  The  current  list  of  more  important  limitations  are  summarized  in 
Table  A- 2.  Less  important  limitations,  t3rpically  those  pertaining  to  output 
options,  are  specified  in  the  option  menus  shown  in  Part  IV  of  the  main  text. 


Run-Time  Optimization 


A5 .  Run  times  are  highly  machine  dependent  as  might  be  expected  and  shown 
above.  For  a  given  machine,  the  run  times  can  also  vary  considerably  depending 
on  the  values  of  certain  input  parameters.  Some  of  the  more  time-consuming 
parameters  are:  the  number  of  FFT  terms,  the  number  of  layers,  and  the 
convergence  criterion  (ERR),  The  number  of  FFT  terms  and  number  of  layers  may 
non-negotiable ,  ERR  is  easily  adjusted;  as  ERR  is  decreased,  the  number  of 
iterations  required  increases  (and  therefore  time  to  solution  increases).  An 
optimal  value  of  ERR  is  probably  about  5  percent,  the  value  used  by  default  in 
WESHAKE. 
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Table  A-1 

Performance  Statistics  for  WESHAKE 


Example  File  for  comparison  with  a  15  layer  system 


Computer  Brand 
and  model 

Processor 

Processor 
Speed  (MHz) 

Execution 
Time  (Sec) 

Ratio 

Time  XT 

IBM  PC -XT 

Intel  8086 

8 

1374 

1 

IBM  PC -AT 

Intel 

In-Board  386 

16 

209 

.15 

Compaq  Deskpro 

Intel  80386 

16 

258 

.19 

Compaq  Deskpro 

Intel  80386 

20 

193 

.14 

IBM  PS/2,  Model  70 

Intel  80386 

20 

206 

.15 

Unisys 

Intel  80386 

20 

160 

.12 

Compu  Add 

Intel  80386 

25 

161 

.12 

Gateway  2^'00 

Intel  80486 

25 

41 

.03 

Gateway  2000 

Intel  80486 

33 

21 

.01 
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Table  A- 2 

Array  Limits  for  WESHAKE 


Description 

Program 

Variable 

Limit 

Maximum  number  of 
(including  rock) 

layers  in  soil  column 

ML 

20 

Maximum  number 
relationships 

of 

modulus  degradation/damping 

NSOILT 

10 

Maximum  number 
(FFT)* 

of 

terms  in  Fast  Fourier  Transform 

MAMAX 

4096 

Number  of  soil 

layers  for  acceleration  output 

15 

*  MAMAX  must  be  aX  greater  than  or  equal  to  the  largest  FFT  used  in  the 
specifications  file  and  at  least  twice  as  large  as  number  of  non-zer^'  terms  in 
the  earthquake  records.  The  array  in  SHAKE!  has  been  fixed  for  a  maximum  of  4096 
points  to  greatly  reduce  fixed  memory  size. 
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APPENDIX  B: 

SHEAR  MODULUS  AND  DAMPING  RATIO  DATA  BASES 


The  modulus  and  damping  data  bases  have  the  following  format: 

Columns  Format  Parameter  (s') 

FIRST  LINE 

1-5  15  NUMBER  OF  VALUES  PLOTTED  ON  RVE  [NV] 

-  maximum  of  20  points  per  curve 

6- 10  15  MULTIPLICATION  FACTOR  IN  PLOTTING  [NPL] 

7- 15  15  CURVE  IDENTIFICATION  NUMBER  [NUM]^ 

16  -  75  A60  IDENTIFICATION  OF  SOIL  PROFILE  [ID] 

SECOND  LINE(S) 

1  -  80  8F10.4  NV  VALUES  OF  SHEAR  STRAIN  (percent)  IN  INCREASING  ORDER 

[R]: 

-  Eight  per  line 

The  second  line  repeats  until  all  values  of  shear  strain  have  been  specified 
THIRD  LINE(S) 

1-80  8  F10.4  NV  VALUES  OF  NORMALIZED  SHEAR  MODULUS  (percent)  OR 

DAMPING  RATIO  (percent),  IN  INCREASING  ORDER, 
CORRRES PONDING  TO  VALUES  OF  SHEAR  STRAIN  IN  SECOND 
LINE(S)  [U]: 

-  Eight  per  line 

The  third  line  repeats  until  all  values  of  modulus  or  damping  have  been 
specified. 

The  second  and  thi.  1  lines  form  a  set  and  a  set  is  created  for  each  material  type 
(i.e.,  these  lines  are  repeated)  in  the  respective  data  bases  (unique  to  NUM) . 


^  Deviation  from  syntax  of  specification  file  and  OPTION  8  of  SHAKE 
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Shear  Modulus  Data  Base 


8  100. 

1  ROCK 

(Schnabel 

1973) 

0.0001 

0.0003 

0.0010 

0.0030 

0.0100 

0.0300 

0.1000 

1.00 

1.00 

0.99 

0.95 

0.90 

0.81 

0.725 

9  100. 

2  GRAVEL.  Average 

(Seed  et  al.  1986) 

0.0001 

0.0003 

o.riio 

0.0030 

0.0100 

0.0300 

0.1000 

1.000 

1.00 

0.050 

0.97 

0.87 

0.73 

0.55 

0.37 

0.20 

9  100. 

3  SAND, 

Lower  Bound  (Seed  & 

Idriss  1970) 

0.0001 

0.0003 

0.0010 

0.0030 

0.0100 

0.0300 

0.1000 

1.0000 

1.000 

0.04 

0.985 

0.93 

0.83 

0.635 

0.425 

0.225 

9  100. 

4  SAND, 

Average  (Seed  &  Idriss  1970) 

0.0001 

0.0003 

0.0010 

0.0030 

0.0100 

0.0300 

0.1000 

1.0000 

1.00 

0.06 

0.98 

0.95 

0.89 

0.73 

0.52 

0.29 

9  100. 

5  SAND, 

Upper  Bound  (Seed  & 

Idriss  1970) 

0.0001 

0.0003 

0.0010 

0.0030 

0.0100 

0.0300 

0.1000 

1.0000 

1.00 

0.085 

1.00 

0.99 

0.96 

0.85 

0.655 

0.37 

9  100. 

6  CLAY 

(Pl-5-10, 

Sun  et  al. 

1988) 

0.0001 

1.0 

0.0003 

0.0010 

0.0030 

0.0100 

0.0300 

0.1000 

1.00 

1.00 

0.975 

0.91 

0.78 

0.565 

0.305 

0.04 

9  100. 

7  CLAY 

(PI-10-20, 

Sun  et  al 

.  1988) 

0.0001 

0.0 

0.0003 

0.0010 

0.0030 

0.0100 

0.0300 

0.1000 

1.00 

1.00 

1.00 

0.96 

0.87 

0.70 

0.41 

0.08 

9  100. 

8  CLAY 

(Pl-20-40, 

Sun  et  al 

.  1988) 

0.0001 

1.0 

0.0003 

0.0010 

0.0030 

0.0100 

0.0300 

0.1000 

1.00 

1.00 

1.00 

0.97 

0.90 

0.77 

0.52 

0.14 

9  100. 

9  CLAY 

(PI-40-80, 

Sun  et  al 

.  1988) 

0.0001 

i.O 

0.0003 

0.0010 

0.0030 

0.0100 

0.0300 

0.1000 

1.00 

0.20 

1.00 

1.00 

0.985 

0.92 

0.815 

0.62 

9  100. 

10  CLAY 

(PI>80,  Sun  et  al .  1988) 

0.0001 

1.0 

0.0003 

0.0010 

0.0030 

0.0100 

0.0300 

0.1000 

1.00 

0.33 

1.00 

1.00 

0.985 

0.94 

0.86 

0.71 

9  100. 

11  Mexico  City  Clay  (Sun  et  ^ 

al.  1988) 

0.0001 

1.0 

0.0003 

0.0010 

0.0030 

0.0100 

0.0300 

0.1000 

1.00 

0.46 

1.00 

1.00 

1.000 

0.995 

0.975 

0.920 

1.0000 

0.55 

0.3000 

0.10 

0.3000 

0.11 

0.3000 

0.14 

0.3000 

0.19 

0.3000 

0.14 

0.3000 

0.20 

0.3000 

0.30 

0.3000 

0.41 

0.3000 

0.53 

0.3000 

0.8000 
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5  5.0  1  ROCK  (Schnabel  1973) 

0.0001  0.0010  0.0100  0.1000  1.0000 

0.40  0.80  1.50  3.00  4.60 


9  5.0 

0.0001 

2  GRAVEL  (Average,  Seed  et  al  1986) 

0.0003  0.0010  0.0030  0.0100  0.0300 

0.1000 

0.3000 

1.0000 

0.8 

1.0 

1.9  3.0  5.4  9.6 

15.4 

20.8 

24.6 

9  5.0 

0.0001 

3  SAND 
0.0003 

(Average,  Seed  &  Idriss  1970) 

0.0010  0.0030  0.0100  0.0300 

0.1000 

0.3000 

1 . 0000 
0.8 

1.0 

1.9  3.0  5.4  9.6 

15.4 

20.8 

24.6 

9  5.0 

0.0001 

4  SAND 
0.0003 

(Lower  Bound,  Seed  &  Idriss  1970) 

0.0010  0.0030  0.0100  0.0300 

0 . 1000 

0.2780 

1.0000 

0.3 

0.4 

0.7  1.4  2.7  5.0 

9.8 

15.5 

20.7 

9  5.0 

0.0001 

5  SAND 
0.0003 

(Upper  Bound,  Seed  &  Idriss  1970) 

0.0010  0.0030  0.0100  0.0315 

0.1000 

0.3000 

1.0000 

0.7 

1.2 

2.7  5.5  9.9  14.8 

21.0 

25.5 

27.9 

9  5.0 

0.0001 

6  CLAY 
0.0003 

(Average,  Seed  &  Idriss  1970) 

0.0010  0.0030  0.0100  0.0300 

0.1000 

0.4000 

1.0000 

2.5 

2.5 

2.5  3.5  4.5  6.5 

9.0 

13.5 

20.5 

9  5.0 

0.0001 

7  CLAY 
0.0003 

(Lower  Bound,  Seed  Idriss  1970) 

0.0010  0.0030  0.0100  0.0300 

0.1000 

0.4000 

1.0000 

1.3 

1.3 

1.3  1.5  1.7  3.5 

4.0 

6.5 

12.3 

9  5.0 

0.0001 

8  CLAY 
0.0003 

(Upper  Bound,  Seed  &  Idriss  1970) 

0.0010  0.0030  0.0100  0.0300 

0.0780 

0.3000 

1.0000 

4.0 

4.0 

4.0  5.0  7.5  11.0 

16.0 

21.8 

APPENDIX  C: 
EARTHQUAKE  DATA  BASE 


Cl.  The  earthquake  data  base  (EARTHQ)  is  described  and  characteristics  of 
each  earthquake  motion  is  presented.  Note  that  earthquake  records  to  be  used 
with  WESHAKE  that  are  not  included  in  the  data  base  need  to  comply  with  the 
specification  file  format  (refer  to  Appendix  E)  and  not  the  data  base  format. 
The  differences  between  these  two  formats  can  be  seen  by  comparing  the  format 
below  with  that  presented  beginning  on  page  E4.  The  earthquake  data  base  has  the 
following  format: 


Columns  Format 

Parameter (s) 

FIRST  LINE 

1  - 

5 

15 

NUMBER  OF  VALUES  IN  EARTHQUAKE  MOTION 

[NV] 

6  - 

10 

15 

NUMBER  OF  TERMS  IN  FFT  [MA] 

-  must  be  a  power  of  2  and  <  MAMAX 

-  should  be  >  2  *  NV 

11  - 

20 

F10.3 

TIME  STEP  FOR  MEASUREMENT  [DT] 

21  - 

25 

15 

EARTHQUAKE  DATA  BASE  NUMBER  [INUME]^ 

26  - 

85 

A60 

EARTHQUAKE  TITLE  [EQTITLE] 

SECOND  LINE^ 

1  - 

80 

A80 

CHARACTERISTICS  OF  EARTHQUAKE 

THIRD  LINE 

1  - 

10 

F10.3 

MULTIPLICATION  FACTOR  FOR  ACCELERATION 

[XF] 

11  - 

20 

F10.3 

MAXIMUM  ACCELERATION  VALUE  TO  BE  USED 

[XMAX] 

21  - 

30 

F10.3 

MAXIMUM  (CUTOFF)  FREQUENCY  [FMAX] 

FOURTH  AND  SUBSEQUENT  LINES 

1  -  72  8(1X,F8.6)  NV  VALUES  OF  ACCELERATION  (g's)  [XR] 

-  eight  per  line 

73  -  79  17  LINE  NUMBER  [K] 


C2 .  Note  that  the  number  of  non-zero  values  in  the  earthquake  record  must 
be  less  than  or  equal  to  half  of  the  maximum  number  of  terms  for  the  FFT  (MA) . 
The  records  in  the  data  base  all  conform  to  this  requirement.  Users  providing 
additional  records  must  ensure  that  this  requirement  is  met,  however.  The  proper 
format  is  shown  in  Appendix  C  along  with  the  existing  data  base  and  plots  of 
records.  Also,  the  maximum  number  of  terms  for  the  FFT  can  not  exceed  4096. 
This  value  was  fixed  in  WESHAKE  and  can  be  modified  is  necessary. 


^  Deviation  from  syntax  of  specification  file  and  OPTION  1  of  SHAKE 
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C3.  At  the  present  time,  there  are  22  earthquakes  in  the  data  base. 
Important  information  about  earthquakes  in  the  data  base  was  presented  in  Tables 
III-l  and  III-2  of  the  main  text.  The  data  for  these  records  are  presented  on 
the  following  pages  in  the  form  of  figures.  Numerical  values  can  be  extracted 
from  the  data  base.  The  first  figure  for  each  earthquake  is  a  plot  of  the 
variation  of  acceleration  with  time.  The  second  figure  for  each  earhtquake  is 
a  plot  of  pseudo -response  velocity  spectra  using  a  tripartite  format  (assuming 
5  percent  damping)  at  6  levels  of  system  damping  (2,  5,  7,  10,  12  and  15 
percent) . 
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Figure  Cl.  Golden  Gate  record  of  1957  San  Francisco,  California,  aftershock 


Pseudo  Relative  Response  Velocity,  (in. /sec) 


Frequency,  f  (Hz) 

10'  10’  10°  10”’ 
[TT  I  I  1 — I \ - 1 - |"l  M  I  I~l  I - 1 - 1  M  1  I  I  I  I  I  I 


Period,  T  (sec) 


Figure  C2 .  Tripartite  presentation  of  pseudo-velocity  spectra  for 
Golden  Gate  record  of  1957  San  Francisco,  California,  aftershock 
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Figure  C3.  Cholame -Temblor  record  of  1966  Parkfield,  California,  earthquake 


Pseudo  Relative  Response  Velocity,  S  ,,  (in./sec) 


Frequency,  f  (Hz) 

10'  10'  10°  10“’ 
|Ti"i  rr  I — I 1 - 1 1  I  M  I  I — I 1 - [Mill  I — I 1 - 1 


10“'  10“’  10°  io’ 

Period,  T  (sec) 


Figure  C4.  Tripartite  presentation  of  pseudo-velocity  spectra  for 
Cholame -Temblor  record  of  1966  Parkfield,  California,  earthquake 
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Pseudo  Relative  Response  Velocity,  S  ^  (in./sec) 


Magnitude  =  6.5  Rock  Outcrop 


O 
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Figure  C7 .  Lake  Hughes  Array  No.  4  record  of  1971  San  Fernando,  California,  earthquake 


Pseudo  Relative  Response  Velocity,  S  ,,  (in. /sec) 


Frequency,  f  (Hz) 

10"  id  10°  10“' 

P'n  I  I — I  I  I - 1 1 1  I  I  I — I — 1 1 - 1 1  I'l  I  1 — I — I 1 - 1 
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Figure  C8.  Tripartite  presentation  of  pseudo -velocity  spectra  for 
for  Lake  Hughes  Array  No.  4  record  of  1971  San  Fernando, 
California,  earthquake 
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Figure  C9.  Sitka  Magnetic  Observatory  record  of  1972  Alaskan  earthquake 


Pseudo  Relative  Response  Velocity,  S  ,,  (in./sec) 


Frequency,  f  (Hz) 

10"  10'  10°  10“' 
I  I  I  I  I  I — I - 1 - 1 - 1  I  I  I  I  I  I - 1 - 1 - 1  T  M  I  1 — I - 1 - 1  I 


Period,  T  (sec) 


Figure  CIO.  Tripartite  presentation  of  pseudo-velocity  spectra  for 
Sitka  Magnetic  Observatory  record  of  1972  Alaskan  earthquake 
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Figure  Cll.  Gilroy  No.  1  record  of  1974  Hollister,  California,  earthquake 


Pseudo  Relative  Response  Velocity,  S  ^  Cin./sec) 


Frequency,  f  (Hz) 

lO'  id  10°  10"' 

1 1  I  IT  r  I — I 1 - pTi  I  I  1 — I - 1 - I'l  M  1  I  I  I  :  I 


10  ■  10 
Period,  T  (sec) 


Figure  C12 .  Tripartite  presentation  of  pseudo-velocity  spectra  for 
Gilroy  No.  1  record  of  1974  Hollister,  California,  earthquake 
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record  of  1979  Coyote  Lake,  California,  earthquake 


Pseudo  Relative  Response  Velocity,  S  ^  (in./sec) 


10^ 

ri  rr  I  I  I — r 


Frequency,  f  (Hz) 

10’  10°  10"’ 
T  I  I  1  U  I  I - 1 - 1 - 1  1  M  I  I  1 - 1 - 1 - 1 


10""  10“’  10°  io’ 

Period,  T  (sec) 


Figure  C14.  Tripartite  presentation  of  pseudo-velocity  spectra  for 
Gilroy  No,  1  record  of  1979  Coyote  Lake,  California,  earthquake 
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Figure  C15.  Superstition  Mountain  record  of  1979  Imperial  Valley,  California,  earthquake 


Pseudo  Relative  Response  Velocity,  (in./sec) 


Frequency,  f  (Hz) 

lO"  10’  10°  10“’ 

I M  I  I  I  I  I  I  1 1  n  I  I  I  I  I  1 1  I  I  T  "i  I  I  I  I 


io“"  10“'  10°  id 
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Figure  C16.  Tripartite  presentation  of  pseudo-velocity  spectra  for 
Superstition  Mountain  record  of  1979  Imperial  Valley, 
California,  earthquake 


o 

o’ 


d  6  d  c)  Q  o  d  d  d 

'll 


(5)  D  ‘uOipj0p3DV 


C20 


Figure  C17.  Superstition  Mountain  record  of  1981  El  Centro,  California,  earthquake 


Pseudo  Relative  Response  Velocity,  S  ^  (in./sec) 


(-requency,  f  (Hz) 

10"  id  10°  10“' 
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Figure  C18.  Tripartite  presentation  of  pseudo-velocity  spectra  for 
Superstition  Mountain  record  of  1981  El  Centro, 

California,  earthquake 
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Figure  C19.  Gilroy  No.  1  record  of  1984  Morgan  Hill,  California,  earthquake 


Pseudo  Relative  Response  Velocity,  S  ^  (in./sec) 


Frequency,  f  (Hz) 

10"  id  10°  10'’ 
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Figure  C20.  Tripartite  presentation  of  pseudo-velocity  spectra  for 
Gilroy  No.  1  record  of  1984  Morgan  Hill,  California,  earthquake 
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Figure  C21.  Iverson  site  record  of  1985  Nahannl,  Northwest  Territories,  aftershock 


Pseudo  Relative  Response  Velocity,  S„  (in. /sec) 
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Figure  C23.  Slide  Mountain  record  of  1985  Nahanni,  Northwest  Territories,  aftershock 


Pseudo  Relative  Response  Velocity,  (in./sec) 


Frequency,  f  (Hz) 
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Figure  C2A.  Tripartite  presentation  of  pseudo-velocity  spectra  for 
Slide  Mountain  record  of  1985  Nahanni,  Northwest 
Territories,  aftershock 
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Hollister  Airport  record  of  1989  Loma  Prieta.  California,  earthquake 


Pseudo  Relative  Response  Velocity,  S  ^  (in./sec) 


Frequency,  f  (Hz) 

10"  id  10°  10“’ 
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Figure  C26.  Tripartite  presentation  of  pseudo-velocity  spectra  for 
Hollister  Airport  record  of  1989  Loma  Prieta, 

California,  earthquake 
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Figure  C27.  Synthetic  MCE  "Record  A"  used  for  Folsom  Dam 
seismic  stability  study  (Bolt  and  Seed  1983) 


Pseudo  Relative  Response  Velocity,  (in./sec) 
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Figure  C28.  Tripartite  presentation  of  pseudo-velocity 
spectra  for  synthetic  MCE  "Record  A"  used  for  Folsom 
Dam  seismic  stability  study  (Bolt  and  Seed  1983) 
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Figure  C29.  Synthetic  MCE  "Record  B"  used  for  Folsom  Dam 
seismic  stability  study  (Bolt  and  Seed  1983) 


Pseudo  Relative  Response  Velocity,  S  ,,  (in. /sec) 


Frequency,  f  (FIz) 

lO"  id  10°  10"’ 
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Figure  C30.  Tripartite  presentation  of  pseudo -velocity 
spectra  for  synthetic  MCE  "Record  B"  used  for  Folsom 
Dam  seismic  stability  study  (Bolt  and  Seed  1983) 
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Figure  C31.  Synthetic  500-year  event  record,  horizontal  1  component, 
for  New  Madrid,  Missouri,  earthquake  (Risk  Engineering,  Inc.  1992) 
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Figure  C32.  Tripartite  presentation  of  pseudo-velocity  spectra  for 
synthetic  500-year  event  record,  horizontal  1  component,  for  New 
Madrid,  Missouri,  earthquake  (Risk  Engineering,  Inc.  1992) 
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3.  Synthetic  500 -year  event  record,  horizontal  2  component, 
Madrid,  Missouri,  earthquake  (Risk  Engineering,  Inc.  1992) 


Pseudo  Relative  Response  Velocity,  S  „  (in. /sec) 


Frequency,  f  (Hz) 
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Figure  C34.  Tripartite  presentation  of  pseudo-velocity  spectra  for 
synthetic  500-year  event  record,  horizontal  2  component,  for  New 
Madrid,  Missouri,  earthquake  (Risk  Engineering,  Inc.  1992) 
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Figure  C35.  Synthetic  1000-year  event  record,  horizontal  2  component, 
for  New  Madrid,  Missouri,  earthquake  (Risk  Engineering,  Inc.  1992) 
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Figure  C36.  Tripartite  presentation  of  pseudo -velocity  spectra  for 
synthetic  1000-year  event  record,  horizontal  1  component,  for  New 
Madrid,  Missouri,  earthquake  (Risk  Engineering,  Inc.  1992) 


C39 


o 


tD 

o 

in 

o 

ID 

O 

r-; 

iT) 

OJ 

q 

fN 

r^ 

6 

d 

d 

d 

6 

d 

d 

I  I  I 


(6)  D  ‘uOj^DJ0iaDDV 


C40 


Figure  C37.  Synthetic  1000-year  event  record,  horizontal  2  component, 
for  New  Madrid,  Missouri,  earthquake  (Risk  Engineering,  Inc.  1992) 
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Figure  C38 .  Tripartite  presentation  of  pseudo -velocity  spectra  for 
synthetic  1000-year  event  record,  horizontal  2  component,  for  New 
Madrid,  Missouri,  earthquake  (Risk  Engineering,  Inc.  1992) 
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Figure  C40.  Tripartite  presentation  of  pseudo -velocity  spectra  for 
synthetic  5000-year  event  record,  horizontal  1  component,  for  New 
Madrid,  Missouri,  earthquake  (Risk  Engineering,  Inc.  1992) 
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Figure  C41.  Synthetic  5000-year  event  record,  horizontal  2  component 
for  New  Madrid,  Missouri,  earthquake  (Risk  Engineering,  Inc.  1992) 
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Figure  C42.  Tripartite  presentation  of  pseudo -velocity  spectra  for 
synthetic  5000-year  event  record,  horizontal  2  component,  for  New 
Madrid,  Missouri,  earthquake  (Risk  Engineering,  Inc.  1992) 
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Figure  C44.  Tripartite  presentation  of  pseudo -velocity 
spectra  for  synthetic  MCE  record  used  for 
Ririe  Dam  Stability  Study  (Seed  1987) 


APPENDIX  D: 


MOTION:  ACCELEROGRAM  PLOTTING  PROGRAM 


Dl.  The  computer  plotting  program  called  MOTION  can  be  used  to  visualize 
earthquake  motions  that  are  used  in  WESHAKE  program.  Currently,  earthquake 
motions  must  be  stored  in  the  earthquake  data  base  called  "EARTHQ".  If  the  user 
ha.«;  access  to  better  plotting  software,  WESHAKE  still  prints  the  earthquake 
motion  to  an  output  file  called  "EQIN"  which  can  be  used  as  in  input  to  their 
software  package  (some  modifications  to  the  input  file  might  have  to  be  generated 
by  the  user  for  the  various  software  packages) .  As  the  earthquake  data  base  is 
updated  each  year,  these  updates  will  be  automatically  incorporated  into  the 
MOTION  software  package. 

D2 .  MOTION  involves  two  basic  stages;  identification  and  plotting.  The 
identification  stage  involves  the  selection  from  a  menu  sample  earthquake  motions 
(identical  with  the  selection  of  motions  in  WESHAKE) .  MOTION  allows  the  user  to 
pick  as  many  many  of  the  sample  earthquake  in  the  data  base  that  will  be  required 
for  the  analaysis.  MOTION  will  allow  the  user  to  continually  pick  each  earthquake 
motion.  MOTION  will  also  generate  a  print  out  of  the  graph.  This  option  is 
limited  to  an  Epson  printer  or  an  emulation  of  an  Epson  printer. 

D3 .  In  order  to  execute  MOTION,  the  user  types  MOTION  at  the  C:\  prompt. 
The  first  screen  following  the  execution  of  MOTION  is  (next  page) : 


******************** 

WESHAKE 

******************** 

EARTHQUAKE  PLOT  MENU 
******************** 

WRITTEN  BY: 

DR.  DAVID  C.  WALLACE 
APPLIED  COMPUTER  SCIENCE  DEPARTMENT 
ILLINOIS  STATE  UNIVERSITY 
NORMAL,  ILLINOIS  61761 

*********************************************************** 

THE  COMPUTER  WILL  GENERATE  A  PLOT  OF  YOUR  EARTHQUAKE  MOTION 
*********************************************************** 

WHEN  YOU  ARE  FINISHED  VIEWING  THE  DISPLAY, 

CAN  HIT  THE  ENTER  KEY  TO  EXIT  THE  PROGRAM. 
*********************************************************** 


HIT  ANY  KEY  TO  CONTINUE 


The  next  screen  will  be  the  selection  screen.  The  selection  screen  lists  the 
various  earthquake  motions  which  can  be  plotted  on  the  screen. 
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EARTHQUAKE  DATA  BASE 

********************************************************* 


* 

NO. 

MEASURED  RECORD 

NO. 

SYNTHETIC 

RECORD 

* 

* 

* 

1 

GOLDEN  GATE  1957 

14 

FOLSOM  RECORD  "A" 

* 

* 

2 

PARKFIELD  1966 

15 

FOLSOM  RECORD  "B" 

* 

* 

3 

CASTAIC  RIDGE  1971 

16 

NEW  MADRID 

500-YR  HI 

* 

* 

4 

LAKE  HUGHES  #  4  1971 

17 

NEW  MADRID 

500-YR  H2 

* 

* 

5 

SITKA  1972 

18 

NEW  MADRID 

1000 -YR 

HI 

* 

* 

6 

GILROY  #1  1974 

19 

NEW  MADRID 

1000-YR 

H2 

* 

* 

7 

GILROY  #1  1979 

20 

NEW  MADRID 

5000-YR 

HI 

* 

* 

8 

SUPERSTITION  1979 

21 

NEW  MADRID 

5000-YR 

H2 

* 

* 

9 

SUPERSTITION  1981 

22 

RIRIE  DAM 

* 

* 

10 

GILROY  #1  1984 

* 

* 

11 

IVERSON  1985 

* 

* 

12 

SLIDE  MT  1985 

* 

* 

13 

HOLLISTER  AIRPORT  1989 

* 

******************************************************** 


ENTER  EARTHQUAKE  MOTION  NUMBER: 

1 


Once  the  earthquake  motion  is  selected,  the  user  can  choose  to  generate  a  print 
out  of  the  graph  on  the  screen.  The  following  screen  prompts  the  user  for 
print  out  of  the  plotted  graph. 


******************************************************** 
**  DO  YOU  WANT  A  PRINT  OUT  OF  THE  EARTHQUAKE  MOTION  ?  ** 
**  YOU  NEED  AN  EPSON  PRINTER  OR  AN  EMULATION  OF  AN  ** 
**  EPSON  PRINTER  FOR  BEST  RESULTS.  ** 

******************************************************** 
ENTER  0  FOR  NO 
ENTER  1  FOR  YES 

0 


Once  the  user  responds  to  the  last  screen,  the  program  will  generate  the  graph 
of  the  earthquake  motion  on  the  screen.  If  the  user  enters  1  for  a  print  out 
of  the  motion,  the  screen  graph  of  the  earthquake  motion  will  be  basically  two 
colors  which  are  easily  generated  for  the  EPSON  print.  If  the  user  enters  0  for 
no  print  outs,  the  graph  should  be  a  muti-color  graph.  Figure  D1  shows  the  type 
of  graph  generated  by  MOTION.  The  message  in  the  lower  right  hand  corner  of  the 
graph  tells  the  user  to  hit  the  enter  key  to  continue. 
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Figure  Dl.  Example  plot  from  HOTIOH 


The  computer  will  then  prompt; 


DO  YOU  WANT  TO  PLOT  MORE  EARTHQUAKE  MOTIONS? 
ENTER  0  FOR  NO 
ENTER  1  FOR  YES 

1 


If  the  user  responds  yes  to  the  prompt,  the  program  will  start  over  again  by 
generating  the  menu  for  earthquake  motions  again. 
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APPENDIX  E; 

SPECIFICATION  FILE  FORMAT 


I.  PROJECT  TILE 


FIRST  LINE 

1-5  15  MAXIMUM  NUMBER  OF  TERMS  FOR  FOURIER  TRANSFORM  [MAMAX] 

6  -  65  A60  PROJECT  TITLE  [PTITLE]^ 


II.  MANDATORY  ACTIONS 

A.  READ  SOIL  COLUMN  (OPTION  2  IN  SHAKE) 
FIRST  LINE 

1-5  15  OPTION  NUMBER  [KK]  -  2 

SECOND  LINE 


1  - 

5 

15 

SOIL  COLUMN  NUMBER  [MSOIL] 

6  - 

10 

15 

NUMBER  OF  SOIL  LAYERS  [ML] 

11  - 

15 

15 

NUMBER  OF  SOIL  LAYER  ABOVE  WHICH 
[MWL] 

IS  THE  WATER  TABLE 

16  - 

25 

FIO.O 

UNIT  WEIGHT  OF  PORE  LIQUID  (kef) 

[WW] 

26  - 

61 

6A6 

IDENTIFICATION  FOR  THE  SOIL  PROFILE  [IDNT] 

THIRD  LINE(S)^ 

1  - 

5 

15 

LAYER  NUMBER  [K] 

6  - 

10 

15 

TYPE  OF  SOIL  LAYER  [TP] 

-  Corresponds  layer  with  TP'th 

set  of  material 

properties  defined  in  the  next  option 
(SHAKE  OPTION  8) 


11 

-  15 

15 

NUMBER  OF  SUBLAYERS  IN  LAYER  K  [NLN] 

16 

-  25 

FIO.O 

THICKNESS  OF  LAYER  (ft)  [H] 

26 

-  35 

FIO.O 

COEFFICIENT  OF  LATERAL  EARTH  PRESSURE 

[SKO] 

36 

-  45 

FIO.O 

DAMPING  RATIO  (percent)  [BL] 

46 

-  55 

FIO.O 

MOIST  UNIT  WEIGHT  FOR  LAYER  (kef)  [W] 

56 

-  65 

FIO.O 

INITIAL  SHEAR  WAVE  VELOCITY  OR  Kj  [RVl]^ 

66 

-  75 

FIO.O 

MAXIMUM  SHEAR  WAVE  VELOCITY  OR  (Kj)^ 

[RV2] 

76 

-  80 

F5.0 

DAMPING  MODIFICATION  FACTOR  [BF] 

81 

11 

CODE  INDICATING  INPUT  TYPE  [ INKEY] 

-  0  -  SHEAR  WAVE. 

-  1  - 

The  third  line  is  repeated  for  each  soil  layer,  including  rock.  For  rock,  the 
parameters:  H,  SKO,  and  INKEY  are  not  applicable. 
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Deviation  from  SHAKE 


B.  READ  MATERIAL  PROPERTIES  (OPTION  8  IN  SHAKE) 

FIRST  LINE 

1-5  15  OPTION  NUMBER  [KK]  -  8 

SECOND  LINE 

1-5  15  NUMBER  OF  SOIL  TYPES  FOR  PLOTTING  [NST] ; 

-  Maximum  of  4 

6-10  15  PLOT  OPTION  [NPL] : 

-  0:  No  plot 

-  1:  Plot  in  OUTPUT 

11  -  15  15  NUMBER  OF  STRAIN  UNITS  IN  EACH  LOG  CYCLE  [NPL] 

16-25  FIO.O  MAXIMUM  VALUE  OF  ORDIANTE  [SC] : 

-  0:  Maximum  value  of  data 

THIRD  LINE 

1-5  15  NUMBER  OF  VALUES  PLOTTED  ON  RVE  [NV] 

-  Maximum  of  20  points  per  curve 

6-10  F5.0  MULTIPLICATION  FACTOR  IN  PLOTTING  [FPL] 

12  -  77  11A6  DESCRIPTION  OF  SOIL  PROFILE  [ID] 

FOURTH  LINE(S) 

1  -  80  8F10.3  NV  VALUES  OF  SHEAR  STRAIN  (percent)  IN  INCREASING 

ORDER  [S] 

-  Eight  per  line 

The  fourth  line  repeats  until  all  values  of  shear  strain  have  been  specified 
FIFTH  LINE(S) 

1  -  80  8F10.3  NV  VALUES  OF  NORMALIZED  SHEAR  MODULUS  (percent)  OR 

DAMPING  RATIO  (percent),  IN  INCREASING  ORDER, 
CORRRESPONDING  TO  VALUES  OF  SHEAR  STRAIN  IN  SECOND 
LINE(S)  [Y] 

-  Eight  per  line 

The  fifth  line  repeats  until  all  values  of  modulus  or  damping  have  been 
specified. 

The  fourth  and  fifth  lines  form  a  set  and  a  set  is  created  for  each  material 
type  (i.e.,  these  lines  are  repeated)  in  the  respective  data  bases  (unique  to 
NUM)  . 
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C.  SELECT  EARTHQUAKE  RECORD  (OPTION  1  IN  SHAKE) 


FIRST 

LINE 

1  - 

5 

15 

OPTION  NUMBER  (KK)  -  1 

SECOND  LINE 

1  - 

5 

15 

NUMBER  OF  VALUES  IN  EARTHQUAKE  MOTION 

[NV] 

6  - 

10 

15 

NUMBER  OF  TERMS  IN  FFT  [MA] : 

-  Must  be  a  power  of  2  and  <  MAMAX 

-  Should  be  >  2  *  NV 

11  - 

20 

F10.3 

TIME  STEP  FOR  MEASUREMENT  [DT] 

22  - 

50 

5A6 

EARTHQUAKE  TITLE  [EQTITLE] 

THIRD 

LINE^ 

1  - 

10 

FIO.O 

MULTIPLICATION  FACTOR  FOR  ACCELERATION 

[XF] 

11  - 

20 

FIO.O 

MAXIMUM  ACCELERATION  VALUE  TO  BE  USED 

[XMAX] 

21  - 

30 

FIO.O 

MAXIMUM  (CUTOFF)  FREQUENCY  [FMAX] 

31  - 

35 

15 

ECHO  CODE  [OUTKEY]:^ 

-  0:  Do  not  echo  accelerations  in  OUTPUT 

-  1;  Echo  accelerations  in  OUTPUT 


FOURTH  AND  SUBSEQUENT  LINES 

1  -  72  8(1X,F8.6)  NV  VALUES  OF  ACCELERATION  (g's)  [XR] : 

-  Eight  per  line 

73  -  79  17  LINE  NUMBER  [K] 


D.  ASSIGN  OBJECT  MOTION  (OPTION  3  IN  SHAKE) 


FIRST  LINE 

1-5  15  OPTION  NUMBER  {KKJ  -  3 

SECOND  LINE 


1-5  15 

6-10  15 


NUMBER  OF  (SUB)LAYER  WHERE  OBJECT  MOTION  IS  ASSIGNED 
[IN] 

TYPE  OF  (SUB) LAYER  (INT) : 

-  0:  Outcropping 

-  1:  (Sub) layer  within  profile 


EU 


E.  OBTAIN  STRAIN  COMPATIBLE  PROPERTIES  (OPTION  4  IN  SHAKE) 
FIRST  LINE 

1-5  15  OPTION  NUMBER  [KK]  -  4 

SECOND  LINE 


1  - 

5 

15 

PUNCH  OPTION  [KS] : 

-  0:  Do  not  write  to  PUNCH  file 

-  1:  Write  to  PUNCH  file 

6  - 

10 

15 

MAXIMUM  NUMBER  OF  ITERATIONS  [ITMAX] 

11  - 

20 

FIO.O 

MAXIMUM  ACCEPTABLE  DIFFERENCE  BETWEEN  THE  LAST -US ED 
MODULUS  AND  DAMPING  VALUES  AND  THE  STRAIN  COMPATIBLE 
VALUES  (percent)  [ERR] 

21  - 

30 

FIO.O 

RATIO  BETWEEN  EFFECTIVE  STRAIN  AND  MAXIMUM  STRAIN 
[PRMUL] : 

0.65  reconunended 


III.  OPTIONAL  ACTIONS 
A.  COMPUTE  MOTION  IN  SPECIFIED  LAYERS 
FIRST  LINE 

1-5  15  OPTION  NUMBER  [KK]  -  5 

SECOND  LINE 

1  -  80  1515  (SUB)LAYERS  FOR  COMPUTATION  OF  MOTION  [LL5]: 

•  Maximum  15 

THIRD  LINE 

1  -  80  1515  TYPE  OF  MOTION  CORRESPONDING  TO  LL5  [LT5] : 

-  Maximum  15 

-  0 ;  Outcropping 

-  1:  (Sub) layer  within  soil  profile 

FOURTH  LINE 

1  -  80  1515  OUPUT  OPTION  CORRESPONDING  TO  LL5  and  LT5  [LP5]: 

-  0:  Maximum  accelerations  only  to  OUTPUT  and  ACCEL 

-  1;  Acceleration  time  history  to  PUNCH 
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B.  PRINT  OR  PUNCH  OBJECT  MOTION 
FIRST  LINE 

1-5  15  OPTION  NUMBER  [KK]  -  6 

SECOND  LINE 

1-5  15  SELECT  MODE  OF  OUTPUT  [K2]: 

-  0:  Maximum  acceleration  only 

-  1:  Write  to  PUNCH  file 

-  2;  Write  to  PUNCH  and  OUTPUT  files 


C.  CHANGE  OBJECT  MOTION 
FIRST  LINE 


1  - 

5 

15 

OPTION  NUMBER  [KK]  -  7 

SECOND  LINE 

1  - 

5 

15 

NUMBER  OF  SUBLAYER  [LLl] : 

-  0:  OBJECT  MOTION  ORIGINALLY  ASSIGNED  IS  RETAINED 

6  - 

10 

15 

TYPE  OF  (SUB) LAYER  [LTl]; 

11  - 

20 

FIO.O 

-  0;  Outcropping 

-  1:  (Sub) layer  within  soil  profile 

MULTIPLICATION  FACTOR  FOR  ACCELERATION  VALUES  [XF] : 

21- 

30 

FIO.O 

-  1.0:  No  change 

NEW  TIME  STEP  [DTNEW] 

D. 

COMPUTE  RESPONSE 

SPECTRA 

FIRST 

LINE 

1  - 

5 

15 

OPTION  NUMBER  [KK]  -  9 

SECOND  LINE 

1  - 

5 

15 

SUBLAYER  NUMBER  [LLl]: 

-  0:  Object  Motion 

6  - 

10 

15 

TYPE  OF  SUBLAYER  [LTl]: 

-  0:  Outcropping 

-  1;  (Sub) layer  within  soil  profile 
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THIRD  LINE 


1  - 


15 


NUMBER  OF  DAMPING  VALUES  TO  BE  USED  [ND] : 
-  Maximum  of  6 


6  - 

10 

15 

PUNCH  OPTION  [KP]: 

-  0:  No  write  to  PUNCH  file 

-  1:  Write  to  PUNCH  file 

11  - 

15 

15 

PARAMETER  OPTION  [KAV] ; 

-  0:  Spectral  velocity 

-  1 ;  Spectral  acceleration 

-  2:  Spectral  velocity  and  acceleration 

16  - 

20 

15 

PLOT  OPTION  [KPL] ; 

-  0:  Store  plot  for  later  (combined)  plotting 

-  1:  Plots  of  all  spectra  calculated  to  this  point 

21  - 

25 

15 

SITE  PERIODS  FOR  COMPUTATIONS  [KPER] : 

-  KPER  -  0 

9  LINEAR  STEPS  from  0,1  to  1.0  sec 


5 

steps 

from 

1.0 

to 

2.0 

sec 

4 

steps 

from 

2.0 

to 

4.0 

sec 

KPER 

-  1 

18 

steps 

from 

0.1 

to 

1.0 

sec 

10 

steps 

from 

1.0 

to 

2.0 

sec 

8 

steps 

from 

2.0 

to 

4.0 

sec 

KPER 

-  2 

38 

steps 

from 

0.05 

to 

1.0 

sec 

20 

steps 

from 

1.0 

to 

2.0 

sec 

30 

steps 

from 

2.0 

to 

5.0 

sec 

KPER 

-  3 

FOURTH  LINE 
1  -  60  6F10.0 


LOG  INCREMENTS  with  10  steps  per  log  unit  from 
0.05  TO  5.0 
-  KPER  -  4 

LOG  INCREMENTS  with  10  steps  per  log  unit  from 
0.05  TO  10.0 


ND  VALUES  OF  CRITICAL  DAMPING  RATIOS  (decimal)  [ZLD] 


E.  INCREASE  TIME  INTERVAL 
FIRST  LINE 

1-5  15  OPTION  NUMBER  [KK]  -  10 

SECOND  LINE 

1-5  15  FACTOR  FOR  INCREASING  TIME  INTERVAL  [IFR]: 

-  Must  be  a  power  of  2 


E7 


F.  DECREASE  TIME  INTERVAL 


FIRST  LINE 

1-5  15  OPTION  NUMBER  [KK]  -  11 

SECOND  LINE 

1-5  15  FACTOR  FOR  DECREASING  TIME  INTERVAL  [IFR]; 

-  Must  be  a  power  of  2 


G.  CALCULATE  FOURIER  SPECTRUM  OF  OBJECT  MOTION 
FIRST  LINE 

1-5  15  OPTION  NUMBER  [KK]  -  12 

SECOND  LINE 

1-5  15  PLOTTING  OPTION  [Kl]: 

-  0:  Store  for  later  (combined)  plot 

-  1;  Plot  all  stored  spectra 

6-10  15  NUMBER  OF  TIMES  THE  SPECTRUM  IS  TO  BE  SMOOTHED  [NSW] 

11  -  15  15  NUMBER  OF  VALUES  TO  BE  PLOTTED  [N] : 

-  Maximum  of  2049 


H.  CALCULATE  FOURIER  SPECTRUM  OF  COMPUTED  MOTION 
FIRST  LINE 

1-5  15  OPTION  NUMBER  [KK]  -  13 

SECOND  LINE 

1-5  15  (SUB)LAYER  NUMBER  [LL] 

6-10  15  TYPE  OF  (SUB)LAYER  [LT] 

-  0:  Outcropping 

-  1:  (Sub)layer  within  soil  profile 

11  -  15  15  PLOT  OPTION  [LP]: 

-  0;  Store  for  later  (combined)  plot 

-  1:  Plot  all  stored  spectra 

16  -  20  15  NUMBER  OF  TIMES  THE  SPECTRUM  IS  TO  BE  SMOOTHED  [LNSW] 

21  -  25  15  NUMBER  OF  VALUES  TO  BE  PLOTTED  [LLL] : 

-  Maximum  of  2049 
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I.  PLOT  TIME  HISTORY  OF  OBJECT  MOTION 
FIRST  LINE 

1-5  15  OPTION  NUMBER  [KK]  -  14 

SECOND  LINE 

1-5  15  NUMBER  OF  VALUES  SKIPPED  IN  PLOTTING  [NSKIP] : 

-  Every  NSKIP  values  skipped 

6-10  15  NUMBER  OF  VALUES  TO  BE  PLOTTED  [NN] : 

-  Maximum  of  2049 


J.  COMPUTE  AMPLIFICATION  SPECTRUM 
FIRST  LINE 


1-5  15  OPTION  NUMBER  [KK]  -  15 

SECOND  LINE 


1  - 

5 

15 

NUMBER  OF  FIRST  LAYER  [LIN] 

6  - 

10 

15 

FIRST  LAYER  TYPE  [LINT): 

-  0:  Outcropping 

-  1:  (Sub)layer  within  soil  profile 

11  - 

15 

15 

NUMBER  OF  SECOND  LAYER  [LOUT] : 

16  - 

20 

15 

SECOND  LAYER  TYPE  [LOTP]: 

-  0:  Outcropping 

-  1:  (Sub)layer  within  soil  profile 

21  - 

25 

15 

PLOTTING  OPTION  [KP] : 

-  0:  Store  for  later  plotting 

-  1:  Plot  all  stored  data 

26  - 

30 

F5.0 

NUMBER  OF  FREQUENCY  STEPS  [DFA] : 

-  Amplification  factor  is  computed  for  first  200 
frequencies  at  interval  DFA  (Hz)  beginning  at  0 

32  - 

78 

8A6 

DESCRIPTION  [IDAMP] 
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K.  COMPUTE  STRESS  OR  STRAIN  HISTORY  IN  MIDDLE  SUBLAYERS 


FIRST  LINE 

1-5  15  OPTION  NUMBER  [KK]  -  16 

SECOND  LIME 


1  - 

5 

15 

FIRST  (SUB) LAYER  NUMBER  [LLL] 

6  - 

10 

15 

SELECT  TYPE  OF  RESPONSE  [LLGS]: 

-  0:  STRAIN 

-  1:  STRESS 

11  - 

15 

15 

PUNCH  OPTION  [LLPCH]: 

-  0:  No  write  to  PUNCH  file 

-  1:  Write  to  PUNCH  file 

16  - 

20 

15 

PLOT  OPTION  [LLPL]: 

-  0:  No  plot  in  OUTPUT 

-  1:  Plot  in  OUTPUT 

21  - 

25 

15 

NUMBER  OF  VALUES  TO  BE  PLOTTED  [LNV] : 

-  Maximum  of  2049 

26  - 

35 

FIO.O 

SCALE  FOR  PLOTTING  [SK] : 

-  (i.e.,  maximum  value  of  ordinate) 

-  0:  maximum  of  data 

37  - 

65 

THIRD 

5A6 

LINE 

DESCRIPTION  [ID] 

1  - 

5 

15 

SECOND  (SUB) LAYER  NUMBER  [LLL] 

6  - 

10 

15 

SELECT  TYPE  OF  RESPONSE  [LLGS]: 

-  0:  STRAIN 

-  1:  STRESS 

11  - 

15 

15 

PUNCH  OPTION  [LLPCH]: 

-  0:  No  write  to  PUNCH  file 

-  1:  Write  to  PUNCH  file 

16  - 

20 

15 

PLOT  OPTION  [LLPL]: 

-  0:  No  plot  in  OUTPUT 

-  1:  Plot  in  OUTPUT 

21  - 

25 

15 

NUMBER  OF  VALUES  TO  BE  PLOTTED  [LNV]; 
-  Maximum  of  2049 

26  - 

35 

FIO.O 

SCALE  FOR  PLOTTING  [SK]: 

-  (i.e.,  maximum  value  of  ordinate) 

-  0;  maximum  of  data 

37  - 

65 

5A6 

DESCRIPTION  [ID] 

NOTE:  LEAVE  THIRD  LINE  BLANK  IF  ONLY  ONE  RESPONSE  IS  TO  BE  COMPUTED 

IV.  END  OF  INPUT  FILE 

1-5  15  OPTION  NUMBER  [KK]  -  0 
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APPENDIX  F: 
VALIDATION  OF  WESHAKE 


1024  VALIDATION  PROBLEM:  Example  problem  from  Schnabel,  Lysmer,  Seed  (1972) 

2  OPTION  2:  READ  SOIL  COLUMN  DATA 


0 

9 

9 

.06240 

Example 

1 

2 

1 

7.0 

.43 

.050 

.120 

61.3 

43.1.0001 

2 

1 

1 

13.0 

.43 

.100 

.100 

234. 

430.1.0000 

3 

1 

1 

10.0 

.43 

.030 

.100 

567. 

745.1.0000 

4 

1 

1 

12.0 

.43 

.030 

.100 

567. 

920.1.0000 

5 

2 

1 

20.0 

.43 

.050 

.125 

33.6 

76.1.0001 

6 

1 

1 

18.0 

.43 

.050 

.125 

507. 

910.1.0000 

7 

1 

1 

20.0 

.43 

.030 

.125 

717. 

1090.1.0000 

8 

1 

1 

20.0 

.43 

.030 

.123 

802. 

1155.1.0000 

9 

3 

1 

.030 

.150 

8000. 

8000.1.0000 

8  OPTION  8:  READ  MATERIAL  PROPERTIES 
3  0  2  100.  STRAIN  COMPATIBLE  PROPERTIES 


9100.0 

CLAY  (Schnabel,  Lysmer, 

,  4  Seed 

1972)  Shear 

Modulus 

.0001 

.000316  .0010 

.00316 

.0100 

.0316 

.1000 

.3160 

1.0000 

1 . 0000 

.9130  .7610 

.3630 

.4000 

.2610 

.1320 

.0760 

.0370 

8  3.0 

CLAY  (Schnabel,  Lysmer, 

1  4  Seed 

1972)  Damping 

.0001 

.0010  .00316 

.0100 

.0316 

.1000 

.3160 

1.0000 

2.0000 

2.3000  3.5000 

4 . 7300 

6.3000 

9.2500 

13.7500 

20.0000 

9100.0 

SAND  (Schnabel,  Lysmer, 

,  4  Seed 

1972)  Shear 

Modulus 

.0001 

.000316  .0010 

.00316 

.0100 

.3316 

.1000 

.3160 

1.0000 

1.0000 

.9840  .9340 

.8260 

.6360 

.4430 

.2460 

.1150 

.0490 

8  3.0  SAND  (Schnabal,  Lysoicr  &  S«*d  1972)  Dan^lnc 

.0001  .0010  .0030  .0100  .0300  .1000  .3000  1.0000 

.8000  1.6000  3.1200  S.8000  9. 3000  13.4000  20.9000  23.0000 

8100.0  ROCK  (Schiubal,  Lysoiar,  (  S««d  1972)  Shear  Modulus 
.0001  .0003  .0010  .0030  .0100  .0300  .1000  l.OQOO 

1.0000  1.0000  .9880  .9330  .9000  .8100  .7230  .3300 

3  3.0  ROCK  (Schnabel,  Lysmer,  (  Seed  1972)  OaiiN>lns 

.0001  .0010  .0100  .1000  1.0000 

.4000  .8000  1.3000  3.0000  4.6000 

1  OPTION  1:  READ  EARTHQUAKE  RECORD 
800  1024  .020  PASADENA  1932 

1.000  0.02  23.  1 


002493  - 

.003313 

-.002623 

-.002326 

-.001317 

-.000128 

.001274 

.002382 

1 

002625 

.002518 

. 002390 

.002313 

. 002369 

. 002302 

.  002551 

.002491 

2 

002333 

.002478 

.002316 

.002463 

.002500 

.002451 

.002486 

.002408 

3 

001994 

.001425 

.000923 

.000443 

.000423 

.000420 

.000066  - 

.000829 

4 

001525  - 

.002386 

-.002381 

-.001979 

-.001647 

-.001658  - 

.001639  - 

.001643 

5 

001099  - 

.000310 

.000131 

-.000088 

-.000763 

-.001374  - 

.002079  - 

.002893 

6 

003379  - 

.003281 

-.002718 

-.002233 

-.001509 

.000499 

.002231 

.002673 

7 

003059 

.003649 

.004454 

.005218 

.003970 

.006832 

.  006573 

.004815 

8 

003340 

.000761 

-.001550 

-.001932 

-.001298 

-.000392 

.000267 

.000936 

9 

001844 

.001362 

.000454 

-.000321 

-.001442 

-.002423  - 

.003165  - 

.003689 

10 

004324  ' 

.003168 

- . 003828 

-.005713 

-.003336 

-.005068  - 

.003841  - 

.002303 

11 

001126  - 

.000120 

.001399 

. 002900 

.004300 

.004992 

.004919 

.004081 

12 

002831 

.001863 

.001481 

.001031 

.001121 

.001382 

.001816 

.001063 

13 

001154  - 

.003149 

-.005312 

-.006364 

-.007083 

-.007836  - 

.007772  - 

.006703 

14 
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.005669 

-.004535 

-.003444 

-.002324 

-.001229 

-.000185 

.000867 

.002097 

15 

.003343 

.004508 

.004360 

.003868 

.004611 

. 005548 

.006484 

.007403 

16 

.008347 

.009253 

.010469 

.012118 

.012267 

.011130 

.010349 

.008871 

17 

.006792 

.005192 

.004711 

.003808 

.005427 

.011079 

.011“94 

.012237 

18 

.012588 

.013780 

.015396 

.016909 

.018496 

.020026 

.021497 

.021543 

19 

.021695 

.019953 

.015352 

.013395 

.013656 

.013217 

.009669 

.005633 

20 

.004387 

.003030 

.002845 

.003635 

.004411 

.004443 

.004606 

.003431 

.000276 

-.001939 

-.001932 

-.006423 

-.009774 

-.011147 

- . 012292 

-.012660 

22 

.013513 

-.014858 

-.017037 

-.018811 

-.020927 

-.022735 

-.024635 

-.025974 

23 

.027901 

-.026475 

-.019656 

-.014708 

-.007272 

.002294 

.003546 

.006332 

24 

.004549 

.003636 

.002636 

.002806 

.001202 

-.000432 

- . 002572 

-.001484 

Zj 

.000616 

.000291 

-.000486 

-.006645 

-.014722 

-.012778 

-.033575 

-.042831 

26 

.043009 

-.045807 

-.043294 

-.041669 

-.037745 

- . 025278 

-.019201 

.001040 

27 

.013648 

.012866 

.015824 

.016378 

.018696 

. 019668 

.021342 

.020955 

28 

.021309 

.022849 

.026435 

.029246 

.031983 

.033228 

.034739 

.031762 

29 

.027738 

.024080 

.020339 

.017203 

.011562 

.007591 

-.015690 

-.028448 

30 

.029791 

-.038844 

-.042214 

-.043695 

-.040997 

-.038664 

-.033649 

-.030902 

31 

020213 

-.001768 

.002023 

.007762 

.006125 

.004710 

. 007523 

.009468 

32 

.010966 

.011913 

.012652 

.020084 

.028702 

.032657 

.038837 

.040949 

33 

041798 

.038336 

.032960 

.031175 

.029425 

.026144 

.025072 

.022677 

34 

021590 

.016643 

.017778 

.015500 

.019219 

-.005005 

-.018653 

-.035037 

35 

.047254 

-.042567 

-.045124 

-.039513 

-.034399 

-.020552 

-.004749 

-.002274 

36 

000997 

-.001284 

-.000624 

-.000064 

.000558 

.000123 

-.001028 

-.000921 

37 

000648 

.000488 

.003457 

.006631 

.007229 

.005943 

.005730 

.005723 

38 

005276 

.003131 

.003588 

.004385 

.001724 

-.000949 

.000235 

.003089 

39 

004174 

.002623 

.003442 

.003680 

.004803 

.005154 

.005346 

.007711 

40 

010554 

.009195 

.008490 

.005691 

.005005 

.000425 

-.003943 

-.006070 

41 

007428 

-.009936 

-.010439 

-.012590 

-.011388 

-.019622 

-.039587 

-.041574 

42 

046392 

-.048356 

-.052906 

-.052793 

-.049356 

-.045835 

-.043181 

-.027839 

43 

014535 

-.013859 

.025816 

.043108 

.043491 

. 047994 

.047682 

.051035 

44 

051602 

.054637 

.054828 

.057238 

. 055020 

.049601 

.043950 

.040517 

45 

037062 

.034520 

.031703 

.024568 

.014474 

.011676 

.004249 

-.007235 

46 

026973 

-.041443 

-.042865 

-.045273 

-.045146 

-.044333 

-.040249 

-.033307 

47 

017015 

-.004583 

-.004544 

-.003638 

-.004924 

-.004706 

-.005664 

-.005841 

48 

007728 

-.008074 

-.010774 

-.012404 

-.014434 

-.009161 

-.007826 

.003369 

49 

024767 

.025711 

.029998 

.032436 

.034883 

.034808 

.035693 

.035950 

50 

031252 

.022314 

.016231 

.008574 

-.003059 

-.007785 

-.009396 

-.009446 

51 

008794 

-.008981 

-.009000 

-.009324 

-.009375 

-.009651 

-.009853 

-.011662 

52 

013760 

-.015836 

-.017792 

-.018124 

-.018907 

-.016268 

-.010307 

-.003454 

53 

002631 

.004109 

.004371 

.001234 

-.001186 

-.005053 

- . 008826 

-.012269 

54 

013740 

-.015331 

-.017610 

-.022831 

-.024390 

- . 026574 

-.024601 

-.022072 

55 

018882 

-.015641 

-.011956 

-.012234 

-.014489 

-.016093 

-.017898 

-.021668 

56 

027007 

-.027043 

-.017194 

-.007552 

.001476 

.017401 

.023309 

.023177 

57 

027040 

.028719 

.032583 

.031626 

.030859 

.028315 

.027878 

.023949 

58 

023372 

.013997 

-.008389 

-.017724 

-.021337 

-.033693 

-.037941 

-.036294 

59 

036374 

- .035359 

-.034194 

-.030404 

-.027147 

-.023444 

-.020092 

-.016374 

60 

012404 

- .008319 

-.004306 

-.000756 

.002788 

. 006267 

.010007 

.013599 

61 

017748 

.020897 

.018017 

.012843 

.010150 

.006306 

.004310 

-.002774 

62 

011569 

-.013108 

-.012718 

-.009093 

-.004293 

.000696 

.001663 

.003167 

63 

003833 

.005218 

.005565 

.005307 

.004555 

.004324 

. 002268 

-.002006 

64 

005485 

-.008136 

-.010942 

-.013972 

-.017321 

- . 020956 

- . 023507 

-.019832 

65 

016115 

-.014444 

-.012891 

-.011308 

- . 009298 

-.007287 

-.006066 

-.004843 

66 

004426 

-.005197 

-.005685 

- . 006885 

-.009521 

-.010996 

-.009946 

-.009179 

67 

007934 

-.005983 

-.004401 

-.005072 

- . 005589 

-.006038 

-.002838 

.000617 

68 
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.003857 

.006366 

.007463 

.005531 

.004975 

.002051 

-.001933 

-.008681 

69 

.011511 

-.011291 

-.011264 

-.011423 

-.011180 

-.012156 

-.013166 

-.01*751 

70 

.015533 

-.016972 

-.018909 

-.024914 

-.024367 

-.021703 

-.017011 

-.016185 

71 

.0079*1 

.017308 

.015199 

.017725 

.017247 

. 018997 

.  018<;77 

.019598 

72 

.018284 

.018230 

.017224 

.017375 

.016698 

.015711 

.010446 

. 008182 

73 

.009113 

.010604 

.011417 

.011690 

.011490 

.011663 

.011506 

.011620 

7* 

.011817 

.012871 

.013767 

.013843 

.009655 

. 005485 

.001648 

-.001211 

75 

.005169 

-.009233 

-.013879 

-.018128 

-.015787 

-.010113 

-.006007 

-.002833 

76 

.000194 

.002856 

.005626 

.008930 

.011888 

.014931 

.016787 

.018965 

77 

.020808 

.022943 

.025658 

. 026834 

.023618 

.021884 

. 018231 

.015899 

78 

.006415 

-.003205 

-.005577 

-.014988 

-.018955 

-.021957 

-.024742 

-.027639 

79 

. 030063 

- . 028838 

-.026851 

-.025145 

-.026489 

-.029640 

-.032318 

-.035157 

80 

.035179 

-.034898 

-.033905 

-.032181 

-.030237 

-.028657 

-.025223 

-.019982 

81 

.014925 

- . 009185 

-.003261 

.001808 

.005464 

.009307 

.012679 

.014940 

82 

.017126 

.019247 

.020912 

.022423 

.024020 

.025809 

.029569 

.031829 

83 

.030313 

.028309 

. 026858 

.022232 

.013520 

.010163 

.006899 

.003967 

84 

.000597 

-.000495 

.000120 

.000476 

.000989 

.001265 

.000967 

.000524 

85 

.000241 

-.000063 

-.000264 

-.000584 

.000541 

.002474 

.004356 

.005954 

86 

.007249 

.008512 

.009877 

.010134 

.009723 

.009107 

.007355 

.005186 

87 

.004931 

.006719 

.008156 

.009798 

.012428 

.016216 

.018755 

.016836 

88 

.  014459 

.011738 

. 008580 

.005413 

.002914 

.001680 

.000421 

-.000937 

89 

.002750 

-.004659 

-.006680 

-.012302 

-.015390 

-.015706 

-.016604 

-.017175 

90 

.017957 

-.018356 

-.014196 

-.010422 

-.009050 

-.006920 

-.005341 

-.002633 

91 

.000002 

.003054 

.004365 

.004209 

.004327 

.003693 

.001851 

.000085 

92 

.001814 

-.003043 

-.003466 

-.004068 

-.004504 

-.003832 

-.001625 

-.000235 

93 

.002278 

.005874 

.014722 

.025498 

.026123 

.028565 

.029956 

.031018 

94 

.029316 

.028463 

.026981 

.026664 

.025729 

.025371 

.024097 

.023324 

95 

.021865 

.021047 

.018546 

.014894 

.012006 

.010343 

.008053 

.005314 

96 

.001757 

-.002601 

-.007363 

-.010712 

-.0154f9 

-.021659 

-.021587 

-.019652 

97 

.018365 

-.016570 

-.015225 

-.013493 

-.011212 

-.007042 

-.006514 

-.007858 

98 

.009230 

-.010467 

-.012825 

-.015684 

-.016584 

-.010786 

-.005233 

.000579 

99 

.005328 

.009647 

.012139 

.015901 

.018338 

. 021552 

.019896 

.021718 

100 

3  OPTION  3;  ASSIGN  OBJECT  MOTION 

9  0 

h  OPTION  4;  OBTAIN  STRAIN -COMPATIBLE  PROPERTIES 

0  10  5.00  .65 

5  OPTION  5;  ACCELERATION  RECORDS 

1234567899 
1111111110 
0000000000 

9  OPTION  9:  RESPONSE  SPECTRUM 

1  0 

10  2  11 
0.05 

0  OPTION  0;  END  OF  INPUT 
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VALIDATION  PROBLEM:  Example  problem  from  Schnabel,  Lysmer 


MAX.  NUMBER  OF  TERMS  IN  FOURIER  TRANSFORM  -  1024 

NECESSARY  LENGTH  OF  BLANK  COMMON  X  -  6419 


1******  OPTION  2  ***  READ  SOIL  PROFILE 


MSOIL  -  0 

ML  -  9 

MWL  -  9 

WW  -  .0624 

IDNT  -  Example 

NEW  SOIL  PROFILE  NO.  0  IDENTIFICATION  -  Example 


SHEAR/K2  FACTOR  WAVE  VELOCITY  INPUT  BY  LAYER 


NUMBER  OF  LAYERS 

9 

DEPTH 

TO 

BEDROCK 

120.00 

NUMBER  OF  FIRST  SUBMERGED 

LAYER 

9 

DEPTH 

TO 

WATER  LEVEL 

120.00 

UNIT  WEIGHT  OF  WATER  - 

.0624 

kef 

F5 


SHEAR  MODULUS  SAND  SQ.ROOT  REL.  7.0  .0  7.0 

DAMPING  SAND 

SHEAR  MODULUS  CLAY  13.0  7.0  20.0 

DAMPING  CLAY 
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1****** 


1  ***  READ  INPUT  MOTION 


EARTHQUAKE  -  PASADENA  1952 

800  ACCELERATION  VALUES  AT  TIME  INTERVAL  .0200 


THE  VALUES  ARE  LISTED  ROW  BY  ROW  AS  READ  FROM  CARDS 
TRAILING  ZEROS  ARE  ADDED  TO  GIVE  A  TOTAL  OF  1024  VALUES 


-.00249 

-.00331 

-.00262 

-.00253 

-.00132 

-.00013 

.00127 

.00238 

1 

.00262 

.00252 

.00259 

.00251 

. 00257 

. 00250 

.00255 

.00249 

2 

.00253 

.00248 

. 00252 

.00246 

.00250 

.00245 

. 00249 

.00241 

3 

.00199 

.00143 

.00093 

.00045 

.00043 

.00042 

.00007 

-.00083 

4 

-.00152 

-.00239 

-.00238 

-.00198 

-.OOl*"* 

-.00166 

-.00164 

-.00164 

5 

-.00110 

-.00051 

.00013 

-.00009 

-.00076 

-.00137 

-.00208 

- . 00289 

6 

-.00358 

-.00328 

-.00272 

-.00225 

-.00151 

.00050 

. 00223 

.00268 

7 

.00306 

.00365 

.00445 

. 00522 

.00597 

. 00683 

.00657 

.00482 

8 

.00334 

.00076 

-.00155 

-.00193 

-.00130 

-.00059 

.00027 

.00094 

9 

.00184 

.00156 

.00045 

-.00052 

-.00144 

-.00243 

-.00316 

-.00369 

10 

-.00432 

-.00517 

- . 00583 

-.00571 

-.00534 

-.00507 

-.00384 

-.00231 

11 

-.00113 

-.00012 

.00140 

.00290 

.00450 

.00499 

. 00492 

.00408 

12 

.00285 

.00186 

.00148 

.00103 

.00112 

.00138 

.00182 

.00106 

13 

-.00115 

-.00315 

-.00531 

-.00656 

-.00708 

-.00784 

-.00777 

-.00671 

14 

-.00567 

-.00453 

-.00344 

-.00232 

-.00123 

-.00018 

.00087 

.00210 

15 

.00334 

.00451 

.00436 

.00387 

.00461 

.00555 

.00648 

.00740 

16 

.00835 

. 00925 

.01047 

.01212 

.01227 

.01113 

.01035 

. 00887 

17 

.00679 

.00519 

.00471 

.00381 

.00543 

.01108 

.01199 

.01224 

18 

.01259 

.01378 

.01540 

.01691 

.01850 

.02003 

.02150 

.02154 

19 

.02169 

.01995 

.01535 

.01340 

.01366 

.01322 

.00967 

.00563 

20 

.00439 

.00303 

.00285 

.00363 

.00441 

.00444 

.00461 

.00343 

21 

.00028 

-.00194 

-.00193 

-.00642 

-.00977 

-.01115 

-.01229 

-.01266 

22 

- .01351 

-.01486 

-.01704 

-.01881 

-.02093 

-.02273 

-.02464 

-.02597 

23 

-.02790 

-.02647 

-.01966 

-.01471 

-.00727 

.00229 

.00355 

.00633 

24 

.00455 

.00364 

.00264 

.00281 

.00120 

-.00043 

-.00257 

-.00148 

25 

-.00062 

.00029 

-.00049 

-.00665 

-.01472 

-.01278 

-.03357 

-.04283 

26 

-.04301 

-.04581 

-.04329 

-.04167 

-.03774 

- . 02528 

-.01920 

.00104 

27 

.01365 

.01287 

.01582 

.01638 

.01870 

.01967 

.02134 

.02096 

28 

.02131 

.02285 

.02644 

.02925 

.03198 

.03323 

.03474 

.03176 

29 

.02774 

.02408 

.02034 

.01720 

.01156 

.00759 

-.01569 

- . 02845 

30 

-.02979 

-.03884 

-.04221 

-.04369 

-.04100 

-.03866 

-.03365 

-.03090 

31 

-.02021 

-.00177 

.00202 

.00776 

.00612 

.00471 

. 00752 

.00947 

32 

.01097 

.01191 

.01265 

.02008 

.02870 

.03266 

.03884 

.04095 

33 

.04180 

.03834 

.03296 

.03118 

.02943 

.02614 

, 02507 

.02268 

34 

.02159 

.01664 

.01778 

.01550 

.01922 

-.00501 

- . 01865 

-.03504 

35 

-.04725 

-.04257 

-.04512 

-.03951 

-.03440 

-.02055 

-.00475 

-.00227 

36 

-.00100 

-.00128 

-.00062 

-.00006 

. 00056 

.00012 

-.00103 

-.00092 

37 

- .00065 

.00049 

.00346 

.00663 

.00723 

.00594 

.00573 

.00572 

38 

.00528 

.00313 

.00359 

.00439 

.00172 

-.00095 

. 00023 

.00309 

39 

.00417 

. 00262 

.00344 

.00368 

.00480 

.00515 

. 00535 

.00771 

40 

.01055 

.00919 

.00849 

.00569 

.00501 

.00043 

-.00394 

-.00607 

41 

-.00743 

-.00994 

-.01044 

-.01259 

-.01139 

-.01962 

-.03959 

-.04157 

42 

-.04639 

-.04836 

-.05291 

-.05279 

-.04936 

-.04583 

-.04318 

-.02784 

43 

- .01453 

-.01386 

.02582 

.04311 

.04349 

.04799 

.04768 

.05103 

44 

.05160 

.05464 

.05483 

.05724 

.05502 

.04960 

.04395 

.04052 

45 

.03706 

.03452 

.03170 

.02457 

.01447 

.01168 

.00425 

-.00723 

46 

-.02697 

-.04144 

-.04287 

-.04527 

-.04515 

-.04433 

-.04025 

-.03331 

47 

-.01702 

-.00458 

- .00454 

-.00364 

-.00492 

-.00471 

-.00566 

- . 00584 

48 

-.00773 

-.00807 

-.01077 

-.01240 

-.01443 

-.00916 

-.00783 

.00337 

49 

.02477 

.02571 

.03000 

.03244 

.03488 

.03481 

.03569 

.03595 

SO 

.03125 

.02231 

.01623 

.00857 

-.00306 

-.00778 

-.00940 

-.00945 

51 

-.00879 

- . 00898 

-.00900 

-.00932 

-.00938 

-.00965 

-.00985 

-.01166 

52 

-.01376 

-.01584 

-.01779 

-.01812 

-.01891 

-.01627 

-.01031 

-.00345 

53 

.00263 

.00411 

.00437 

.00123 

-.00119 

-.00505 

-.00883 

-.01227 

54 

- .01374 

-.01533 

-.01761 

- . 02283 

-.02439 

-.02657 

-.02460 

-.02207 

55 

-.01888 

-.01564 

-.01196 

-.01223 

-.01449 

-.01609 

-.01790 

-.02167 

56 

-.02701 

-.02704 

-.01719 

-.00755 

.00148 

.01740 

. 02331 

.02318 

57 

F8 


.02704 

. 02872 

.03258 

.03163 

.03086 

.02832 

.02788 

.02395 

58 

.02337 

.01400 

-.00839 

-.01772 

-.02134 

-.03369 

-.03794 

-.03629 

59 

-.03637 

-.03536 

-.03419 

-.03040 

-.02715 

-.02344 

-.02009 

-.01637 

SO 

-.01240 

-.00832 

-.00431 

-.00076 

.00279 

. 00627 

.01001 

.01360 

51 

.01775 

. 02090 

.01802 

.01284 

.01015 

.00631 

.00431 

-.00277 

62 

-.01157 

-.01311 

-.01272 

-.00909 

-.00429 

.00070 

.00166 

.00317 

63 

.00383 

. 00522 

.00557 

.00531 

.00455 

.00432 

.00227 

-.00201 

64 

-.00549 

-.00814 

-.01094 

-.01397 

-.01732 

-.02096 

-.02351 

-.01983 

65 

-.01612 

-.01444 

-.01289 

-.01131 

-.00930 

- . 00729 

-.00607 

-.00484 

66 

-.00443 

-.00520 

-.00569 

-.00689 

-.00952 

-.01100 

-.00995 

-.00918 

67 

-.00793 

- . 00598 

-.00440 

-.00507 

-.00559 

-.00604 

- . 00284 

.00062 

68 

.00386 

.00637 

.00746 

. 00553 

.00497 

. 00205 

-.00193 

-.00868 

69 

-.01151 

-.01129 

-.01126 

-.01142 

-.01118 

-.01216 

-.01317 

-.01475 

70 

-.01553 

-.01697 

-.01891 

-.02491 

-.02437 

-.02170 

-.01701 

-.01619 

71 

.00794 

.01731 

.01520 

.01773 

.01725 

.01900 

.01898 

.01960 

72 

.01828 

.01823 

.01722 

.01737 

.01670 

.01571 

.01045 

.00818 

73 

.00911 

.01060 

.01142 

.01169 

.01149 

.01166 

.01151 

.01162 

74 

.01182 

.01287 

.01377 

.01384 

.00965 

.00549 

.00165 

-.00121 

75 

-.00517 

-.00923 

-.01388 

-.01813 

-.01579 

-.01011 

-.00601 

- . 00283 

76 

- .00019 

.00286 

.00563 

.00893 

.01189 

.01493 

.01679 

.01897 

77 

.02081 

.02294 

.02566 

.02683 

.02362 

.02188 

.01823 

.01590 

78 

.00642 

-.00320 

- . 00558 

-.01499 

-.01895 

-.02196 

-.02474 

-.02764 

79 

-.03006 

- . 02884 

-.02685 

-.02514 

-.02649 

-.02964 

-.03232 

-.03516 

80 

-.03518 

-.03490 

-.03390 

-.03218 

-.03024 

-.02866 

- . 02522 

-.01998 

81 

-.01493 

-.00919 

-.00326 

.00181 

.00546 

.00931 

. 01268 

.01494 

82 

.01713 

.01925 

.02091 

.02242 

.02402 

. 02581 

. 02957 

.03183 

83 

.03031 

.02831 

.02686 

. 02223 

.01352 

.01016 

.00690 

.00397 

84 

.00060 

-.00049 

.00012 

.00048 

.00099 

. 00126 

.00097 

.00052 

85 

.00024 

-.0?006 

- .00026 

-.00058 

.00054 

.00247 

.00436 

.00595 

86 

.00725 

.00851 

.00988 

.01013 

.00972 

.00911 

.00735 

.00519 

87 

.00493 

.00672 

.00816 

.00980 

.01243 

.01622 

.01876 

.01684 

88 

.01446 

.01174 

. 00858 

.00541 

.00291 

.00168 

.00042 

-.00094 

89 

-.00275 

-.00466 

- . 00668 

-.01230 

-.01539 

-.01571 

-.01660 

-.01718 

90 

-.01796 

-.01836 

-.01420 

-.01042 

- . 00905 

-.00692 

-.00534 

-.00263 

91 

.00000 

.00305 

.00436 

.00421 

.00433 

. 00369 

.00185 

.00008 

92 

-.00181 

-.00304 

-.00347 

-.00407 

-.00450 

-.00383 

-.00163 

-.00023 

93 

.00228 

.00587 

.01472 

.02550 

.02612 

. 02857 

.02996 

.03102 

94 

.02932 

.02846 

.02698 

.  02666 

.02573 

.  02537 

.02410 

.02332 

95 

.02186 

.02105 

.01855 

.01489 

.01201 

.01034 

. 00805 

.00531 

96 

.00176 

-.00260 

-.00736 

-.01071 

-.01549 

-.02166 

-.02159 

-.01965 

97 

-.01836 

-.01657 

-.01522 

-.01349 

-.01121 

-.00704 

-.00651 

-.00786 

98 

-.00923 

-.01047 

-.01283 

-.01568 

-.01658 

-.01079 

-.00523 

. 00058 

99 

.00533 

.00965 

.01214 

.01590 

.01834 

. 02155 

.01990 

.02172 

100 

MAXIMUM  ACCELERATION  -  .05724 

AT  TIME  -  7.10  SEC 


THE  VALUES  WILL  BE  MULTIPLIED  BY  A  FACTOR  -  . 349 

TO  GIVE  NEW  MAXIMUM  ACCELERATION  -  .02000 

MEAN  SQUARE  FREQUENCY  -  1.58  C/SEC. 

MAX  ACCELERATION  -  .02000  FOR  FREQUENCIES  REMOVED  ABOVE  25.00 

C/SEC. 
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3  ***  read  where  object  motion  is  given 


OBJECT  MOTION  IN  LAYER  NUMBER  9  OUTCROPPING 


1******  4  ***  OBTAIN  STRAIN  COMPATIBLE  SOIL  PROPERTIES 


MAXIMUM  NUMBER  OF  ITERATIONS  -  10 
MAXIMUM  ERROR  IN  PERCENT  -  5.00 
FACTOR  FOR  EFFECTIVE  STRAIN  IN  TIME  DOMAIN  -  .65 


EARTHQUAKE  -  PASADENA  1952 
SOIL  PROFILE  -  Example 
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THE  CALCULATION  HAS  BEEN  CARRIED  OUT  IN  THE  TIME  DOMAIN  WITH  EFF.  STRAIN  -  .65*  MAX.  STRAIN 
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SOIL  PROFILE  -  Exarap 


THE  CALCULATION  HAS  BEEN  CARRIED  OUT  IN  THE  TIME  DOMAIN  WITH  EFF.  STRAIN  -  .65*  MAX.  STRAIN 
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THE  CALCULATION  HAS  BEEN  CARRIED  OUT  IN  THE  TIME  DOMAIN  WITH  EFF.  STRAIN  -  .65*  MAX.  STRAIN 
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THE  CALCULATION  HAS  BEEN  CARRIED  OUT  IN  THE  TIME  DOMAIN  WITH  EFF,  STRAIN  -  .65*  MAX.  STRAIN 
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THE  CALCULATION  HAS  BEEN  CARRIED  OUT  IN  THE  TIME  DOMAIN  WITH  EFF.  STRAIN  -  .65*  MAX.  STRAIN 
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MAXIMUM  AMPLIFICATION  -  12.74 
FOR  FREQUENCY  -  1.18  C/SEC. 

PERIOD  -  .84  SEC. 


COMPUTE  MOTION  IN  NEW  SUBLAYERS 


c 

cC 


d  <s 
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1****** 


9  ***  COMPUTE  RESPONSE  SPECTRUM 


COMPUTE  RESPONSE  SPECTRUM  IN  LAYER  1 

RESPONSE  SPECTRUM  ANALYSIS  FOR  LAYER  NUMBER  1 
CALCULATED  FOR  DAMPING  .050 

TIMES  AT  WHICH  MAX.  SPECTRAL  VALUES  OCCUR 
TD  -  TIME  FOR  MAX.  RELATIVE  DISP. 

TV  -  TIME  FOR  MAX.  RELATIVE  VEL. 

TA  -  TIME  FOR  MAX.  ABSOLUTE  ACC. 
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DAMPING  RATIO 


.05 


PERIOD 


TD 

.00 

7.6800 

.10 

7.6600 

.15 

7.6R00 

.20 

7.72/1 

.25 

7.3000 

.30 

7.3200 

.35 

8.5200 

.40 

7.6600 

.45 

13 . 0000 

.50 

7.7000 

.55 

7 . 7400 

.60 

7.7800 

.65 

8.5600 

.70 

6.4200 

.75 

6.4800 

.80 

8.6800 

.85 

8.7600 

.90 

8.8200 

.95 

8.8800 

'  00 

8  4800 

1.10 

8 . 1000 

1.20 

8 . 1000 

1.30 

8.1000 

1.40 

7.6600 

1.50 

7.6800 

1.60 

15.5600 

1.70 

15.6000 

1.80 

15.6600 

1.90 

7.6600 

2.00 

7.6600 

2.25 

13.4200 

2.50 

7 . 2400 

2.75 

7.2800 

3.00 

7 . 3000 

3.25 

5.9600 

3.50 

8.1400 

3.75 

6.7800 

4.00 

15.6600 

TIMES  FOR  MAXIMA  (SEC) 


TV 

TA 

UESHAKE 

SHAKE* 

7.1800 

7.68 

7.66 

8 . 3400 

7.66 

7.66 

7 . 4400 

7.66 

7.64 

7 . 8000 

7.70 

7.70  1 

6.0200 

7.30 

7.28 

8 . 3400 

7.32 

7.30 

8.4000 

8.52 

8.50 

7.5200 

7.66 

7.62 

12.8800 

13.00 

12.94 

7 . 5400 

7.68 

7.66 

7 . 9000 

7.74 

7.72 

7.9400 

7.76 

8.08 

8.0200 

8.56 

8.54 

6 . 2400 

6.40 

6.38 

6.3000 

6.46 

6.44 

8.8800 

8.66 

8.62 

8.9600 

8.74 

8.70 

8.6200 

8.82 

8.78 

8.6600 

8.86 

8.40 

8.7000 

8.48 

8.44 

8.3200 

8.08 

8.04 

8.3200 

8.08 

8.06 

8.3200 

8.08 

8.06 

7.8800 

7.64 

7.60 

7.4600 

7.66 

7.64 

7.4800 

15.54 

15.50 

7.4800 

15.58 

15.54 

7.8800 

15.64 

15.62 

7.8800 

7.62 

7.60 

7.9000 

7.64 

7.62 

7.4800 

13.40 

13.38 

7.4800 

7.20 

7.18 

7.5000 

7.24 

7.22 

8 . 3400 

7.24 

7.22 

8.3400 

5.90 

5.88 

7.8800 

8.08 

6.68 

7.8800 

6.74 

6.70 

7.8800 

15.60 

7.60 

*  Results  using  SHAKE  do  not  include  correction  to  constitutitive  model  by 
Udaka  and  Lysmer  (1973). 
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SPECTR/'L  VALUES -- 

PASADENA  1952  Example  DAMPING  RATIO 
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The  results  using  SHAKE  do  not  reflect  update  to  constitutive  model  by  Udaka  and  Lysmer  (1973). 


VALUES  IN  PERIOD  RANGE  .1  TO  2.5  SEC. 


AREA  OF  ACC.  RESPONSE  SPECTRUM  -  .344 
AREA  OF  VEL.  RESPONSE  SPECTRUM  -  1.877 
MAX.  ACCELERATION  RESPONSE  VALUE  -  .553 
MAX.  VELOCITY  RESPONSE  VALUE  -  2.366 
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CURVE  1 


5.00  %  DAMPING 


1 


ABSISSA 

CURVE  1 

.001 

.102 

.100 

.107 

.150 

.109 

.200 

.116 

.250 

.121 

.300 

.136 

.350 

.153 

.400 

.171 

.450 

.197 

.500 

.227 

.550 

.267 

.600 

.320 

.650 

.342 

.700 

.345 

.750 

.335 

.800 

.504 

.850 

.553 

.900 

.473 

.950 

.367 

1.000 

.283 

1.100 

.158 

1.200 

.097 

1.300 

.086 

1.400 

.074 

1.500 

.063 

1.600 

.050 

1.700 

.048 

.800 

.040 

1 . 900 

.034 

2 . 000 

.038 

.  250 

.035 

.  500 

.024 

.1.750 

.021 

3 . 000 

.015 

3.250 

.010 

3.500 

.008 

3.750 

.007 

♦ .  000 

.007 

F24 


PLOT  OF  VELOCITY  SPECTRA 
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PERIOD  IN  SECONDS 


CURVE  1  5.00  %  DAMPING 


AfiSISSA 

CURVE  1 

.001 

.000 

.100 

.008 

.150 

.022 

.200 

.036 

.250 

.066 

.300 

.102 

.350 

.160 

.400 

.248 

.450 

.368 

.500 

.429 

.550 

.623 

.600 

.898 

.650 

.984 

.700 

1.169 

.750 

1.232 

.800 

1.999 

.850 

2.366 

.900 

2.235 

.950 

1.939 

1.000 

1.587 

1.100 

1.098 

1.200 

.830 

1.300 

.804 

1.400 

.770 

1.500 

.697 

1.600 

.632 

1.700 

.531 

1.800 

.553 

1.900 

.607 

2.000 

.618 

2.250 

.681 

2.500 

.666 

2.750 

.560 

3.000 

.502 

3.250 

.462 

3.500 

.413 

3.750 

.426 

4.000 

.420 

Total  execution  time  -  120.91  secs 
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APPENDIX  G: 
FLOWCHART  OF  WESHAKE 


APPENDIX  H: 

EXAMPLE  SPECIFICATION  FILE 


4096  EXAMPLE  PROBLEM  FOR  WESHAKE  USER'S  MANUAL 
2  OPTION  2;  READ  SOIL  COLUMN  DATA 


0 

9 

4 

.06240 

EIGHT 

LAYERS 

OVERLYING 

ROCK 

1 

3 

1 

8.0 

.45 

.050 

.120 

43. 

43.1.0001 

2 

4 

1 

12.0 

.00 

.050 

.100 

430. 

430.1.0000 

3 

5 

1 

10.0 

.00 

.050 

.100 

745. 

745.1.0000 

A 

5 

1 

11.0 

.00 

.050 

.100 

920. 

920.1.0000 

5 

2 

1 

20.0 

.45 

.050 

.125 

76. 

76.1.0001 

6 

6 

1 

I'ii.O 

.00 

.050 

.125 

910. 

910.1.0000 

7 

7 

1 

20.0 

.00 

.050 

.125 

1090. 

1090.1.0000 

8 

8 

1 

20.0 

.00 

.050 

.125 

1155. 

1155.1.0000 

9 

1 

1 

.020 

.150 

8000. 

8000.1.0000 

8  OPTION  8:  READ  MATERIAL  PROPERTIES 


8 

0  2 

100. 

STRAIN 

COMPATIBLE  PROPERTIES 

8100.0 

ROCK  (Schnabal  1973) 

.0001 

.0003 

.0010 

.0030 

.0100 

.0300 

.1000 

1.0000 

1.0000 

1.0000 

.9900 

.9500 

.9000 

.8100 

.7250 

.5500 

5  5.0 

ROCK  (Schiuibcl  1973) 

.0001 

.0010 

.0100 

.1000 

1.0000 

.4000 

.8000 

1.5000 

3.0000 

4.6000 

9100.0 

SANDj  l.ow«r 

Bound  <S«od  &  IdrLss 

1970) 

.0001 

.0003 

.0010 

.0030 

.0100 

.0300 

.1000 

.3000 

1.0000 

1.0000 

.9850 

.9300 

.8300 

.6350 

.4250 

.2250 

.1100 

.0400 

9  5.0 

SARD,  Lomr 

Bound  (Scad  4  Idrlaa 

1970) 

.0001 

.0003 

.0010 

.0030 

.0100 

.0300 

.1000 

.3000 

1.0000 

.3000 

.4000 

.7000 

1.4000 

2.7000 

5.0000 

9.8000 

15.0000 

20.7000 

9100.0 

SAND,  Av«rat«  (S««d  4 

Idclta  1970) 

.0001 

.0003 

.0010 

.0030 

.0100 

.0300 

.1000 

.3000 

1.0000 

1.0000 

.9800 

.9500 

.8900 

.7300 

.5200 

.2900 

.1400 

.0600 

9  5.0 

SAND>  Av«r«g«  (Seed  & 

Idrlss  1970) 

.0001 

.0003 

.0010 

.0030 

.0100 

.0300 

.1000 

.3000 

1.0000 

.8000 

1 . 0000 

1.9000 

3.0000 

5.4000 

9.6000 

15.4000 

20.8000 

24.6000 

9100.0 

CLAY/SILT  (PI-5-10,  Sun  et  al .  1988) 

.0001 

.0003 

.0010 

.0030 

.0100 

.0300 

.1000 

.3000 

1.0000 

1.0000 

1.0000 

.9750 

.9100 

.7800 

.5650 

.3050 

.1400 

.0400 

9  5.0 

CLAY,  Lowar 

Bound  (Sood  4  Idrias 

1970) 

.0001 

.0003 

.0010 

.0030 

.0100 

.0300 

.1000 

.3000 

1.0000 

1.3000 

1.3000 

l.30r 

1.5000 

1.7000 

2.5000 

4.0000 

6.5000 

12.3000 

H3 


9100.0 

CLAY /SILT 

(PI-10-20, 

Sun  «t  al . 

1988) 

.0001 

.0003 

.0010 

.0030 

.0100 

.0300 

.1000 

.3000 

1.0000 

1.0000 

1.0000 

1.0000 

.9600 

.8700 

.7000 

.4100 

.2000 

.0800 

9  5.0 

Cl^Y,  Average  (Seed 

&  Idrlss  1970) 

.0001 

.0003 

.0010 

.0030 

.0100 

.0300 

.1000 

.3000 

1.0000 

2.5000 

2.5000 

2.5000 

3.2000 

4.5000 

6 . 5000 

9.0000 

13.5000 

20.5000 

9100.0 

CLAY /SILT 

(PI-20-40, 

Sun  et  al . 

1988) 

.0001 

.0003 

.0010 

.0030 

.0100 

.0300 

.1000 

.3000 

1.0000 

1.0000 

1.0000 

1.0000 

.9700 

.9000 

.7700 

.5200 

.3000 

.1400 

9  5.0 

CLAY,  Average  (Seed 

4  Idrias  1970) 

.0001 

.0003 

.0010 

.0030 

.0100 

.0300 

.1000 

.3000 

1.0000 

2.5000 

2.5000 

2.5000 

3.2000 

4.5000 

6.5000 

9.0000 

13.5000 

20.5000 

9100.0 

CLAY/ SILT 

(PI-40-80, 

Sun  et  al . 

1988) 

.0001 

.0003 

.0010 

.0030 

.0100 

.0300 

.1000 

.3000 

1.0000 

1.0000 

1.0000 

1 . 0000 

.9850 

.9200 

.8150 

.6200 

.4100 

.2000 

9  5.0 

CLAY,  Average  (Seed 

4  Idrlas  1970) 

.0001 

.0003 

.0010 

.0030 

.0100 

.0300 

.1000 

.3000 

1.0000 

2.5000 

2.5000 

2.5000 

3.2000 

4.5000 

6.5000 

9.0000 

13.5000 

20.5000 

9100.0 

CLAY/SILT 

(PI>80,  Sun  et  al.  1988) 

.0001 

.0003 

.0010 

.0030 

.0100 

.0300 

.1000 

.3000 

1.0000 

1.0000 

1.0000 

1.0000 

.9850 

.9400 

.8600 

.7100 

.5300 

.3300 

9  5.0 

CLAY  Upper 

Bound  (Seed  Idrlss 

1970) 

.0001 

.0003 

.0010 

.0030 

.0100 

.0300 

.1000 

.3000 

1.0000 

4.0000 

4.0000 

4.0000 

5.0000 

7.5000 

11.0000 

16.0000 

21.8000 

27.0000 

1 

OPTION  1:  READ  ACCELERATION 

VALUES 

FROM 

"EARTHQ 

2048  4096  0 

. 02  Alaska  EQ 

(7/30/72)  Sitka  Record 

] 

..0 

0. 

25.0 

0 

0.00000 

-0.00434 

0.00860 

0.00540 

-0.00565 

-0.00944 

-0.00369 

-0.00669 

1 

-0.00336 

-0.00111 

0.00358 

0.00303 

-0.00323 

-0.00907 

-0.01522 

-0.01029 

2 

-0.00706 

-0.00194 

0.00135 

0.00191 

0.00743 

0.00043 

-0.00657 

0.00527 

3 

-0.01948 

-0.00854 

-0.01807 

-0.01060 

-0.00396 

0.00315 

0.01088 

0.01252 

4 

0.00859 

0.00066 

-0.00698 

-0.01661 

-0.01454 

-0.00959 

-0.00047 

0.01193 

5 

u. 02522 

0.02534 

0.02107 

-0.00808 

-0.01286 

-0.00702 

-0.01510 

-0.00409 

6 

0.00969 

0.00253 

-0.00687 

-0.01304 

-0.00716 

0.00005 

0.00802 

0.00310 

7 

-0.00076 

0.00385 

-0.00101 

0.00063 

-0.00544 

-0.00056 

0.00393 

0.00958 

8 

0.01092 

0.00872 

0.00348 

-0.00373 

-0.00307 

-0.00171 

0.00996 

0.01039 

9 

0  00378 

0.00255 

0.00594 

0.00241 

0.01074 

0.00452 

-0.00165 

-0.01075 

10 

-0.01783 

-0.00081 

0.01251 

0.01050 

0.00821 

0.01259 

0.01636 

0.01832 

11 

O.OOS21 
-0.02*20 
0.00012 
0.05*49 
■0.00108 
-0.02062 
-0.00322 
0.00005 
0.03332 
■0.03177 
'0.0*843 
-0.00*22 
■0.01152 
0.027*5 
0.02713 
-0.01686 
1>.C0*39 
0.00957 
0.01372 
0.0009* 
0.002*1 
0.00*99 
0.00590 
0.00979 
0.01818 
0.0137* 
0.01*62 
0.03302 
0  1*. 
0.^.301 
0.00811 
0.01*90 
0.03136 
■0.01759 
0.02537 
0.000*4 
0.0006* 
0.01889 
0.00581 
0.0002* 
0.00271 
0.0039* 
■0.01959 
0.00606 
0.01610 
0.00*89 
0.01377 
0.01*4* 
0.02*39 
0.00761 
0.02696 
0.0323* 
0.03831 
0.05187 


0.00569 

0.02808 

-0.00870 

0.01070 

-0.03383 

-0.03*49 

-0.020*0 

-0.01511 

-0.00*55 

-0.01033 

0.00853 

-0.0061* 

0.01078 

0.00331 

0.0103* 

-0.01236 

0.00*85 

0.01188 

-0.0221* 

-0.00631 

0.01099 

-0.02736 

0.00388 

-0.0082* 

-0.01190 

0.01379 

0.00183 

0.02109 

-0.03076 

-0.02333 

0.015*7 

0.00*79 

0.02286 

-0.0090* 

-0.01385 

-0.00860 

0.012*5 

-0.02958 

-0.01801 

-0.00585 

0.0087* 

0.01529 

-0.01197 

0.00099 

0.02566 

0.00926 

-0.01568 

0.00361 

-0.01251 

'0.01182 

0.022*0 

-0.00119 

-0.02617 

-0.05828 


0.00351 
0.0*292 
-0.01*17 
-0.02C89 
-0.0*305 
-0.05223 
-0.036*7 
-0.01056 
-0.01688 
0.02196 
0.02995 
0.00187 
0.00*64 
-0.01182 
-0.00558 
0.01295 
-0.00509 
0.01098 
-0.0272* 
0.00213 
0.01575 
-0.00*09 
0.01022 
-0.0168* 
0.00168 
0.00*31 
-0.00603 
0.00173 
-0  "  .47 
-0.01350 
0.00337 
0.01760 
0.01396 
-0.00*02 
-0.0051* 
-0.02157 
-0.001*3 
-0.0277* 
-0.01608 
-0.00319 
0. J1993 
0.005*7 
0.00786 
-0.00909 
0.01704 
0.001*0 
-0.02369 
0.01523 
0.00092 
-0.02*5* 
0 . 00582 
-0.0*61* 
0.01885 
-0.04*87 


-0.01063 
0.05900 
-0.02311 
-0.0*311 
-0.00785 
-0.06713 
-0.00725 
-0.02773 
-0.00372 
0.02633 
0.0*202 
0.00510 
-0.00813 
-0.00357 
0.00056 
0.0*609 
-0.03677 
0.00701 
-0.01255 
0.02019 
0.00035 
0.01117 
0.0CC29 
0.00068 
0.01*22 
-0.00  >11 
-0.01695 
-0.01092 
-0.00982 
0.00751 
-0.0063* 
0.0250* 
0.01*45 
-0.00582 
0.00377 
-0.02092 
'0.00622 
-0.01258 
0.00*52 
0.0015* 
0 . 02882 
-0.00612 
0.017*9 
-0.02959 
0.01297 
-0.00020 
-0.01927 
0.00286 
0.00103 
-0.01078 
0.01180 
-0.0*519 
0.080*4 
-0.01279 


-0.03369 
0 . 05520 
0.00280 
-0.0*94* 
0.01159 
-0.03791 
0 . 02286 
-0.0*252 
0.02375 
-0.00*78 
0.05*4* 
0.00357 
-0.00155 
0.0011* 
-0.01203 
0.01396 
-0.0*885 
-0 . 00082 
0.00013 
0.0*131 
-0.01550 
-0.01109 
0.00519 
0.01020 
0.00*52 
-0.00300 
-0.012*4 
0.00355 
0.00613 
0.01*15 
-0.00365 
0.00820 
0.00588 
-0.00260 
-0.00156 
-0.00729 
-0.00229 
0.02002 
0.01662 
0.00773 
0.02*28 
0.01018 
-0.0061* 
-0.0*610 
0.01980 
0.01116 
0.00197 
-0.02275 
0.01293 
0.02007 
0.00672 
-0.03707 
0.09012 
0.00973 


-0.0*15* 

0.0*862 

0.03229 

-0.01*7* 

0.02328 

0.01826 

0.02*32 

-0.02969 

0.04392 

-0.03830 

0.0*226 

-0.0076* 

0.01288 

0.00*61 

-0.02606 

0.00772 

-0.021*8 

-0.01205 

-0.00*18 

0.01*11 

-0.00539 

-0.0161* 

0.01*08 

0.00*42 

-0.00203 

0.00*76 

-0.00316 

0.01733 

-0.00193 

0.00*83 

0.011*7 

-0.0109* 

-0.00117 

0.00196 

-0.00361 

0.01*71 

-0.01037 

0.029*0 

0.02330 

-0.00708 

0.009*0 

0.01381 

-0.03373 

-0.0*088 

0.0101* 

0.01507 

0.027*4 

-0.0*107 

0.03030 

0.01272 

0.016*9 

-0.02089 

0.06607 

-0.01*3* 


-0.0*725 
0.03961 
0 . 0*878 
0.01639 
0.0*201 
0.01*50 
0.02051 
0.00359 
0.011*4 
-0.05378 
-0.011*6 
-0.01967 
0.00717 
0 . 02062 
-0.020*3 
0.02510 
0.0015* 
-0.02089 
-0.00666 
-C. 00789 
0.0198* 
-0.00198 
0.00196 
0.00309 
-0.00786 
0.0087* 
0.00906 
0.01532 
-0.00031 
'0.00811 
0.00627 
0.00392 
-0.00988 
-0.00326 
0.00225 
-0.00015 
-0.01635 
0.038*3 
0.01926 
-0.00*9* 
-0.00351 
-0.00033 
'0.0*386 
-0.02783 
0.00*40 
0.00220 
0.0*390 
-0.03273 
0.01*05 
0.00832 
0.02120 
0.00107 
0.0275* 
-0.0*368 


-0.05131 
0.02367 
0.05*55 
0.03284 
0.01*61 
0.00062 
0.01305 
0.03776 
-P.02161 
-0.06298 
-0.01*36 
-0.01591 
0.02209 
0.03680 
-0.00679 
0.01280 
0.00619 
-0.009*2 
0.00568 
-0.00313 
0.01635 
0.0061* 
-0.00296 
-0.01187 
0.002*1 
-0.00625 
0.015*2 
0.0107* 
0.00152 
0.00230 
-0.01083 
0 . 02532 
-0.01560 
-0.012*1 
0.00130 
-0.01126 
-0.01162 
0.02579 
0.00635 
0.00160 
-0.01617 
-0.01355 
-0.012*7 
-0.0009* 
0.0066* 
-0.01522 
0.02878 
-0.01835 
-0.00277 
0.02089 
0.030*8 
-0.00696 
-0.03119 
-0.02830 


12 

13 

1* 

15 

16 

17 

18 

19 

20 
21 
22 
23 
2* 

25 

26 

27 

28 

29 

30 

31 

32 

33 
3* 

35 

36 

37 

38 

39 
*0 
*1 
*2 
*3 
** 
*5 
*6 
*7 
*8 
*9 

50 

51 

52 

53 
5* 

55 

56 

57 

58 

59 

60 
61 
62 
63 
6* 
65 


H5 


0.01788  0.033S8  0.0092*  -0.00*73  -0.01*37  -0.02338  -0.00393  0.00820  66 

-0.00180  -0.02392  -0.01361  -0.0021*  0.01*76  0.01*71  0.01820  0.009*7  67 

0.0161*  0.02600  0.01371  0.00223  -0.01201  -0.02*89  -0.01608  0.00081  68 

0.00730  0.00*27  -0.01*67  -0.0090*  0.00908  0.01007  -0.00970  -0.006*7  69 

0.00731  -0.01211  -0.02071  -0.01328  -u. 02376  -0.029*7  0.001*3  0.01660  70 

0.01238  0.02222  0.01102  -0.02331  -0.02830  -0.01223  0.01122  0.02306  71 

0.00341  -0.01089  -0.002*3  0.013*8  0.00811  -0.00038  0.0087*  -0.00061  72 

-0.00681  -0....43'’  0.010*3  0.00382  -0.00667  -0.01692  -0.0038*  -0.00*5*  73 

-0.01270  -0.00729  0.00693  0.01899  0.00332  -0.01323  -0.01803  -0.0033*  7* 

0.00988  0.00626  0.00619  0.01939  0.03712  0.01838  -0.00837  -0.011*9  73 

0.00396  0.01*72  0.00917  0.01693  0.00907  0.01768  0.02319  0.02032  76 

0.01729  0.01970  0.010*6  -0.00827  -0.00660  0.00362  0.0083*  -0.00078  77 

0.00927  0.01710  0.01262  0.01088  0.02193  0.03039  0.00867  -0.01616  78 

-0.01259  -0.00820  -0.01671  -0.00971  -0.00055  -0.00930  -0.013*8  0.01078  79 

0.01920  0.001**  -0.00133  0.00*79  -0.01090  -0.02128  0.00336  0.04*38  80 

0.03903  0.028*8  0.02266  0.00809  -0.001*8  0.00338  0.01329  0.01168  81 

0.00763  0.01111  -0.00232  -0.0329*  -0.0*362  -0.0*408  -0.05987  -0.07305  82 

-0.06930  -0.01932  0.06129  0.069*3  0.03233  0.0333*  0.00823  -0.03993  83 

-0.03839  -0.01*79  0.00379  0.0011*  0.00330  0.008*5  -0.00170  -0.01027  8* 

-0.01305  -0.02255  -0.02061  -0.02808  -0.02691  -0.00266  0.01325  0.00*93  83 

-0.00383  -0.01081  0.00103  -0.00*30  -0.00062  -0.00193  -0.01392  -0.02768  86 

-0.03897  -0.0393*  -0.01737  0.00296  0.00386  -0.0080*  0.00*62  0.00713  87 

-0.01*93  -0.02255  -0.00*67  0.00127  -0.00793  -0.00719  -0.01665  0.00*87  88 

0.02578  0.03937  0.01163  -0.02163  -0.03282  -0.021*0  0.00*10  0.029*7  89 

0.00922  -0.00739  0.01118  0.02992  0.023*6  0.00331  0.00671  -0.00118  90 

0.00363  0.00355  -0.00893  -0.01157  -0.01369  -0.01951  -0.00*73  0.01539  91 

0.01076  0.01138  0.01672  0.00396  -0.00257  0.000*8  0.00939  0.00306  92 

-0.00770  -0.00222  -0.01167  -0.01913  -0.0063*  -0.00170  0.00*55  0.00712  93 

-0.00183  0.01029  0.00928  -0.01023  -0.03021  -0.0281*  -0.02520  -0.00*82  9* 

0.01621  -0.00706  -0.021*0  -0.01215  0.018..8  0.03323  0.0*160  0.0*932  93 

0.0338*  0.01602  0.01*20  0.021*8  0.00067  -0.01882  -0.03701  -0.02130  96 

0.00159  -0.01039  -0.0203*  -0.02035  -0.02096  -0.03021  -0.01721  0.00120  97 

0.00999  0.03996  0.03998  0.03830  -0.00*92  -0.02001  0.00320  0.02389  98 

0.00*97  -0.00600  -0.0176*  -0.02316  -0.02193  -0.00326  0.01606  0.0088*  99 

0.00004  0.00*68  0.01070  -0.00262  -0.00310  0.00693  0.00391  -0.01337  100 

-0.0186*  -0.00*51  0.01710  0.032*3  0.01531  -0.00271  0.01398  0.03107  101 

0.02225  0.00668  -0.00579  -0.01110  0.000*1  0.00267  -0.01021  -0.0187*  102 

-0.00005  0.02106  0.0109*  -0.00310  -0.00377  -0.00199  0.00007  -0.0026*  103 

0.00*10  0.01021  -0.0010*  -0.00563  0.00179  0.009*0  -0.000*3  -0.00336  10* 

0.00295  0.01777  0.02610  0.01180  -0.00021  -0.01087  -0.02012  -0.02598  105 

-0.00*08  0.01120  -0.00623  -0.00*60  -0.00328  -0.0061*  -0.01131  -0.01637  106 

-0.00265  0.01229  0.00288  -0.01200  -0.01732  -0.00575  0.00976  0.01836  107 

-0.00139  0.00336  -0.0017*  0.000*0  -0.00911  0.00089  0.01*99  0.02777  108 

0.01212  -0.00123  0.00050  -0.00230  -0.00115  0.00*99  0.01707  0.00707  109 

0.00979  0.01802  0.00866  -0.00113  -0.0075*  -0.01593  -0.0213*  -0.01153  110 

-0.00258  0.00*66  0.02393  0.03739  0.01*61  -0.00099  -0.00361  -0.00*31  111 

-C.0032>  <1  00660  -0.00998  -0.00790  0.00162  0.00033  -0.00901  -0.0133*  112 

-0.02*81  -0.01925  -0.00260  -0.00136  -0.00778  0.00029  -0.00380  -0.01267  113 

-0.020*0  -0.00911  -0.01991  -0.01799  -0.0106*  0.00973  0.00196  0.0020*  11* 

-0.00660  -0.01336  -0.001*1  0.01197  0.02***  0.02398  0.00893  0.00091  113 

-0.00723  -0.010*8  -0.01393  -0.0192*  -0.00975  -0.00233  -0.00383  -0.00612  116 

-0.00923  -0.01318  -0.01172  -0.00663  0.00208  0.00322  0.010*3  0.0113*  117 

0.307*6  0.00*85  0.01*37  0.01779  0.010*4  0.01*3*  0.01263  0.00298  118 

0.00381  0.012*3  0.00021  -0.01136  -0.02783  -0.03070  -0.02931  -0.02010  119 
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-0.00*16  0.008*0  0.00670  0.01730  0.03012  0.02617  0.02023  0.01081  120 

-0.00*17  -0.00796  -0.00138  0.006*9  0.013*8  -0.00262  -0.02233  -0.0180*  121 

-0.00719  0.01786  -0.00060  -0.00602  0.00219  -0.00869  -0.01331  -0.01373  122 

-0.00*6*  0.00706  0.0117*  0.0003*  -0.00718  0.00630  0.00663  0.00377  123 

0.00272  0.00699  0.00*13  0.00*91  0.003*7  -0.00313  -0.00107  0.01029  12* 

0.00917  0.00173  -0.00861  0.00103  0.009*6  -0.00*33  -0.01179  -0.010**  123 

0.00108  0.01373  0.00313  -0.00336  -0.00869  -0.00823  -0.00390  -0.01303  126 

-0.017*9  -0.01077  -0.00*22  -0.00*69  -0.01339  -0.01739  -0.00293  0.009*9  127 

0.00762  0.00523  -0.00012  -0.00810  -0.00398  0.00033  0.01293  0.01819  128 

0.00800  -0.00*02  -0.00993  -0.00809  0.00630  -0.0003*  -0.00866  -0.00029  129 

0.01*37  0.00786  0.01239  0.02087  0.01806  0.01357  0.01173  0.008*1  130 

0.00167  0.00807  0.01183  0.00321  -0.00817  -0.00706  0.00138  -0.00155  131 

0.00729  0.00030  -0.00792  -0.00127  0.01071  0.00333  0.00092  0.00*18  132 

0.00683  0.01*07  0.013*7  0.01118  0.00739  0.00137  -0.003*6  -0.00068  133 

0.00*98  0.00313  -0.00399  -0.0123*  -0.00801  0.00013  0.00773  0.0011*  13* 

0.00293  -0.00637  -0.01397  -0.01907  -0.02*69  -0. 01**3  0.00*39  0.01*19  133 

0.01323  0.00331  -0.0109*  -0.01070  -0.01037  -0.01861  -0.01630  -0.003*6  136 

0.00370  0.01113  0.01*72  0.000*7  -0.0038*  0.00039  0.00703  0.00023  137 

0.001*8  0.006*9  -0.00782  -0.00223  0.01217  0.00813  -0.00778  -0.00269  138 

0.01192  0.00161  -0.00220  0.00712  0.00298  -0.0031*  0.00369  0.00963  139 

0.00*63  0.00098  0.01031  0.0067*  0.00290  0.00612  0.01218  0.0089*  1*0 

0.00269  -0.00*78  -0.01337  -0.0099*  -0.00*09  0.00*30  -0.00371  -0.0173*  1*1 

-0.02829  -0.02100  -0.017*3  -0.02372  -0.00931  -0.00389  -0.00122  0.00683  1*2 

0.0112*  0.00*61  0.00003  0.01231  0.01193  0.00*99  -0.00288  -0.003*6  1*3 

-0.00*17  -0.00*33  -0.00923  -0.00*88  -0.00203  -0.00793  -0.01602  -0.01367  1** 

-0.00783  0.00073  0.002*7  0.00332  0.00881  0.01613  -0.00013  -0.007*6  1*3 

-0.00116  -0.00038  -0.00739  0.00*08  0.01*93  0.00833  0.00733  0.00777  1*6 

-0.00071  -0.007*3  -0.00817  -0.00267  0.00723  0.00277  -0.00133  0.0036*  1*7 

0.001*0  -0.00366  -0.00300  0.00377  0.00697  0.00*17  0.00776  0.00667  1*8 

0.003*7  -0.00393  -0.0060*  -0.00667  -0.00916  -0. 01**1  -0.01089  -0.00789  1*9 

-0.00331  -0.00096  0.00333  0.01282  0.02221  0.01700  0.00333  -0.006*3  150 

-0.00810  -0.0032*  -0.00812  -0.01182  -0.01137  -0.00318  -0.01382  -0.00378  131 

0.0030*  -0.00362  -0.0133*  -0.00662  0.00006  -0.00010  0.00371  0.01393  132 

0.01*62  0.00720  0.00318  -0.0000*  -0.00*82  -0.00*68  0.0020*  0.00198  133 

-0.002*8  0.00020  0.00301  0.00366  0.00292  -0.00116  -0.0007*  0.00369  13* 

0.00309  0.0008*  0.00*33  0.00283  -0.001*9  0.00130  0.00629  0.0063*  133 

0.0013*  0.00009  0.00788  -0. 0^9*3  -0.00186  0.0067*  0.01803  0.0109*  136 

0.00796  0.00370  -0.00333  -0.01033  -0.01017  -0.01023  -0.00701  -0.00323  137 

0.00117  0.00667  -0.00035  -0.01017  -0.01683  -0.01599  -0.00368  0.00526  138 

0.01373  0.01238  0.00*99  0.00882  0.01371  0.00870  0.00*69  -0.00032  139 

0.000*6  0.00361  0.00819  0.00272  -0.00303  -0.00032  0.00037  -0.00*36  160 

-0.00197  -0.00303  -0.00329  -0.00361  -0.00960  -0.00337  -0.00028  -0.00372  161 

-0.00303  -0.00*33  -0.00093  0.00110  0.00263  0.01113  0.01007  0.0030*  162 

-0.00*91  -0.008*3  -0.00120  0.00839  0.01086  0.01663  0.02560  0.02103  163 

0.01*06  0.01299  0.00660  -0.00022  0.00966  0.00089  -0.00783  -0.0167*  16* 

-0.00880  -0.00620  -0.0093*  -0.00900  -0.00339  -0.00172  -0.0015*  -0.00139  163 

-0.00233  -0.00293  -0.00298  -0.00379  -0.00*88  -0.00*83  -0.00619  -0.0028*  166 

0.0026*  0.00027  0.000*6  -0.00399  -0.00866  -0.00*81  -0.00110  -0.00*75  167 

-0.01372  -0.00871  0.0006*  -0.001*3  -0.00382  0.00011  0.00778  0.000*9  168 

-0.008*2  -0.00086  0.00937  0.00399  -0.00191  -0.00*65  -0.0002*  0.00593  169 

0.01396  0.00737  -0.00318  -0.00572  -0.00376  -0.00872  -0.01399  -0.00976  170 

0.00127  -0.0006*  0.00*12  0.01191  0.00*72  -0.00737  -0.00015  0.00802  171 

-0.00227  -0.003*0  -0.000*0  -0.0023*  -0.0023*  0.00959  0.0193*  0.01*81  172 

0.01338  0.0117*  0.0071*  0.00*31  0.00713  0.000*1  -0.00738  -0.00123  173 
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0.00326  0.00031  -0.01406  -0.01610  -0.00978  -0.00691  -0.01393  -0.01908  174 

-0.01067  -0.00362  0.00773  0.00369  0.00843  0.01329  0.00978  0.00003  175 

-0.00444  0.00691  0.01310  0.00232  0.01163  0.00343  -0.00324  0.00427  176 

0.01460  0.00533  -0.00439  0.00164  -0.00869  -0.01688  -0.01292  -0.00923  177 

-0.01491  -0.00730  0.01233  0.00462  -0.00211  0.00310  0.00652  0.00391  178 

-0.00674  -0.00387  0.00429  0.00288  -0.00264  0.00252  0.00834  0.00018  179 

-0.00639  0.00110  -0.00003  -0.00632  0.00328  -0.00140  -0.01767  -0.01005  180 

0.00233  0.00410  0.00229  0.01794  0.01037  -0.01380  -0.02630  -0.00572  181 

0.01804  0.01007  -0.00375  -0.00191  0.00433  -0.00434  -0.01267  -0.02232  182 

-0.00863  0.00292  0.00633  0.01642  0.00621  -0.00904  -0.01407  0.00341  183 

0.01094  0.00423  -0.01288  -0.01305  -0.01069  -0.00231  -0.00348  -0.00155  184 

0.00349  0.00367  0.00618  0.00650  0.00551  0.00083  -0.00506  0.00128  185 

0.01664  0.00907  -0.00772  -0.01304  -0.00309  0.00273  -0.00175  -0.00238  186 

0.00014  0.00405  0.00676  0.00696  -0.00069  -0.00173  0.00307  0.00024  187 

0.00461  -0.00160  -0.00588  -0.00888  -0.00027  0.00441  0.00029  -0.00213  188 

-0.00195  0.00085  0.00575  0.01111  0.00865  0.00504  0.00339  0.00263  189 

0.00865  -0.00472  -0.01937  -0.01073  0.00773  0.01535  0.00532  -0.00272  190 

-0.00626  0.00019  0.00218  0.00178  0.00861  0.01788  0.02463  0.01591  191 

0.00537  0.00333  0.01023  0.01439  0.00877  0.00134  0.00161  -0.00218  192 

-0.00871  -0.01528  -0.00942  0.00272  0.00657  -0.00581  -0.02075  -0.01980  193 

-0.01109  -0.00358  0.00258  -0.00117  -0.00444  0.00247  -0.00274  -0.00953  194 

-0.00832  -0.00303  0.00896  0.01694  0.01604  0.00847  -0.00449  -0.01003  195 

-0.01478  -0.00969  -0.00734  -0.00392  0.00215  0.00973  0.00837  0.00574  196 

0.01087  0.00418  -0.00381  -0.00441  -0.01057  -0.01705  -0.01061  -0.00392  197 

-0.00857  -0.01133  -0.00758  -0.00265  0.00192  -0.00093  -0.00734  -0.00193  198 

0.00078  0.00111  0.00130  0.00223  0.00368  0.01050  0.01233  0.00451  199 

-0.00419  -0.00129  0.00963  0.00933  0.00408  0.00636  0.00824  0.00264  200 

-0.00814  -0.01375  -0.01036  -0.00782  -0.00014  0.00614  0.00113  -0.00744  201 

-0.00987  -0.00432  0.00332  -0.00196  -0.00883  -0.00171  0.00008  -0.00105  202 

-0.00458  0.00326  0.00637  0.00175  0.01363  0.01516  0.00405  0.00112  203 

-0.00042  -0.00296  -0.00334  0.00346  0.00734  0.00117  0.00074  0.00736  204 

0.01305  0.00937  0.00492  0.00476  0.00452  -0.00131  -0.00478  -0.00833  205 

-0.00594  -0.00493  -0.00572  -0.00769  -0.00519  -0.00077  0.00435  0.00944  206 

0.00797  0.00689  0.00597  -0.00162  -0.00001  -0.00187  -0.00538  -0.00212  207 

0.00196  -0.00355  -0.00234  -0.00053  0.00061  -0.00275  -0.00431  0.00156  208 

0.00099  -0.00087  -0.00534  -0.00042  0.00440  0.00125  -0.00179  -0.00029  209 

0.00420  0.00873  0.00232  0.00236  0.00402  0.00218  0.00652  0.00714  210 

0.00041  -0.00780  -0.00370  0.00135  0.00254  0.00296  0.00201  0.00315  211 

0.00806  0.00585  0.00240  0.00201  0.00141  -0.00149  -0.00462  -0.00621  212 

-0.00528  -0.00608  -0.00850  -0.00667  -0.00066  -0.00324  -0.00071  -0.00217  213 

-0.00459  -0.00280  -0.00246  -0.00295  0,00012  0.00214  -0.00317  -0.00417  214 

-0.00561  -0.00801  -0.00409  -0.00225  -0.00371  0.00263  -0.003C3  -0.00208  215 

0.00017  -0.00240  -0.00030  -0.00081  -0.00396  -0.00724  0.00198  0.00630  216 

0.00356  0.00399  0.00396  0.00366  0.00043  -0.00190  0.00343  0.00061  217 

-0.00385  -0.00298  0.00076  -0.00179  -0.00166  -0.00192  0.00027  0.00409  218 

0.00450  0.00087  -0.00585  -0.00549  -0.00501  -0.00232  0.00212  0.00004  219 

-0.00178  0.00006  0.00018  -0.00189  -0.00588  -0.01025  -0.00791  0.00046  220 

-0.00421  -0.00756  -0.00048  0.00860  0.00606  0.00799  0.00519  -0.00399  221 

-0.00573  -0.00082  0.00731  0.00202  0.00415  0.00491  0.00502  0.00119  222 

-0.00330  -0.00503  -0.00221  -0.00285  -0.00586  -0.00953  -0.01040  -0.00525  223 

-0.00008  0.00083  0.00940  0.01439  0.00796  0.00377  -0.00201  0.00464  224 

0.01212  0.01037  0.00155  -0.00968  -0.00331  -0.00181  0.00343  0.00781  225 

0.00150  -0.00103  0.00140  -0.00006  -0.00545  -0.00038  0.00170  0.00386  226 

0.00453  0.00381  0.00080  -0.00276  0.00166  0.00204  0.00140  0.00343  227 
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0.00560 

0.00546 

0.00180 

-0.00141 

-0.00027 

-0.00087 

-0.00379 

-0.01018 

228 

-0.00799 

-0.00204 

-0.00492 

0.00122 

-0.00428 

-0.00634 

-0.00138 

0.00291 

229 

0.00336 

-0.00305 

-0.01595 

-0.00003 

0.00940 

0.00570 

0.00315 

0.00755 

230 

-0.00227 

-0.00019 

0.00464 

0.00433 

-0.00018 

-0.00015 

0.00351 

0.00793 

231 

0.00340 

-0.00091 

-0.00134 

-0.00063 

0.00099 

0.00257 

0.00320 

0.00227 

232 

-0.00150 

-0.00044 

0.00353 

0.00802 

0.00292 

-0.00261 

-0.00374 

-0.00288 

233 

-0.00360 

0.00369 

0.00984 

0.00620 

0.00638 

0.00612 

0.00371 

-0.00288 

234 

-0.00680 

-0.00078 

0.00961 

-0.00069 

0.00312 

0.00151 

0.00078 

-0.00128 

235 

-0.00528 

-0.01092 

-0.00871 

0.00074 

-0.00467 

-0.00113 

0 . 00680 

-0.00210 

236 

-0.00959 

-0.00009 

0.00430 

0.00009 

0.00923 

0.01201 

0.00479 

-0.00273 

237 

-0.00122 

-0.00366 

-0.00120 

0.01083 

-0.00002 

-0.00782 

-0.00656 

-0.00167 

238 

0.00085 

-0.00290 

-0.00291 

-0.00542 

-0.00715 

-0.00277 

0.00182 

0.00738 

239 

0.00267 

-0.00792 

-0.00918 

-0.00088 

0.00111 

-0.00308 

-0.01068 

-0.00549 

240 

0.00238 

0.00541 

0.00617 

0.00160 

-0.00530 

-0.01071 

-0.00747 

-0.00244 

241 

0.00112 

-0.00105 

-0.00316 

-0.00082 

-0.00390 

-0.00567 

0.00074 

0.00271 

242 

-0.00343 

-0.00642 

-0.00055 

-0.00409 

-0.01229 

-0.00478 

0.00025 

0 . 00281 

243 

0.00196 

0.00996 

0.00687 

-0.00482 

0.00244 

0.00702 

-0 . 00025 

-0.00578 

244 

0.00554 

0.01694 

0.00300 

-0.00745 

-0.00379 

-0.00143 

-0.00443 

-0.00423 

245 

0.00206 

0.00680 

-0.00042 

0.00647 

0.00482 

-0.00147 

-0.00288 

-0.00179 

246 

0.00209 

0.01019 

0.00571 

0.00180 

0.00593 

0.01130 

0.00726 

0.00596 

247 

0.00827 

0.00250 

-0.00167 

0.00186 

-0.00006 

-0.00163 

-0.00471 

-0.00738 

248 

-0.00999 

-0.00573 

-0.00147 

0.00298 

0.01015 

0.00287 

-0.00733 

-0.00010 

249 

0.00366 

0.00740 

0.00463 

-0.00107 

-0.00869 

-0.00025 

0.00574 

0.00011 

250 

0.00095 

0.00234 

0.00194 

-0.00064 

-0.00351 

-0.00212 

0.00033 

0.00165 

251 

0.00600 

0.00984 

0.00584 

-0.00597 

0.00197 

0.00268 

-0.00498 

-0.00370 

252 

0.00307 

0.00438 

-0.00118 

-0.00108 

0.00030 

-0.00243 

-0.00631 

-0.00363 

253 

-0.00317 

-0.00467 

-0.00851 

-0.00520 

0.00064 

0.00181 

-0.00349 

-0.00195 

254 

-0.00074 

-0.00335 

-0.00189 

-0.00637 

-0.00593 

-0.00091 

0.00204 

-0.00005 

255 

-0.00517 

-0.00781 

-0.00562 

-0.00368 

-0.00413 

-0.00607 

-0.00206 

0.00118 

256 

3  OPTION  3:  ASSIGN  OBJECT  MOTION 

9  0 

4  OPTION  4:  OBTAIN  STRAIN -COMPATIBLE  PROPERTIES 

0  20  .05  .65 

5  OPTION  5:  COMPUTE  MOTION  IN  SPECIFIC  LAYERS 
12345678 
11111111 
00000000 

9  OPTION  9:  COMPUTE  RESPONSE  SPECTRA 

1  1 

5  0  111 

.020  .050  .070  .100  .120 

16  OPTION  16;  COMPUTE  STRESS/STRAIN  HISTORY 

5111  512  .OOOSand  and  gravel  layer 

0  END  OF  INPUT 
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INPUT  EARTHQUAKE  MOTION; 


TIME (sec) 

ACCELERATION  (j 

.000 

.00000 

.020 

- .00434 

.040 

.00860 

.060 

.00540 

.080 

- .00565 

.100 

-  .00944 

.120 

-.00369 

.140 

- .00669 

.160 

- .00336 

.180 

-.00111 

.200 

.00358 

.220 

.00303 

.240 

-.00323 

.260 

- .00907 

.280 

-.01522 

.300 

-  .01029 

.320 

-  .00706 

.340 

-  .00194 

.360 

.00135 

.380 

.00191 

.400 

.00743 

(intermediate  lines  not  shovm) 


40.480 
40.500 
40.520 
40.540 
40.560 
40.580 
40.600 
40.620 
40.640 
40.660 
40.680 
40.700 
40.720 
40.740 
40.760 
40.780 
40.800 
40.820 
40.840 
40.860 
40.880 
40.900 
40.920 
40 . 940 


-.00317 
-.00467 
- .00851 
- .00520 
.00064 
.00181 
-.00349 
- .00195 
-.00074 
-.00335 
- .00189 
-.00637 
- .00593 
- .00091 
.00204 
- .00005 
- .00517 
- .00781 
- .00562 
- .00368 
- .00413 
- .00607 
-.00206 
.00118 
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AMAX 


(TOP  OF  LAYER) 

MAXIMUM 

DEPTH 

ACCELERATION 

(ft) 

(g) 

.0. 

.257 

8.0. 

.222 

20.0, 

.146 

30.0, 

.143 

41.0, 

.115 

61.0, 

.115 

80.0, 

.094 

100.0, 

.082 

"STRESS 

1 

MID -DEPTH 

EFF.  SHEAR 

EFF.  SHEAR 

LAYER 

(ft) 

STRESS  (psf) 

STRAIN  (%) 

1, 

4.0, 

79., 

17E-01 

2, 

14.0, 

196.  , 

13E-I-00 

3. 

25.0, 

216.  , 

16E-01 

4, 

35.5, 

259.  , 

12E-01 

5, 

51.0, 

339.  , 

16E-01 

6. 

70.5, 

422.  , 

15E-01 

7, 

90.0, 

447.  , 

llE-01 

8, 

110.0, 

441.  , 

90E-02 

14 


ACCSPEC 


ABSOLUTE  ACCELERATION  SPECTRA  (g)  : 

period  .  DAMPING  RATIOS 


(sec) 

.02 

.05 

.07 

.10 

.12 

.001, 

.257, 

.257, 

.257, 

.257, 

.257, 

.100, 

.477, 

.415, 

.403, 

.389, 

.382, 

.150, 

.803, 

.526, 

.431, 

.408, 

.401, 

.200, 

.562, 

.493, 

.460, 

.422, 

.401, 

.250, 

.687, 

.54^, 

.531, 

.504, 

.482, 

.300, 

1.337, 

.94/, 

.798, 

.676, 

.612, 

.350, 

.797, 

.566, 

.540, 

.506. 

.476, 

.400, 

.763, 

.576, 

.515, 

.450, 

.423, 

.450, 

.928, 

.575, 

.499, 

.434, 

.406, 

.500, 

.895, 

.534, 

.455, 

.388, 

.354. 

.550, 

.797, 

.433, 

.347, 

.283, 

.261, 

.600, 

.510, 

.324, 

.281, 

.240, 

.221. 

.650, 

.412, 

.280, 

.237. 

.197, 

.183, 

.700, 

.276, 

.218, 

.206, 

.192, 

.184, 

.750, 

.328, 

.232, 

.200, 

.176, 

.165, 

.800, 

.182, 

.140, 

.136, 

.132, 

.130, 

.850, 

.167, 

.116, 

.106, 

.101, 

.102, 

.900, 

.127, 

.091, 

.087, 

.087, 

.089, 

.950, 

.150, 

.105, 

.090, 

.079, 

.078, 

1.000, 

.149, 

.103, 

.089, 

.078, 

.074, 

1.100, 

.087, 

.069, 

.067, 

.062, 

.059, 

1.200, 

.087, 

.071, 

.065, 

.060, 

.058. 

1.300, 

.060, 

.062, 

.060, 

.057, 

.055, 

1.400, 

.076, 

.063, 

.057, 

.052, 

.049, 

1.500, 

.066, 

.047. 

.043, 

.040, 

.039. 

1.600, 

.054, 

.037, 

.033, 

.031, 

.033, 

1.700, 

.043, 

.032, 

.030, 

.030. 

.032, 

1.800, 

.036, 

.032, 

.029, 

.030, 

.031, 

1.900, 

.033, 

.031, 

.030, 

.030, 

.031, 

2.000, 

.050, 

.037, 

.033, 

.030, 

.029, 

2.250, 

.039, 

.034, 

.032, 

.030, 

.029, 

2.500, 

.037, 

.033, 

.032, 

.030, 

.029. 

2.750, 

.041, 

.033, 

.029, 

.027, 

.026. 

3.000, 

.022, 

.022, 

.022, 

.021, 

.020, 

3.250, 

.026, 

.023, 

.021, 

.019, 

.018. 

3.500, 

.029, 

.023, 

.020, 

.017, 

.016, 

3.750, 

.024, 

.019, 

.017, 

.015, 

.015, 

4.000, 

.016, 

.015, 

.015, 

.014, 

.013. 
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OUTPUT 


EXAMPLE  PROBLEM  FOR  WESHAKE  USER'S  MANUAL 


MAX.  NUMBER  OF  TERMS  IN  FOURIER  TRANSFORM  -  4096 

NECESSARY  LENGTH  OF  BLANK  COMMON  X  -  25619 


*******  OPTION  2  ***  READ  SOIL  PROFILE 


MSOIL  -  0 

ML  -  9 

MWL  -  4 

WW  -  .0624 

IDNT  -  EIGHT  LAYERS  OVERLYING  ROCK 

NEW  SOIL  PROFILE  NO.  0  IDENTIFICATION  -  EIGHT  LAYERS  OVERLYING  ROCK 


SHEAR/K2  FACTOR  WAVE  VELOCITY  INPUT  BY  LAYER 


NUMBER  OF  LAYERS  9  DEPTH  TO  BEDROCK  120.00 

NUMBER  OF  FIRST  SUBMERGED  LAYER  4  DEPTH  TO  V-aTER  LEVEL  30.00 

UNIT  WEIGHT  OF  WATER  -  .0624  kef 
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SAND,  Average  (Seed  &  Idriss  1970) 
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Damping 
Esc . 

.  Small 

Vs 

Strain 

K2max 

Gmax 

. Inil 

Vs 

:ial  Est. 
G 

G/Gmax 

Layer 

(-) 

(fps) 

(-) 

(ksf) 

(fps) 

(ksf) 

(-) 

******************************************************************************* 

1 

.050 

448. 

43. 

750. 

449. 

750. 

1.00 

2 

.050 

430. 

25. 

575. 

430. 

575. 

1.00 

3 

.050 

745. 

58. 

1725. 

745. 

1725. 

1.00 

4 

.050 

920. 

79. 

2631. 

920. 

2631. 

1.00 

5 

.050 

1004. 

76. 

3920, 

1005. 

3920. 

1.00 

6 

.050 

910. 

76. 

3218. 

910. 

3218. 

1.00 

7 

.050 

1090. 

98. 

4616. 

1091. 

4616. 

1.00 

8 

.050 

1155. 

101. 

5183. 

1156. 

5183. 

1.00 

9 

.020 

8000. 

298415. 

8004. 

298415. 

1.00 

PERIOD 

-  .53 

FROM  AVERAGE 

SHEAR  VELOCITY  - 

905.  FT/SEC 

MAXIMUM  AMPLIFICATION  -  15.39 
FOR  FREQUENCY  -  2.19  C/SEC. 

PERIOD  -  .46  SEC. 


*******  OPTION  8  ***  READ  RELATION  BETWEEN  SOIL  PROPERTIES  AND  STRAIN 


*******  OPTION  1  ***  read  input  MOTION 


EARTHQUAKE  -  Alaska  EQ  (7/30/72)  Sitka  Record 

2048  ACCELERATION  VALUES  AT  TIME  INTERVAL  .0200 

THE  VALUES  ARE  LISTED  ROW  BY  ROW  AS  READ  FROM  CARDS 
TRAILING  ZEROS  ARE  ADDED  TO  GIVE  A  TOTAL  OF  4096  VALUES 

MAXIMUM  ACCELERATION  -  .09012 

AT  TIME  -  10.16  SEC 

THE  VALUES  WILL  BE  MULTIPLIED  BY  A  FACTOR  -  1.000 

TO  GIVE  NEW  MAXIMUM  ACCELERATION  -  .09012 

MEAN  SQUARE  FREQUENCY  -  5.65  C/SEC. 


MAX  ACCELERATION  -  .09015  FOR  FREQUENCIES  REMOVED  ABOVE  25.00 

C/SEC. 


*******  OPTION  3  ***  READ  WHERE  OBJECT  MOTION  IS  GIVEN 


OBJECT  MOTION  IN  LAYER  NUMBER  9  OUTCROPPING 


18 


*******  OPTION  4  ***  OBTAIN  STRAIN  COMPATIBLE  SOIL  PROPERTIES 


MAXIMUM  NUMBER  OF  ITERATIONS  -  20 
MAXIMUM  ERROR  IN  PERCENT  -  .05 
FACTOR  FOR  EFFECTIVE  STRAIN  IN  TIME  DOMAIN  -  .65 


EARTHQUAKE 
SOIL  PROFILE 


Alaska  EQ  (7/30/72)  Sitka  Record 
EIGHT  LAYERS  OVERLYING  ROCK 


19 


H 

CO 


U4 

U4 

Cd 

X 

H 

H-l 

3 

X 


A« 


Q 

W  h 
M  M 

O 


flO 

«0 


•H  Ok  CO 

tn  «n  flO 

r-  r-  «M 


r«»  rj  «n 

m  r-  kO 

«H  M 


Cd 

X  q 

M  Cd 
H  ^ 


O'<^*^CMO00v£>VD 
fOr^COCsJCNi-Hi-Hi— 1 


Ofomrofncor^ro 


(N 

v> 


«H  kO 

flO  m  r«> 
CM  <n  «H 


Ok  00  O 

cn  ko  «-! 

CM  CM  -9 


i 

b) 


•»  •M  «n 


CO 

CO 

o 

C^ 

00 

o\ 

Os 

Id 

rH 

CO 

rv 

Os 

oi 

(S4 

• 

• 

• 

• 

• 

• 

H 

CO 

•d’ 

iH 

00 

m 

so 

CO 

Ck 

o\ 

00 

VD 

CM 

fH 

CO 

<»• 

00 

ON 

cn  in  «o  00 


3 

Z 


«n  ^  00  «n 


Ct« 

b« 

b) 


X 

H 

a 


CO 

CO 

Cd 


oooinoinoo 


^^inmiHOOO 
iH  CM  fo  ir>  a\  r-J 


X 

u 

z 

B 

H 

O 

o 

CQ 

o 

b* 

H 

X 

l-t 

b9 

Sb 

•H 

X 

H 

Q 

z 

3 

z 

X 

o 

u 

U 

b« 

I-I 

hJ 

«n 

H 

< 

>2 

u 

Cd 

H 

u 

X 

g 

>« 

fM 

M 

H 

o  m 

lO  if) 


o 

o 


if)  fit  io  00 


O 

a 

u 

£ 

M 

H 

z 


v> 

u 

3 


u 

l-t 

X 

H 


u 

Ok 


Oi 

U 


oooooooo 


OOCMOf-IOONOO 
1-1  t-l  1-1  CM  l-t  CM  CM 


ro-tfirvir)CMvcr~oo 


ftcMm-^rinvorKoo 


T3 

O 


a 

a 

z 

c; 

M 

00 

VO 

00 

VO 

'S- 

o 

0) 

Cd 

9 

m 

m 

m 

m 

m 

m 

m 

m 

CM 

in 

Os 

CM 

m 

o 

CM 

m 

0^ 

X 

o 

o 

o 

o 

o 

o 

o 

o 

X  H 

rs 

eo 

VO 

o 

rv 

1— 1 

o 

X 

HH 

& 

H  Z 

fH 

VO 

CM 

CM 

*H 

CM 

CM 

CM 

CO 

o 

H 

CO  CJ) 

o 

o 

o 

o 

o 

o 

o 

o 

X 

o 

9 

X 

• 

■ 

• 

• 

• 

4J 

U 

><  X 

•rH 

X 

CO 

o 

CM  i-J 
r-»  QJ 
\  W 
O  > 
n  O 

r--  c/5 

oi 

U 

►t 


u 

to 


H 
X 
tn  o 
to  i-i 

i-t  w 


u 
■-) 
U  t-l 

%'s 

X  X 
O'  X 

gK 

ss 

b]  c/5 


110 


THE  CALCULATION  HAS  BEEN  CARRIED  OUT  IN  THE  TIME  DOMAIN  WITH  EFF.  STRAIN  -  .65*  MAX;  STRAIN 


06  H 

I  g 


O 

U)  bk 
t/>  (A 

9  ^ 


hi 

(A 


S  Z 

s  s 


a 

ui 

M 

a 


Q 

3 

u 


b« 

b) 


H 

Z 

u 

06 

Ok 


r«.  W 


00 

«-l  ov 
o  M 


^  «  <M  O 

«n  r<»  A 

ro  f-t  A 

^  C4  <M  M 


<M  A  O  tA  <A  lA  A 

m  lA  A»  <n  A  <n  tH 

e4  A  A  ^  A  A  A 

<S  A  «A  O  rk  A 

^  ^  ^  A  A  *A 

A  <M  A  «-«  «H  A  aH 


M  M 


e  A 

A  A 
O  O 


r»  ^ 
r»  A 
CN  A 
A  A 
O  O 


A  r«.  <S 


Cx] 

£  U 
M  b3 
H  W 


cn 

CO 

Cx] 


s^OO^CMCMOOOvO 

porocNiCMCMCsif-ii-^ 


CO 

CO 

<1- 

o 

Ov 

Ov 

m 

po 

u 

VO 

CM 

o 

m 

en 

o 

0^ 

(u 

• 

• 

• 

» 

* 

H 

CO 

o 

VO 

in 

00 

cn 

CM 

CO 

Cl4 

o 

eo 

o 

CM 

00 

00 

rM 

CM 

MJ* 

m 

p*«* 

■•  • 

•o 

u 

o 

z 

o 

M 

m 

in 

00 

00 

ON 

vd 

VO 

CM 

0) 

•cC 

VO 

'd' 

1— < 

m 

00 

oi 

H 

ov 

o 

CM 

-d 

vO 

00 

vO 

bi 

H 

z 

o 

cn 

CM 

CM 

CM 

to 

u 

CO 

u 

o 

rM 

o 

o 

O 

o 

o 

O 

jk 

o 

QC« 

• 

• 

» 

> 

• 

4J 

04 

iH 

CO 

C3 

z 

l-H 

>< 

CM 

hJ 

X 

0^ 

H 

o 

o 

o 

m 

o 

in 

o 

o 

\  U3 

Ck 

H 

• 

• 

• 

• 

* 

o 

> 

U 

£ 

-d 

in 

in 

o 

o 

o 

m 

o 

a 

CM 

«n 

m 

OV 

'v 

CO 

s-^ 

cd 

OOOOOOOO 


ooc^^0^-^00^00 

,-1  ^  CM  rH  CN  CS 


U 

O'  5 
u  J 

ni  H 
^  X 
W  CD 

(0  »-i 

!-•  U 
< 


fr.  H 
M  Ok 
O 


u 

?: 

H 


S 

>* 

3 


•H  M  m 


O 

a 

(d 

z: 

l-H 

H 

Z 


N  <o  CO 


VI 

bl 

3 

< 

> 


X 

H 


U 

CH 


5  S 


u 

>< 

:s 


ro.tfmir>{Mvor~oo 


c-tcsro-tfinvor^oo 


Ill 


NOTE:  ITERATIONS  3  THROUGH  8  NOT  SHOWN!! 


ITERATION  NUMBER  9 

THE  CALCULATION  HAS  BEEN  CARRIED  OUT  IN  THE  TIME  DOMAIN  WITH  EFF.  STRAIN  -  .65*  MAX.  STRAIN 
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ITERATION  NUMBER  10 
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ITERATION  NUMBER  11 
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*******  OPTION  9  ***  COMPUTE  RESPONSE  SPECTRUM 


COMPUTE  RESPONSE  SPECTRUM  IN  LAYER  1 

RESPONSE  SPECTRUM  ANALYSIS  FOR  LAYER  NUMBER  1 

CALCULATED  FOR  DAMPING  .020  .050  .070  .100 

TIMES  AT  WHICH  MAX.  SPECTRAL  VALUES  OCCUR 
TD  -  TIME  FOR  MAX.  RELATIVE  DISP. 

TV  -  TIME  FOR  MAX.  RELATIVE  VEL. 

TA  -  TIME  FOR  MAX.  ABSOLUTE  ACC. 

DAMPING  RATIO  -  .02 

PERIOD  TIMES  FOR  MAXIMA 


TD 

TV 

TA 

.00 

10.3400 

10.3000 

10.3400 

.10 

10.3400 

10.9600 

10.3400 

.15 

3 . 6400 

3.9000 

3.6400 

.20 

13.4200 

13.4600 

13.4200 

.25 

13.4400 

2.8200 

13.4400 

.30 

11.1600 

11.2400 

11.1600 

.35 

4.0800 

4.0000 

4.0800 

.40 

16.0600 

15.9600 

16.0600 

.45 

15.8600 

15.9800 

15.8600 

.50 

20.2600 

20.1400 

20.2600 

.55 

19.1400 

19.2600 

19.1400 

.60 

16.5200 

16.3800 

16.5200 

.65 

15.5200 

15.6600 

15.5200 

.70 

10.3000 

10.4200 

10.3000 

.75 

13.4000 

13.5400 

13.4000 

.80 

13.2000 

13.3200 

13.2000 

.85 

11.3800 

10.2800 

11.3800 

.90 

7.4400 

15.6600 

7.4200 

.95 

13.4200 

10.2800 

13.4200 

1.00 

9.2800 

10.9800 

9.2800 

1.10 

22.8600 

23.0400 

22.8600 

1.20 

9.0000 

9.1800 

9.0000 

1.30 

9.0200 

28.0200 

9.0200 

1.40 

9 . 0400 

13.5400 

9.0400 

1.50 

10.8000 

13.3200 

10.7800 

1.60 

22.7000 

10.5800 

22.7000 

1.70 

22.0600 

13.3200 

22.0400 

1.80 

5.0400 

13.3200 

5 . 0400 

1.90 

5.0600 

10.5800 

5.0600 

2.00 

13.4400 

13.3200 

13.4400 

2.25 

14.0800 

15.6600 

14.0600 

2.50 

14.0200 

13.3200 

14.0000 

2.75 

15.5800 

13.3200 

15.5600 

3.00 

15.9800 

15.4400 

15.9600 

3.25 

16.0800 

15.4400 

16.0600 

3.50 

16.5200 

15.4600 

16.4800 

3.75 

18.9000 

15.4600 

18.8800 

4.00 

23.3800 

2.2800 

23.3600 

.120 
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SPECTRAL  VALUES -- 

Alaska  EQ  (7/30/72)  Sitka  Record  EIGHT  LAYERS  OVERLYING  ROCK  DAMPING  RATIO 
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TIMES  AT  WHICH  MAX  SPECTRAL  VALUES  OCCUR 
TD  -  TIME  FOR  MAlv.  RELATIVE  DISP. 

TV  -  TIME  FOR  MAX.  RELATIVE  VEL. 

TA  -  TIME  FOR  MAX.  ABSOLUTE  ACC. 

DAMPING  RATIO  -  .05 


PERIOD 

TD 

TIMES  FOR  MAXIMA 
TV 

TA 

.00 

10.3400 

13.3200 

10.3400 

.10 

10.3400 

3.4600 

10 . 3400 

.15 

3 . 6400 

3 . 6000 

3 . 6400 

.20 

13.4200 

10.4400 

13.4200 

.25 

10.5200 

2.8200 

13.4400 

.30 

11.1600 

11.0800 

11.1600 

.35 

3.9200 

10.7000 

3.9000 

.40 

15.8600 

15.9600 

15.8600 

.45 

13.7800 

13.6600 

13.7600 

.50 

13.8200 

13.7000 

13.8200 

.55 

19.1200 

19 . 2400 

19.1000 

.60 

15.9600 

15.8400 

15.9600 

.65 

13.9400 

15.6400 

13.9400 

.70 

14.0000 

13.5200 

13.9800 

.75 

13.4000 

13,3000 

13.4000 

.80 

2.7600 

13.3200 

2.7400 

.85 

11.3600 

10.2800 

11.3600 

.90 

11.4000 

15.6600 

11.4000 

.95 

11.0800 

13.3200 

11.0600 

1.00 

9.2800 

13.3200 

9.2800 

1.10 

9.3200 

9.1600 

9.3200 

1.20 

9.0000 

9.1600 

8.9800 

1.30 

9.0200 

9.1800 

9.0000 

1.40 

9.0400 

10.2800 

9.0200 

1.50 

9 . 0800 

13.3200 

9.0600 

1.60 

10.1600 

13.3200 

10.1400 

1.70 

5 . 0400 

13.3200 

5.0000 

1.80 

5 . 0400 

13.3200 

5.0200 

1.90 

13.4400 

13.3200 

13.4000 

2.00 

13.4600 

13.3200 

13.4200 

2.25 

14.0400 

13.3200 

14 . 0000 

2.50 

14.0400 

13.3200 

14.0000 

2.75 

15.5600 

13.3200 

15.5200 

3.00 

14.4600 

13.3200 

14.3800 

3.25 

14.4800 

15.4400 

14.4400 

3.50 

16.5000 

15.4600 

16.4600 

3.75 

18.8800 

15.4600 

18.8400 

4.00 

18 . 9400 

15.4600 

18.8800 
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SPECTRAL  VALUES- - 

Alaska  EQ  (7/30/72)  Sitka  Record  EIGHT  LAYERS  OVERLYING  ROCK  DAMPING  RATIO 
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1  PLOT  OF  ACCELERATION  SPECTRA 

100  PER  CENT  CORRESPONDS  TO  1.3369 
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CURVE  1  2.00  %  DAMPING 
CURVE  2  5.00  %  DAMPING 
CURVE  3  7.00  %  DAMPING 
CURVE  4  10.00  %  DAMPING 
CURVE  5  12.00  %  DAMPING 


ABSISSA 

CURVE  1 

CURVE  2 

CURVE  3 

CURVE  4 

CURVE  5 

.001 

.257 

.257 

.257 

.257 

.257 

.100 

.477 

.415 

.403 

.389 

.382 

.150 

.803 

.526 

.431 

.408 

.401 

.200 

.562 

.493 

.460 

.422 

.401 

.250 

.687 

.543 

.531 

.504 

.482 

.300 

1.337 

.947 

.798 

.676 

.612 

.350 

.797 

.566 

.540 

.506 

.476 

.400 

.763 

.576 

.515 

.450 

.423 

.450 

.928 

.575 

.499 

.434 

.406 

.500 

.895 

.534 

.455 

.388 

.354 

.550 

.797 

.433 

.347 

.283 

.261 

.600 

.510 

.324 

.281 

.240 

.221 

.650 

.412 

.280 

.237 

.197 

.183 

.700 

.276 

.218 

.206 

.192 

.184 

.750 

.328 

.232 

.200 

.176 

.165 

.800 

.182 

.140 

.136 

.132 

.130 

.850 

.167 

.116 

.106 

.101 

.102 

.900 

.127 

.091 

.087 

.087 

.089 

.950 

.150 

.105 

.090 

.079 

.078 

1.000 

.149 

.103 

.089 

.078 

.074 

1.100 

.087 

.069 

.067 

.062 

.059 

1.200 

.087 

.071 

.065 

.060 

.058 

1.300 

.060 

.062 

.060 

.057 

.055 

1.400 

.076 

.063 

.057 

.052 

.049 

1.500 

.066 

.047 

.043 

.040 

.039 

1.600 

.054 

.037 

.033 

.031 

.033 

1.700 

.043 

.032 

.030 

.030 

.032 

1.800 

.036 

.032 

.029 

.030 

.031 

1.900 

.033 

.031 

.030 

.030 

.031 

2.000 

.050 

.037 

.033 

.030 

.029 

2.250 

.039 

.034 

.032 

.030 

.029 

2.500 

.037 

.033 

.032 

.030 

.029 

2.750 

.041 

.033 

.029 

.027 

.026 

3.000 

.022 

.022 

.022 

.021 

.020 

3.250 

.026 

.023 

.021 

.019 

.018 

3.500 

.029 

.023 

.020 

.017 

.016 

3.750 

.024 

.019 

.017 

.015 

.015 

4.000 

.016 

.015 

.015 

.014 

.013 
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*******  OPTION  16  ***  COMPUTE  STRESS/STRAIN  HISTORY 
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